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ON  THE  STABILITY   OF  THE  MOTION 


SATURN'S    RINGS. 


There  are  some  questions  in  Astronomy,  to  which  we  are  attracted  rather  on  account  of 
their  peculiarity,  as  the  possible  illustration  of  some  unknown  principle,  than  from  any  direct 
advantage  which  their  solution  would  afford  to  mankind.  The  theory  of  the  Moon's  in- 
equalities, though  in  its  first  stages  it  presents  theorems  interesting  to  all  students  of  mechanics, 
has  been  pursued  into  such  intricacies  of  calculation  as  can  be  followed  up  only  by  those  who 
make  the  improvement  of  the  Lunar  Tables  the  object  of  their  lives.  The  value  of  the 
labours  of  these  men  is  recognised  by  all  who  are  aware  of  the  importance  of  such  tables  in 
Practical  Astronomy  and  Navigation.  The  methods  by  which  the  results  are  obtained  are 
admitted  to  be  sound,  and  we  leave  to  professional  astronomers  the  labour  and  the  merit  of 
developing  them. 

The  questions  which  are  suggested  by  the  appearance  of  Saturn's  Rings  cannot,  in  the 
present  state  of  Astronomy,  call  forth  so  great  an  amount  of  labour  among  mathematicians. 
I  am  not  aware  that  any  practical  use  has  been  made  of  Saturn's  Rings,  either  in  Astronomy 
or  in  Navigation.  They  are  too  distant,  and  too  insignificant  in  mass,  to  produce  any 
appreciable  effect  on  the  motion  of  other  parts  of  the  Solar  system  ;  and  for  this  very  reason  it 
is  difficult  to  determine  those  elements  of  their  motion  which  we  obtain  so  accurately  in  the 
case  of  bodies  of  greater  mechanical  importance. 

But  when  we  contemplate  the  Rings  from  a  purely  scientific  point  of  view,  they  become 
the  most  remarkable  bodies  in  the  heavens,  except,  perhaps,  those  still  less  useful  bodies — the 
spiral  nebulae.  When  we  have  actually  seen  that  great  arch  swung  over  the  equator  of  the 
planet  without  any  visible  connexion,  we  cannot  bring  our  minds  to  rest.  We  cannot  simply 
admit  that  such  is  the  case,  and  describe  it  as  one  of  the  observed  facts  in  nature,  not  ad- 
mitting or  requiring  explanation.  We  must  either  explain  its  motion  on  the  principles  of 
mechanics,  or  admit  that,  in  the  Saturnian  realms,  there  can  be  motion  regulated  by  laws 
which  we  are  unable  to  explain. 
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The  arrangement  of  the  rings  is  represented  in  the  figure  (l)  on  a  scale  of  one  inch  to  a 
hundred  thousand  miles.  S  is  a  section  of  Saturn  through  his  equator,  A,  B  and  G  are  the 
three  rings.  A  and  D  have  been  known  for  200  years.  They  were  mistaken  by  Galileo  for 
protuberances  on  the  planet  itself,  or  perhaps  satellites.  Huyghens  discovered  that  what  he 
saw  was  a  thin  flat  ring  not  touching  the  planet,  and  Ball  discovered  the  division  between  A 
and  B.  Other  divisions  have  been  observed  splitting  these  again  into  concentric  rings,  but 
these  have  not  continued  visible,  the  only  well-established  division  being  one  in  the  middle  of 
A.  The  third  ring  C  was  first  detected  by  Mr  Bond,  at  Cambridge,  U.  S.  on  November  15, 
18.50;  Mr  Dawes,  not  aware  of  Mr  Bond's  discovery,  observed  it  on  November  29th,  and 
Mr  Lassel  a  few  days  later.  It  gives  little  light  compared  with  the  other  rings,  and  is  seen 
where  it  crosses  the  planet  as  an  obscure  belt,  but  it  is  so  transparent  that  the  limb  of  the 
planet  is  visible  through  it,  and  this  without  distortion,  showing  that  the  rays  of  light  have 
not  passed  through  a  transparent  substance,  but  between  the  scattered  particles  of  a  dis- 
continuous stream. 

It  is  difficult  to  estimate  the  thickness  of  the  system  ;  according  to  the  best  estimates  it  is 
not  more  than  100  miles,  the  diameter  of  A  being  176,418  miles;  so  that  on  the  scale  of  our 
figure  the  thickness  would  be  one  thousandth  of  an  inch. 

Such  is  the  scale  on  which  this  magnificent  system  of  concentric  rings  is  constructed ;  we 
have  next  to  account  for  their  continued  existence,  and  to  reconcile  it  with  the  known  laws  of 
motion  and  gravitation,  so  that  by  rejecting  every  hypothesis  which  leads  to  conclusions  at 
variance  with  the  facts,  we  may  learn  more  of  the  nature  of  these  distant  bodies  than  the 
telescope  can  yet  ascertain.  We  must  account  for  the  rings  remaining  suspended  above  the 
planet,  concentric  with  Saturn  and  in  his  equatoreal  plane ;  for  the  flattened  figure  of  the 
section  of  each  ring,  for  the  transparency  of  the  inner  ring,  and  for  the  gradual  approach  of 
the  inner  edge  of  the  ring  to  the  body  of  Saturn  as  deduced  from  all  the  recorded  observations 
by  M.  Otto  Struve  (Sur  les  dimensions  des  Anneaux  de  Saturne — Recueil  de  Memoires 
Astronomiques,  Poulkowa.  15  Nov.  1851).  For  an  account  of  the  general  appearance  of  the 
rings  as  seen  from  the  planet,  see  Lardner  on  the  Uranography  of  Saturn,  Mem.  of  the  Astrono- 
mical Society,  1853.  See  also  the  article  "Saturn"  in  NichoFs  Cyclopcedia  of  the  Physical 
Sciences. 

Our  curiosity  with  respect  to  these  questions  is  rather  stimulated  than  appeased  by  the 
investigations  of  Laplace.  That  great  mathematician,  though  occupied  with  many  questions 
which  more  imperiously  demanded  his  attention,  has  devoted  several  chapters  in  various  parts 
of  his  great  work,  to  points  connected  with  the  Saturnian  System. 

He  has  investigated  the  law  of  attraction  of  a  ring  of  small  section  on  a  point  very  near  it, 
(Mec.  Cel.  Liv.  in.  Chap,  vi.),  and  from  this  he  deduces  the  equation  from  which  the  ratio 
of  the  breadth  to  the  thickness  of  each  ring  is  to  be  found, 

_  JF   p_  _         X(\-l) 

6  ~  3a3  p  ~  (\+  1)  (3\2  +  1)  ' 

where  B  is    the   radius    of   Saturn,  and  p  his   density  ;   a    the  radius   of  the  ring,   and   p'  its 

density ;     and   X  the    ratio  of  the  breadth   of    the   ring  to  its  thickness.      The  equation   for 

determining  \  when  e  is  given  has  one  negative  root  which   must   be  rejected,  and   two  roots 
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which  are  positive  while  e  <  0-0543,  and  impossible  when  e  has  a  greater  value.  At  the 
critical  value  of  e,  X  =  2'5Q4  nearly. 

The  fact  that  X  is  impossible  when  e  is  above  this  value,  shows  that  the  ring  cannot  hold 
together  if  the  ratio  of  the  density  of  the  planet  to  that  of  the  ring  exceeds  a  certain  value. 
This  value  is  estimated  by  Laplace  at  1*3,  assuming  a  =  2i2. 

We  may  easily  follow  the  physical  interpretation  of  this  result,  if  we  observe  that  the 
forces  which  act  on  the  ring  may  be  reduced  to — 

(1)  The  attraction  of  Saturn,  varying  inversely  as  the  square  of  the  distance  from  his 
centre. 

(2)  The  centrifugal  force  of  the  particles  of  the  ring,  acting  outwards,  and  varying 
directly  as  the  distance  from  Saturn's  polar  axis. 

(3)  The  attraction  of  the  ring  itself,  depending  on  its  form  and  density,  and  directed, 
roughly  speaking,  towards  the  centre  of  its  section. 

The  first  of  these  forces  must  balance  the  second  somewhere  near  the  mean  distance  of  the 
ring.  Beyond  this  distance  their  resultant  will  be  outwards,  within  this  distance  it  will  act 
inwards. 

If  the  attraction  of  the  ring  itself  is  not  sufficient  to  balance  these  residual  forces,  the 
outer  and  inner  portions  of  the  ring  will  tend  to  separate,  and  the  ring  will  be  split  up;  and  it 
appears  from  Laplace's  result  that  this  will  be  the  case  if  the  density  of  the  ring  is  less  than 
JL2.  of  that  of  the  planet. 

This  condition  applies  to  all  rings  whether  broad  or  narrow,  of  which  the  parts  are 
separable,  and  of  which  the  outer  and  inner  parts  revolve  with  the  same  angular  velocity. 

Laplace  has  also  shown  (Liv.  v.  Chap,  in.),  that  on  account  of  the  oblateness  of  the  figure 
of  Saturn,  the  planes  of  the  rings  will  follow  that  of  Saturn's  equator  through  every  change 
of  its  position  due  to  the  disturbing  action  of  other  heavenly  bodies. 

Besides  this,  he  proves  most  distinctly  (Liv.  in.  Chap,  vi.),  that  a  solid  uniform  ring 
cannot  possibly  revolve  about  a  central  body  in  a  permanent  manner,  for  the  slightest  dis- 
placement of  the  centre  of  the  rjng  from  the  centre  of  the  planet  would  originate  a  motion 
which  would  never  be  checked,  and  would  inevitably  precipitate  the  ring  upon  the  planet,  not 
necessarily  by  breaking  the  ring,  but  by  the  inside  of  the  ring  falling  on  the  equator  of  the 
planet. 

He  therefore  infers  that  the  rings  are  irregular  solids,  whose  centres  of  gravity  do  not 
coincide  with  their  centres  of  figure.  We  may  draw  the  conclusion  more  formally  as  follows, 
"If  the  rings  were  solid  and  uniform,  their  motion  would  be  unstable,  and  they  would  be 
destroyed.  But  they  are  not  destroyed,  and  their  motion  is  stable;  therefore  they  are;  either 
not  uniform  or  not  solid." 

I    have    not    discovered*    either   in    the   works   of  Laplace    or   in    those   of   more  recent 


"  Since  this  was  written,  Prof.  Challis  has  pointed  out  to  me 

three  important  papers  in  Gould's  Astronomical  Journal: Mr 

G.  P.  Bond  on  the  Ringsof  Saturn  (May  1851)  and  Prof.  B.  Pierce 
of  Harvard  University  on  the  Constitution  of  Saturn's  Rings 


(June  1851),  and  on  the  Adams'  Prize  Problem  for  1856  (Sept. 
1855).  These  American  mathematicians  have  both  considered 
the  conditions  of  statical  equilibrium  of  a  transverse  section  of 
a  ring,  and  have  come  to  the  conclusion  that  the  rings,  if  they 
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mathematicians,  any  investigation  of  the  motion  of  a  ring  either  not  uniform  or  not  solid. 
So  that  in  the  present  state  of  mechanical  science,  we  do  not  know  whether  an  irregular  solid 
ring,  or  a  fluid  or  disconnected  ring  can  revolve  permanently  about  a  central  body  ;  and  the 
Saturnian  svstem  still  remains  an  unregarded  witness  in  heaven  to  some  necessary,  but  as  yet 
unknown,  development  of  the  laws  of  the  universe. 

We  know,  since  it  has  been  demonstrated  by  Laplace,  that  a  uniform  solid  ring  cannot 
revolve  permanently  about  a  planet.  We  propose  in  this  Essay  to  determine  the  amount  and 
nature  of  the  irregularity  which  would  be  required  to  make  a  permanent  rotation  possible. 
We  shall  find  that  the  stability  of  the  motion  of  the  ring  would  be  ensured  by  loading  the 
ring  at  one  point  with  a  heavy  satellite  about  4,1  times  the  weight  of  the  ring,  but  this  load, 
besides  being  inconsistent  with  the  observed  appearance  of  the  rings,  must  be  far  too 
artificially  adjusted  to  agree  with  the  natural  arrangements  observed  elsewhere,  for  a  very 
small  error  in  excess  or  defect  would  render  the  ring  again  unstable. 

We  are  therefore  constrained  to  abandon  the  theory  of  a  solid  ring,  and  to  consider  the 
case  of  a  ring,  the  parts  of  which  are  not  rigidly  connected,  as  in  the  case  of  a  ring  of 
independent  satellites,  or  a  fluid  ring. 

There  is  now  no  danger  of  the  whole  ring  or  any  part  of  it  being  precipitated  on  the  body 
of  the  planet.  Every  particle  of  the  ring  is  now  to  be  regarded  as  a  satellite  of  Saturn, 
disturbed  by  the  attraction  of  a  ring  of  satellites  at  the  same  mean  distance  from  the  planet, 
each  of  which  however  is  subject  to  slight  displacements.  The  mutual  action  of  the  parts  of 
the  rino-  will  be  so  small  compared  with  the  attraction  of  the  planet,  that  no  part  of  the  ring 
can  ever  cease  to  move  round  Saturn  as  a  satellite. 

But  the  question  now  before  us  is  altogether  different  from  that  relating  to  the  solid  ring. 
We  have  now  to  take  account  of  variations  in  the  form  and  arrangement  of  the  parts  of  the 
ring,  as  well  as  its  motion  as  a  whole,  and  we  have  as  yet  no  security  that  these  variations 
may  not  accumulate  till  the  ring  entirely  loses  its  original  form,  and  collapses  into  one  or 
more  satellites,  circulating  round  Saturn.  In  fact  such  a  result  is  one  of  the  leading  doctrines 
of  the  "nebular  theory"  of  the  formation  of  planetary  systems:  and  we  are  familiar  with  the 
actual  breaking  up  of  fluid  rings  under  the  action  of  "capillary"  force,  in  the  beautiful 
experiments  of  M.  Plateau. 

In  this  essay  I  have  shewn  that  such  a  destructive  tendency  actually  exists,  but  that  by 
the  revolution  of  the  ring  it  is  converted   into  the  condition  of  dynamical  stability.      As  the 


move  (ach  as  a  whole,  must  be  very  narrow  compared  with  the 
observed  rings,  so  that  in  reality  there  must  be  a  great  number 
of  them,  each  revolving  with  its  own  velocity.  They  have  also 
entered  on  the  question  of  the  fluidity  of  the  rings,  and  Prof. 
Pierce  has  made  an  investigation  as  to  the  permanence  of  the 
motion  of  an  irregular  solid  ring  and  of  a  fluid  ring.  The 
paper  in  which  these  questions  are  treated  at  large  has  not  (so 
far  as  I  am  aware)  been  published,  and  the  references  to  it  in 
Gould's  Journal  are  intended  to  give  rather  a  popular  account 
of  the  results,  than  an  accurate  outline  of  the  methods  employed. 
In  treating  of  the  attractions  of  an  irregular  ring,  he  makes 
admirable  use  of  the  theory  of  potentials,  but  his  published 


investigation  of  the  motion  of  such  a  body  contains  some  over- 
sights which  are  due  perhaps  rather  to  the  imperfections  of 
popular  language  than  to  any  thing  in  the  mathematical  theory. 
The  only  part  of  the  theory  of  a  fluid  ring  which  he  has  yet 
given  an  account  of,  is  that  in  which  he  considers  the  form  of 
the  ring  at  any  instant  as  an  ellipse ;  corresponding  to  the  case 
where  n  =  m,  and  m=\.  As  I  had  only  a  limited  time  for 
reading  these  papers,  and  as  I  could  not  ascertain  the  methods 
used  in  the  original  investigations,  I  am  unable  at  present  to 
state  how  far  the  results  of  this  essay  agree  with  or  differ  from 
those  obtained  by  Professor  Pierce. 
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scientific  interest  of  Saturn's  Rings  depends  at  present  mainly  on  this  question  of  their 
stability,  I  have  considered  their  motion  rather  as  an  illustration  of  general  principles,  than  as 
a  subject  for  elaborate  calculation,  and  therefore  I  have  confined  myself  to  those  parts  of  the 
subject  which  bear  upon  the  question  of  the  permanence  of  a  given  form  of  motion. 

There  is  a  very  general  and  very  important  problem  in  Dynamics,  the  solution  of  which 
would  contain  all  the  results  of  this  Essay  and  a  great  deal  more.      It  is  this — 

"Having  found  a  particular  solution  of  the  equations  of  motion  of  any  material  system, 
to  determine  whether  a  slight  disturbance  of  the  motion  indicated  by  the  solution  would  cause 
a  small  periodic  variation,  or  a  total  derangement  of  the  motion." 

The  question  may  be  made  to  depend  upon  the  conditions  of  a  maximum  or  a  minimum 
of  a  function  of  many  variables,  but  the  theory  of  the  tests  for  distinguishing  maxima  from 
minima  by  the  Calculus  of  Variations  becomes  so  intricate  when  applied  to  functions  of  several 
variables,  that  I  think  it  doubtful  whether  the  physical  or  the  abstract  problem  will  be  first 
solved. 


PART    I. 

ON  THE  MOTION  OF  A  RIGID  BODY  OF  ANY  FORM  ABOUT  A  SPHERE. 

We  confine  our  attention  for  the  present  to  the  motion  in  the  plane  of  reference,  as  the 
interest  of  our  problem  belongs  to  the  character  of  this  motion,  and  not  to  the  librations,  if 
any,  from  this  plane. 

Let  S  (Fig.  2)  be  the  centre  of  gravity  of  the  sphere,  which  we  may  call  Saturn,  and  R 
that  of  the  rigid  body,  which  we  may  call  the  Ring.      Join  RS,  and  divide  it  in  G  so  that 

SG  :   GR  ::   R  :   S, 
R  and  S  being  the  masses  of  the  Ring  and  of  Saturn  respectively. 

Then  G  will  be  centre  of  gravity  of  the   system,   and   its  position   will   be  unaffected  by 
any  mutual  action  between   the  parts  of  the   system.      Assume  G  as  the   point  to  which  the 
motions  of  the  system  are  to  be  referred.      Draw  GA  in  a  direction  fixed  in  space. 
Let  AGR  =  6,  and  SR  =  r,  then 

O  TO 

GR  = -r,  and  GS  = 


S+R  S  +  R 

so  that  the  positions  of  S  and  R  are  now  determined. 

Let  BRB'  be  a  straight  line  through  R,  fixed  with  respect  to  the  substance  of  the  ring, 
and  let  BRK  =  (p. 

This  determines  the  angular  position  of  the  ring,  so  that  from  the  values  of  r,  6,  and  cf> 
the  configuration  of  the  system  may  be  deduced,  as  far  as  relates  to  the  plane  of  reference. 

We  have  next  to  determine  the  forces  which  act  between  the  ring  and  the  sphere,  and  this 
we  shall  do  by  means  of  the  potential  function  due  to  the  ring,  which  we  shall  call  V. 

The  value  of  V  for  any  point  of  space  S,  depends  on  its  position  relatively   to  the  ring, 
and  it  is  found  from  the  equation 

■where  dm  is  an  element  of  the  mass  of  the  ring,  and  r  is  the  distance  of  that  element  from 
the  given  point,  and  the  summation  is  extended  over  every  element  of  mass  belonging  to  the 
ring.  V  will  then  depend  entirely  upon  the  position  of  the  point  S  relatively  to  the  ring, 
and  may  be  expressed  as  a  function  of  r,  the  distance  of  S  from  R,  the  centre  of  gravity  of 
the  ring,  and  (p,  the  angle  which  the  line  SR  makes  with  the  line  RB,  fixed  in  the  ring. 

A  particle  P,  placed   at  S,   will,   by   the  theory  of  potentials,  experience  a    moving  force 

P  —  in  the  direction  which  tends  to  increase  r,  and  P  -    —  in  a  tangential  direction,  tending 
dr  r  d(p 

to  increase  (p. 
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Now  we  know  that  the  attraction  of  a  sphere  is  the   same  as  that  of  a  particle  of  equal 
mass  placed  at  its  centre.      The  forces  acting  between  the  sphere  and  the  ring  are  therefore 

S tending  to  increase  r,  and  a  tangential  force  S  -  —  ,  applied  at  S  tending  to  increase  d>. 

dr  r  d(p  T 

In  estimating  the  effect  of  this  latter  force  on  the  ring,  we  must  resolve  it  into  a  tangential 

1    dV  .  dV 

force  S  -  —  acting  at  R,  and  a  couple  S  — —  tending  to  increase  d>. 
r  d<p  °  v         d<p  °  T 

We  are  now  able  to  form  the  equations  of  motion  for  the  planet  and  the  ring. 

For  the  planet  t 


dt 


^[(^LVd1\__R_sdV 

dt\\S  +  RI    dt]  S  +  R      d0 

*  t    Rr    N   _s_R^(ddV  _      dV 

For  the  centre  of  gravity  of  the  ring, 

d  |  /    Sr    \"-de\  S  dV 

dt\\s+RJ  ~di]  =  ~irrRsdip (3)' 

d2   (    Sr    \               Sr      (d9y       ^dV 
Rdt>KsTR)-RJ^R\d-t)    =Sdr~ (4)- 

For  the  rotation  of  the  ring  about  its  centre  of  gravity, 

*»£«>  + v-s% v< 

where  k  is  the  radius  of  gyration  of  the  ring  about  its  centre  of  gravity. 

Equations  (3)  and  (4)  are  necessarily  identical  with  (l)  and  (2),  and  show  that  the  orbit 
of  the  centre  of  gravity  of  the  ring  must  be  similar  to  that  of  the  Planet.  Equations  (l) 
and  (3)  are  equations  of  areas,  (2)  and  (4)  are  those  of  the  radius  vector. 

Equations  (3)  (4)  and  (5)   may  be  thus  written, 

„  f     dr  dQ        adrB]       t  dV 

„  f<fr          id0y\                       dV 
R\dF-ridt)\-(R+V-dr  =° ™> 

„,t(<P0      <f0\        adV 

These  are  the  necessary  and  sufficient  data  for  determining  the  motion  of  the  ring,  the 
initial  circumstances  being  given. 
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Prob.  I.      To  find  the  conditions  under  which  a  uniform  motion  of  the  ring  is  possible. 

By  a  uniform  motion  is  here  meant  a  motion  of  uniform  rotation,  during  which  the  position 
of  the  centre  of  the  Planet  with  respect  to  the  ring  does  not  change. 

In  this  case  r  and  <p  are  constant,  and  therefore  V  and  its  differential  coefficients  are  given. 
Equation  (7)  becomes, 

which  shows  that  the  angular  velocity  is  constant,  and  that 

Ld9\*  R  +  SdV       , 

U)  "-    mr  d7  =  a3'>sa* (9)- 

d~9 
Hence,  —  =  0,  and  therefore  by  equation  (8), 
at 

dV  ,     . 

^  =  ° (10)- 

Equations  (9)  and  (10)  are  the  conditions  under  which  the  uniform  motion  is  possible,  and 
if  they  were  exactly  fulfilled,  the  uniform  motion  would  go  on  for  ever  if  not  disturbed. 
But  it  does  not  follow  that  if  these  conditions  were  nearly  fulfilled,  or  that  if  when  accurately 
adjusted,  the  motion  were  slightly  disturbed,  the  motion  would  go  on  for  ever  nearly  uniform. 
The  effect  of  the  disturbance  might  be  either  to  produce  a  periodic  variation  in  the  elements 
of  the  motion,  the  amplitude  of  the  variation  being  small,  or  to  produce  a  displacement 
which  would  increase  indefinitely,  and  derange  the  system  altogether.  In  the  one  case  the 
motion  would  be  dynamically  stable,  and  in  the  other  it  would  be  dynamically  unstable. 
The  investigation  of  these  displacements  while  still  very  small  will  form  the  next  subject  of 
inquiry. 

Prob.  II.     To  find  the  equations  of  the  motion  when  slightly  disturbed. 

Let  r  =  r0,  9  =  wt  and  <p  =  <p0  in  the  case  of  uniform  motion,  and  let 

r  =  r0  +  rt, 

0  =  wt  +  et, 

when    the  motion  is  slightly  disturbed,  where  r  9   and  d)/  are  to  be  treated  as  small  quantities 
of  the  first   order,  and   their    powers    and   products  are  to   be   neglected.     We   may  expand 

dV  dV 

— —  and    — -  by  Taylor's  Theorem, 
dr  d<p    J        J 

dV      dV        d*V  d-V 

+    ~^r  f,  + 


dr        dr         dr      '       drd<p  '' 

dV      dV       d"V  d'V 

—  = 1 r  +    0  , 

d(p       d(p       drd(p    '        dip"    T' 

where  the  values  of  the  differential  coefficients   on   the   right-hand  side  of  the  equations  are 
those  in  which  r0  stands  for  r,  and  <pa  for  (p. 
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d'V      ,       d'V        „     d-T 
Calhng  —  =  Z,   -mj-M,    —  -  A, 

and  taking  account  of  equations  (9)  and  (10),  we  may  write  these  equations, 


dV 

Rr0       t 

Z> 

+  Mfa, 

dr 

tf  +  s"    ' 

dV 

dlp~ 

iHf^  +  Nfa . 

Substituting  these  values  in  equations  (6),  (7),  (8),  and  retaining  all  small  quantities  of  the 
first  order  while  omitting  their  powers  and  products,  we  have  the  following  system  of  linear 
equations  in  r,  dt  and  <plt 

drl 


I         dr,  „  ertf\ 

R  Vr°w  di  +  r°~  de)  +  (jR  +  S)  {Mri  +  N^  =  °     ■•••(»). 

„  (d?6,       dt"  \ 


Rk 


Prob.  III.      To  reduce  the  three  simultaneous  equations  of  motion  to  the  form  of  a  single 
linear  equation. 

d 
Let   us  write  n  instead  of  the  symbol  —  ,  then  arranging   the  equations  in  terms  of  r,   0, 

and  0, ,  they  may  be  written  : 

{2Rr0wn  +  (R  +  S)  M}  r,  +  (Rr0W)  0,  +  (#  +  S)  Nfa  =  0      .      .      .     (u), 
{Rn2- Rw--(R  +  S)L\rl-(2Rrtitvn)9l-(R  +  S)M(pl  =0      .      .      .      (15), 
-  {SM)  r,  +  {Rkhi2)  61  +  {Rk-n-  -  SN)  fa  =  0      .      .     .     (16). 
Here  we  have  three  equations  to  determine  three  quantities  rl   6J   fa ;  but   it  is  evident 
that  only  a  relation  can   be  determined  between    them,   and  that   in   the  process    for  finding 
their  absolute  values,  the  three  quantities  will  vanish  together,  and  leave  the  following  relation 
among   the   coefficients, 

-  \2Rr0wn  +  {R  +  S)  M}  {zRr0wn}  {Rtfrr  -  SN\ 
+  \Rn"  -  Rw2  ~{R  +  S)  L\  {.RAW}  \(R  +  S)  JST\ 

+  {SM)  (Rr*n)  {R  +  S)M-  {SM)  {<2Rr0wn)  {R  +  S)  NK  =  0    .     .     .     (17). 
+  \%Rr&m  +  {R  +  S) M\  \RkQn2}  \{R  +  S) M\ 

-  [Rn2-  R«o2  ~{R  +  S)\   \Rr2n2\  {Rk"n2  -  SN\ 

By   multiplying   up,   and  arranging    by   powers  of  n    and  dividing  by  Rn",  this  equation 
becomes 

An*+Bn"+C  =  0 (18), 

C 
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where 

A  =  R\*k\ 

B  =  SlTVoW  -  R  (R  +  S)  Lrjk-2  -R  \(R  +  S)  k2  +  Sr*}  N  I   .    .      (19). 

C  =  ^  | (i?  +  S)  k~  -  3Sr02\  w1  +  (R  +  S){(R  +  S)  k*  +  Sr0" \  (LN  -  A/2)) 

Here   we   have  a  biquadratic    equation  in  n    which  may  be  treated   as  a  quadratic  in  n", 

it  being  remembered   that  n  stands  for  the  operation  —  . 

Prob.  IV.     To  determine  whether  the  motion  of  the  ring  is  stable  or  unstable,  by  means 
of  the  relations  of  the  coefficients  A,  B,  C. 

The  equations  to  determine  the  forms  of  r,  0,  rf^  are  all  of  the  form 

.  d4u      _  d? v      _  .    % 

A  —  +  B-—  +  Cu  =  o (20), 

and  if  n  be  one  of  the  four  roots  of  equation  (18),  then 

u  =  Cenl 

will  be  one  of  the  four  terms  of  the  solution,  and  the  values  of  rx  Qx  and  <^,  will  differ  only  in 
the  values  of  the  coefficient  C. 

Let  us  inquire  into  the  nature  of  the  solution  in  different  cases. 

(1)  If  n  be  positive,  this  term  would  indicate  a  displacement  which   must  increase  inde- 
finitely, so  as  to  destroy  the  arrangement  of  the  system. 

(2)  If  n  be  negative,  the  disturbance  which  it  belongs  to  would  gradually  die  away. 

(3)  If  n  be  a  pure  impossible  quantity,  of  the  form  ±  a\/-  I,  then  there  will  be  a  term 

in  the  solution  of  the  form  C  cos  {at  +  a),  and  this  would  indicate  a  periodic  variation,  whose 

.         .     .  2tt 
amplitude  is  C.  and  period  —  . 
a 

(4)  If  n  be  of  the  form  b  ±  y/  -  la,  the  first  term   being  positive  and   the  second   im- 
possible, there  will  be  a  term  in   the  solution  of  the  form 

CeH  cos  (at  +  a), 

which  indicates  a  periodic  disturbance,  whose  amplitude  continually  increases  till  it  disarranges 
the  system. 

(5)  If  n  be  of  the  form  -  b  ±  \/~-~l  a,  a  negative   quantity  and  an   impossible  one,   the 
corresponding  term  of  the  solution  is 

Ce~ht  cos  (at  +  a), 
which  indicates  a  periodic  disturbance  whose  amplitude  is  constantly  diminishing. 
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It  is  manifest  that  the  first  and  fourth  cases  are  inconsistent  with  the  permanent  motion 
of  the  system.  Now  since  equation  (18)  contains  only  even  powers  of  n,  it  must  have 
pairs  of  equal  and  opposite  roots,  so  that  every  root  coming  under  the  second  or  fifth  cases, 
implies  the  existence  of  another  root  belonging  to  the  first  or  fourth.  If  such  a  root 
exists,  some  disturbance  may  occur  to  produce  the  kind  of  derangement  corresponding 
to  it,  so  that  the  system  is  not  safe  unless  roots  of  the  first  and  fourth  kinds  are  altogether 
excluded.  This  cannot  be  done  without  excluding  those  of  the  second  and  fifth  kinds,  so 
that,  to  insure  stability,  all  the  four  roots  must  be  of  the  third  kind,  that  is,  pure  impossible 
quantities. 

That  this  may  be  the  case,  both  values  of  n2  must  be  real  and  negative,  and  the  con- 
ditions of  this  are — 

1st.     That  A,  B,  and  C  should  be  of  the  same  sign, 

2ndly.     That  B2>  4, AC. 

When  these  conditions  are  fulfilled,  the  disturbances  will  be  periodic  and  consistent  with 
stability.  When  they  are  not  both  fulfilled,  a  small  disturbance  may  produce  total  derange- 
ment of  the  system. 

Frob.  V.  To  find  the  centre  of  gravity,  the  radius  of  gyration,  and  the  variations  of  the 
potential  near  the  centre  of  a  circular   ring  of  small  but  variable  section. 

Let  a  be  the  radius  of  the  ring,  and  let  9  be  the  angle  subtended  at  the  centre  between 
the  radius  through  the  centre  of  gravity  and  the  line  through  a  given  point  in  the  rino\ 
Then  if  fx  be  the  mass  of  unit  of  length  of  the  ring  near  the  given  point,  fj.  will  be  a 
periodic  function  of  9,  and  may  therefore  be  expanded  by  Fourier's  theorem  in  the  series, 

B  2  2 

M  = \  1  +  2/ cos  9  +  -  g  cos  29  +  -  h  sin  20  +  2  i  cos  (39  +  a)  +  Sec.  I  .     .      (21), 

2ttci  3  3 

where  /,  g,  h,  Sic.  are  arbitrary  coefficients,  and  R  is  the  mass  of  the  ring. 

(l)     The  moment  of  the  ring  about   the  diameter  perpendicular  to  the  prime  radius  is 


Rr0  =   /     yu<r  cos  9(19  =  Baf, 


therefore  the  distance  of  the  centre  of  gravity  from  the  centre  of  the  ring, 

r0  =  of. 

(2)  The  radius  of  gyration  of  the  ring  about  its  centre  in  its  own  plane  is  evidently 
the  radius  of  the  ring  =  a,  but  if  k  be  that   about  the  centre  of  gravity,  we  have 

jfes  +  r02  =  a2; 
.-.    F=a*(l-/2). 

(3)  The  potential  at  any  point  is  found  by  dividing  the   mass  of  each  element  by  its 
distance  from  the  given  point,  and  integrating  over  the  whole  mass. 
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Let  the  given  point  be  near  the  centre  of  the  ring,  and  let  its  position  be  defined  by 
the  co-ordinates  r  and  \J/,  of  which  r  is  small  compared  with  a. 

The  distance  (p)  between  this  point  and  a  point   in    the  ring  is 

i        1  r  1   fr'\ 2      3  /r'\ 2 

-=-  {l  +  -cos  ty-0)  +  -    -      +-(-)   cos2ty-0)  +  &c}. 

pa  a  4  \a  I        4  \a  I 

T 

The  other  terms  contain   powers  of  -  higher  than  the  second. 
We  have  now  to  determine  the  value  of  the  integral, 

V=   I     -add; 

Jo       P 

and  in  multiplying  the  terms  of  (ju)  by  those  of  I  -  ] ,  we  need  retain  only  those  which  contain 

constant  quantities,   for  all  those   which   contain  sines  or  cosines  of  multiples  of  (\|/  -  9)  will 
vanish  when  integrated  between  the  limits.      In  this  way  we  find 

R  r  1  r'2 

V  =  -  \l  +  f- cos\j,  + (1  +gcos  2\J,  +  h  sin  2  v//)  |        .      .      .      (22). 

a  ct  ^h  a 

The  other  terms  containing  higher  powers  of  -  . 

In  order  to  express   V  in  terms  of  r,  and  <£,,  as  we  have   assumed  in  the  former  investi- 
gation, we  must  put 

1 

r  cos  \|/-  =  —  7*!  +  -  ru<pf, 

r  sin  \|/  =  —  ra(pi, 


From  which   we  find 


vA=--f 

(CPV\  R 


d'V\       ..      R 


drd(p/0  2a8 


(24). 


These  results  may  be  confirmed  by  the  following  considerations  applicable  to  any  circular 
ring,  and  not  involving  any  expansion  or  integration.  Let  af  be  the  distance  of  the  centre  of 
gravity  from  the  centre  of  the  ring,  and  let  the  ring  revolve  about  its  centre  with  velocity  w. 
Then  the  force  necessary  to  keep  the  ring  in  that  orbit  will  be  -Rafw". 


OF   SATURN'S   RINGS.  13 


But  let  S  be  a  mass  fixed  at  the  centre  of  the  ring,  then  if 

a     s 

w  —  — , 
o3 

every  portion  of  the  ring  will  be  separately  retained  in  its  orbit  by  the  attraction  of  S, 
so  that  the  whole  ring  will  be  retained  in  its  orbit.  The  resultant  attraction  must  therefore 
pass  through  the  centre  of  gravity,  and  be 

f 

-  Rafw-  =  -  RS^; 

J  a2 

dV  f 

dr  a? 

dV      d*V      d2V 
I  he  equation  — -  +  — —  +  — —  +  iirp  =  0 

dx~       dy        dz' 

is  true  for  any  system  of  matter  attracting  according  to  the  law  of  gravitation.  If  we 
bear  in  mind  that  the  expression  is  identical  in  form  with  that  which  measures  the  total  efflux 

of  fluid  from  a  differential  element  of  volume,  when  — ■  —  —  are  the  rates  at  which  the  fluid 

dx   dy   dz 

passes  through  its  sides,  we  may  easily  form  the  equation  for  any  other  case.  Now  let  the 
position  of  a  point  in  space  be  determined  by  the  co-ordinates  r,  <p  and  z,  where  z  is  measured 
perpendicularly  to  the  plane  of  the  angle  <p.  Then  by  choosing  the  directions  of  the  axes 
■vyz,  so  as  to  coincide  with  those  of  the  radius  vector  r,  the  perpendicular  to  it  in  the  plane 
of  (p,  and  the  normal,  we  shall  have 

dx  =  dr        dz  =  rd<p  dy  =  dz 

dV_dV     dV  _i    dV     dV  _dV 
dx       dr       dy       r    d<p      dz       dz 

The  quantities  of  fluid  passing  through  an  element  of  area  in  each  direction  are 

dV    M  dV    l    J    A  dV     ,*A 

——  rdfbdz, drdz,     — —  rdmdr, 

dr       T  d(p  r  dz 

so  that  the  expression  for  the  whole  efflux  is 

l   dV      d2V      i    d*V     d°V  ,    . 

+  +_—  +  (25), 

r    dr        dr"       r2  d(p?       dz4 

which  is  necessarily  equivalent  to  the  former  expression. 

d'V 
Now  at  the  centre  of  the  ring  -—  may  be  found  by  considering  the  attraction  on  a  point 

just  above  the  centre  at  a  distance  z, 

dx  (a2  +  »2)* 

d'V         R     , 

= when  z  =  0. 

dz'  a6 
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1    dV  R 

Also  we  know  -  — —  = 5  and  r  =  of, 

r  dr  a 

.       d*V        1       d"V         R 

so  that  in  any  circular  ring    -—  +         „   — —  =  2  — (26), 

J  dr*      df-  dip4  or 

an  equation  satisfied  by  the  former  values  of  L  and  N. 

By  referring  to  the  original  expression  for  the  variable  section  of  the  ring,  it  appears  that 
the  effect  of  the  coefficient  f  is  to  make  the  ring  thicker  on  one  side  and  thinner  on  the  other 
in  a  uniformly  graduated  manner.  The  effect  of  g  is  to  thicken  the  ring  at  two  opposite  sides, 
and  diminish  its  section  in  the  parts  between.  The  coefficient  h  indicates  an  inequality  of  the 
same  kind,  only  not  symmetrically  disposed  about  the  diameter  through  the  centre  of  gravity. 

Other  terms  indicating  inequalities  recurring  three  or  more  times  in  the  circumference  of 
the  ring,  have  no  effect  on  the  values  of  L,  M  and  N.  There  is  one  remarkable  case,  however, 
in  which  the  irregularity  consists  of  a  single  heavy  particle  placed  at  a  point  on  the 
circumference  of  the  ring. 

Let  P  be  the  mass  of  the  particle,  and  Q  that  of  the  uniform  ring  on  which  it  is  fixed, 
then  R  =  P  +  Q, 


p 

Q 

+  2as* 

P  + 

Q  (        P\ 

R , 

+  «■); 

~  a3 

2a8 

•'•     g 

3P 
=  ~R~ 

<*r  ■ 

(27). 

Pbob.  VI.  To  determine  the  conditions  of  stability  of  the  motion  in  terms  of  the 
coefficients/,  g,  h,  which  indicate  the  distribution  of  mass  in  the  ring. 

The  quantities  which  enter  into  the  differential  equation  of  motion  (18)  are  R,  S,  k2,  r0,  o?, 
L,  M,  N.  We  must  observe  that  S  is  very  large  compared  with  R,  and  therefore  we  neglect 
R  in  those  terms  in  which  it  is  added  to  S,  and  we  put 

S  =  a  V 

R 
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Substituting  these  values  in  equation  (18)  and  dividing  by  R"a\f2,  we  obtain 

(1  -f)  n*  +  (l  -If*  +  \Pg)  nV  +   (~  -  &P  -  \  g2  -Iv  +  afg]  m*  =  0.      .      .    (28) 

The  condition  of  stability  is  that  this  equation  shall  give  both  values  of  n2  negative,  and  this 
renders  it  necessary  that  all  the  coefficients  should  have  the  same  sign,  and  that  the  square 
of  the  second  should  exceed  four  times  the  product  of  the  first  and  third. 

(1)  Now  if  we  suppose  the  ring  to  be  uniform,  /,  g  and  h  disappear,  and  the  equation 
becomes 

o 
«4+wV  +  -=0 (oq) 

which  gives  impossible  values  to  tr  and  indicates  the  instability  of  a  uniform  ring. 

(2)  If  we  make  g  and  h  =  0,  we  have  the  case  of  a  ring  thicker  at  one  side  than  the 
other,  and  varying  in  section  according  to  the  simple  law  of  sines.  We  must  remember, 
however,  that  /  must  be  less  than  1,  in  order  that  the  section  of  the  ring  at  the  thinnest 
part  may  be  real.      The  equation  becomes 

(1  -f)  n*  +(l-  5-p)  «V  +  (?  -  6f)  <•>*  =  0 (30). 

The  condition  that  the  third  term  should  be  positive  gives 

f*  <  -375. 
The  condition  that  n2  should  be  real  gives 

71  f  -  112/2  +  32  negative, 
which  requires  f*  to  be  between  "37445  and  1*8. 

The  condition  of  stability  is  therefore  that  f2  should  lie  between 
•37445  and  -375, 
but  the  construction   of  the   ring  on   this  principle  requires   that  f  should  be  less  than  -25, 
so  that  it  is  impossible  to  reconcile  this  form  of  the  ring  with  the  conditions  of  stability. 

(3)  Let  us  next  take  the  case  of  a  uniform  ring,  loaded  with  a  heavy  particle  at  a 
point  of  its  circumference.      We  have  then  g  =  Sf,  h  =  0,  and  the  equation  becomes 

Dividing  each  term  by  1  -/,  we  get 

0+/)»4+(l+/-^/!)«V+j{8(l+/)-8/2jft,«  =  0       .  _      (32) 

The  first  condition  gives  f  less  than  "8279. 

The  second  condition  gives  /  greater  than  -815865. 
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Let  us  assume  as  a  particular  case  between  these  limits  f=  "82,  which  makes  the  ratio 
of  the  mass  of  the  particle  to  that  of  the  ring  as  82  to  18,  then  the  equation  becomes 

1-82  n4  +  -8114  nV+  -9696  wl  =  0 (33), 

which  gives  \/-l»=±  "5916  w  or  ±  -3076*. 

These  values  of  n  indicate  variations  of  r,  0,  and  <pv  which  are  compounded  of  two  simple 
periodic  inequalities,  the  period  of  the  one  being  1-69  revolutions,  and  that  of  the  other  3-25l 
revolutions  of  the  ring.  The  relations  between  the  phases  and  amplitudes  of  these  inequalities 
must  be  deduced  from  equations  (14),  (15),  (16),  in  order  that  the  character  of  the  motion 
may  be  completely  determined. 

Equations  (14),  (15),  (16)  may  be  written  as  follows  : 

(4wa)  +  /io;2)  -  +  %fri!Ql  +/(3  -  g)  u?<f)l  =  0 (34), 

j«!Y'J(3+^-MY^,  =  0 (35), 

-fhiO*t±+2(l-f*)n%+{2(l-f*)n>-J*(S-g)<o}<j>1=0       .      (36). 

By  eliminating  one  of  the  variables  between  any  two  of  these  equations,  we  may  determine 
the  relation  between  the   two  remaining  variables.      Assuming  one  of  these  to  be  a  periodic 

function  of  t  of  the  form  A  cos  vt,  and  remembering  that  n  stands  for  the  operation  — ,  we  may 
find  the  form  of  the  other. 

Thus,  eliminating  0,  between  the  first  and  second  equations 

{n3  +  w<ua(5  -  g)  +  hw3\  -  +fwsl  (3  -g)w  --/mi^,  =  0.       .      .      .      (37). 

Assuming  —  =  A  sin  vt,  and  <p1=  Q  cos  (vt  -  fi) 

a 

{-}■*+  im  (5  -g)\  A  cos  vt  +  hw3  A  sin  vt+fw*  (3-g)Q  cos  (vt  -  /8)  +  -fhw*vQ  sin  (vt-fi). 

Equating  vt  to  0,  and  to  — ,  we  get  the  equations 

{v3-  iW  (5  -g)}  A  =fw2Q  \(3  -g)w  cos /3  -  -  Avsin/3>, 

-  hw3  A  =fw*Ql(3  -g)  00  sin/3  +  -  hv  cos  )3>  , 
from  which  to  determine  Q  and  /3. 
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In  all  cases  in  which  the  mass  is  disposed  symmetrically  about  the  diameter  through  the 
centre  of  gravity,  k  =  0  and  the  equations  may  be  greatly  simplified. 

Let  0!  =  Pcos  (vt  +  a),  then  the  second  equation  becomes 

|  v2  +  -  to2  (3  +  g)\  A  sin  t  —  2Pfwv  sin  (vt  +  a), 

whence  a  =  0     P  = ^— H —  A (38). 

2fwv  v     ' 

The  first  equation  becomes 

4Jft>i>cos  vt  -  2Pfv2  cos  vt  +  Qf  (3  -  <jr)  or  cos  (vt  -  /3)  =  0, 

whence  0-0     Q  =  »1~  h^  ^  -  9) 

f(3-g)  ws 

In  the  numerical  example  in  which  a  heavy  particle  was  fixed  to  the  circumference  of  the 
ring,  we  have,  when  f  =  *82, 

v  _  f'5916  P  _  13-21  Q       j  -  1-229 

w  ~  [-3076  A  ~  [572  A  =  \-    797  ' 

so  that  if  we  put  wt  =  90  =  the  mean  anomaly, 

-  =  A  sin  (-59l6  0o-  a)  +  B  sin  (-3076  9,-/3) (40), 

01  =  3-21  A  cos  (-59l6  0„-  a)  +  572  B  cos  (-3O76  0n-/3) (41), 

01=  -  1-229  A  cos  (-591600-  a)-  "797  B  cos  (-3O760O  -  /3)    .      .      .      .      (42). 

These  three  equations  serve  to  determine  rY  0,  and  <£s  when  the  original  motion  is  given. 
They  contain  four  arbitrary  constants  A,  B,  a,  /3-  Now  since  the  original  values  r,  0,  <plt  and 
also  their  first  differential  coefficients  with  respect  to  t  are  arbitrary,  it  would  appear  that  six 
arbitrary  constants  ought  to  enter  into  the  equation.  The  reason  why  they  do  not  is,  that  we 
assume  r0  and  0O  as  the  mean  valties  of  r  and  0  in  the  actual  motion.  These  quantities 
therefore  depend  on  the  original  circumstances,  and  the  two  additional  arbitrary  constants  enter 
into  the  values  of  r0  and  0O.  In  the  analytical  treatment  of  the  problem  the  differential 
equation  in  n  was  originally  of  the  sixth  degree  with  a  solution 

n*  =  0, 

which  implies  the  possibility  of  terms  in  the  solution  of  the  form 

Ct  +  D. 

The  existence  of  such  terms  depends  on  the  previous  equations,  and  we  find  that  a  term  of 
this  form  may  enter  into  the  value  of  0,  and  that  r,  may  contain  a  constant  term,  but  that  in 
both  cases  these  additions  will  be  absorbed  into  the  values  of  0O  and  r0. 

D 
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PART    II. 

ON  THE    MOTION    OF  A  RING,  THE    PARTS   OF    WHICH  ARE  NOT   RIGIDLY    CONNECTED. 


1.  In  the  case  of  the  Ring  of  invariable  form,  we  took  advantage  of  the  principle  that 
the  mutual  actions  of  the  parts  of  any  system  form  at  all  times  a  system  of  forces  in  equi- 
librium, and  we  took  no  account  of  the  attraction  between  one  part  of  the  ring  and  any  other 
part,  since  no  motion  could  result  from  this  kind  of  action.  But  when  we  regard  the  different 
parts  of  the  ring  as  capable  of  independent  motion,  we  must  take  account  of  the  attraction  on 
each  portion  of  the  ring  as  affected  by  the  irregularities  of  the  other  parts,  and  therefore  we 
must  begin  by  investigating  the  statical  part  of  the  problem  in  order  to  determine  the  forces 
that  act  on  any  portion  of  the  ring,  as  depending  on  the  instantaneous  condition  of  the  rest 
of  the  ring. 

In  order  to  bring  the  problem  within  the  reach  of  our  mathematical  methods,  we  limit  it 
to  the  case  in  which  the  ring  is  nearly  circular  and  uniform,  and  has  a  transverse  section 
very  small  compared  with  the  radius  of  the  ring.  By  analysing  the  difficulties  of  the  theory 
of  a  linear  ring,  we  shall  be  better  able  to  appreciate  those  which  occur  in  the  theory  of  the 
actual  rings. 

The  ring  which  we  consider  is  therefore  small  in  section,  and  very  nearly  circular  and 
uniform,  and  revolving  with  nearly  uniform  velocity.  The  variations  from  circular  form, 
uniform  section,  and  uniform  velocity  must  be  expressed  by  a  proper  notation. 

2.  To  express  the  position  of  an  element  of  a  variable  ring  at  a  given  time  in  terms  of 
the  original  position  of  the  element  in  the  ring. 

Let  «S"  (fig.  3)  be  the  central  body,  and  SA  a  direction  fixed  in  space. 

Let  SB  be  a  radius,  revolving  with  the  mean  angular  velocity  (to)  of  the  ring,  so  that 
ASB  =  wt. 

Let  ir  be  an  element  of  the  ring  in  its  actual  position,  and  let  P  be  the  position  it  would 
have  had  if  it  had  moved  uniformly  with  the  mean  velocity  w  and  had  not  been  displaced, 
then  BSP  is  a  constant  angle  =  s,  and  the  value  of  <s  enables  us  to  identify  any  element  of 
the  ring. 

The  element  may  be  removed  from  its  mean  position  P  in  three  different  ways. 

(1)  By  change  of  distance  from  S  by  a  quantity  pir  =  p. 

(2)  By  change  of  angular  position  through  a  space  Pp  =  a. 

(3)  By  displacement  perpendicular  to  the  plane  of  the  paper  by  a  quantity  £. 
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p,  a  and  £  are  all  functions  of  s  and  t.  If  we  could  calculate  the  attractions  on  any 
element  as  depending  on  the  form  of  these  functions,  we  might  determine  the  motion  of  the 
ring  for  any  given  original  disturbance.  We  cannot,  however,  make  any  calculations  of  this 
kind  without  knowing  the  form  of  the  functions,  and  therefore  we  must  adopt  the  following 
method  of  separating  the  original  disturbance  into  others  of  simpler  form,  first  given  in 
Fourier's  Traite  de  Chaleur. 

3.  Let  U  be  a  function  of  s,  it  is  required  to  express  U  in  a  series  of  sines  and  cosines 
of  multiples  of  s  between  the  values  s  =  0  and  s  =  2w. 

Assume  U  =  Ay  cos  s  +  A2  cos  2s  +  &c.  +  Am  cos  ms  +  Aa  cos  ns 

+  Bx  sin  s  +  B.j  cos  2s  +  &c.  +  Bm  sin  ms  +  Ba  sin  ws. 

Multiply  by  cos  msds  and  integrate,  then  all  terms  of  the  form 

/cos  ms  cos  nsds  and  J  cos  ms  sin  nsds 

will  vanish,  if  we  integrate  from  s  =  0  to  s  =  2ir,  and  there  remains 


/2n- 
U  cos  msds  —  irAm ,        /  U  sin  msds  =  ttB,„ 


If  we  can  determine  the  values  of  these  integrals  in  the  given  case,  we  can  find  the  proper 
coefficients  Am  Bm,  &c,  and  the  series  will  then  represent  the  values  of  U  from  s  =  0  to  s  =  2ir, 
whether  those  values  be  continuous  or  discontinuous,  and  when  none  of  those  values  are 
infinite  the  series  will  be  convergent. 

In  this  way  we  may  separate  the  most  complex  disturbances  of  a  ring  into  parts  whose 
form  is  that  of  a  circular  function  of  g  or  its  multiples.  Each  of  these  partial  disturbances 
may  be  investigated  separately,  and  its  effect  on  the  attractions  of  the  ring  ascertained  either 
accurately  or  approximately. 

4.  To  find  the  magnitude  and  direction  of  the  attraction  between  two  elements  of  a 
disturbed  ring. 

Let  P  and  Q  (fig.  4)  be  the  two  elements,  and  let  their  original  positions  be  denoted  by  st 
and  s2,  the  values  of  the  arcs  BP,  BQ  before  displacement.  The  displacement  consists  in  the 
angle  BSP  being  increased  by  cr,  and  BSQ  by  <r2,  while  the  distance  of  P  from  the  centre  is 
increased  by  p{  and  that  of  Q  by  p2.  We  have  to  determine  the  effect  of  these  displacements 
on  the  distance  PQ  and  the  angle  SPQ. 

Let  the  radius  of  the  ring  be  unity,  and  S2  -  Sx  =  29,  then  the  original  value  of  PQ  will 
be  2  sin  9,  and  the  increase  due  to  displacement 

=  G°2  +  p\)  sin  9  +  (o-j  -  a-,)  cos  9. 

We  may  write  the  complete  value  of  PQ  thus, 

PQ=  2  sin0  {l  +  -(/J2  +  /0,)  +  -  (<r2  -  a,)  cot  0}       ....     (I). 

d  2 
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The  original  value  of  the  angle  SPQ  was 0,  and  the  increase  due  to  displacement  is 

-(P2-P1)  cot0  --(ffss-  o"i). 
so  that  we  may  write  the  values  of  sin  SPQ  and  cos  SPQ, 

sin  SPQ  =  cos0  {l  +  -(p2  -pi)  -  -(<r,  -  ff,)  tan  9\     ....     (2), 

cos  SPQ  =  sin0  {l  --(p.,  -pt)cot29  +  -(a>-al)cot9\       .      .      (S). 

If  we  assume  the  masses  of  P  and  Q  each  equal  to  -  R,  where  R  is  the  mass  of  the  ring, 

and  /x  the  number  of  satellites  of  which    it  is  composed,    the  accelerating  effect  of  the  radial 
force  on  P  is 

1       cos  SPQ       1       R      ,  1  1 

^-^Q^  =  -^l^-^+^-i^-^^0--(<r2-ffl)cot0}      (4), 

and  the  tangential  force 

1       sin  SPQ      1  R  cos  9  ,        1  3  ,  1 

The  normal  force  is  -  R  — — =i. 
,u       8sin:i0 

5.      Let  us  substitute  for  p,  a  and  £  their  values  expressed  in  a  series  of  sines  and  cosines 
of  multiples  of  s,  the  terms  involving  ms  being 

px-  A  cos  (wis  +  a),  p2  =  A  cos  (ms  +  a  +  20), 

o-j  =  B  sin  (ms  +  /3),  <r2  =  B  sin  (ms  +  /3  +  20), 

£,  =  C  cos  (ws  +  7),  £2  =  C  cos  (ms  +  7  +  20). 

The  radial  force  now  becomes 

{  1  —  J  cos  (ms  +  a)  (l  +  cos  2»i0)  +  A  sin  (wis  +  a)  sin  2m0  1 

7—7;  ^     +  1 J  cos  (ms  +  a)  (1  -  cos  2»i0)  cot20  -  ^  sin  (ms  +  a)  sin  2m0  cot-  0  >    (6). 

<■     +  ^fi  sin  (ms  +  /3)  (1  -  cos  2m0)  cot  0  -  ^-B  cos  (wis  +  fi)  sin  2m9  cot  0  ' 

The  radial  component  of  the  attraction  of  a  corresponding  particle  on  the  other  side  of 

P  may  be  found  by  changing  the  sign  of  0.     Adding  the  two  together,  we  have  for  the  effect 

of  the  pair 

1        7? 

{l  —  A  cos  (ms  +  a)  (2  cos5  m0  —  sin2  m9  cotc  0) 


M  2  sin  0 


-  B  cos  (ms  +  (S)  -  sin  9.m9  cot  0}    .     .     .      .      (7). 
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Let  us  put 


/ 1  sin2  m9  cos2  9      cos2  m9  \ 


M 


-*C 


sin"0 

sin  2m9  cos  9\ 

4  sin2  9      j 


sin  9 


/sin2  m9  cos2  9       1  sin2  m.9 
N  =  2     r-r^ + 


2    sin  i 


^=2 


(8)*, 


sin30 
( sin'~  m9\ 
{2  sin3  flj 

U  sin  0/ 
where  the  summation  extends  to  all  the  satellites  on   the  same  side  of  P,  that  is,  every  value 


of  9  of  the  form  -  ir ,  where  a?  is  a  whole  number  less  than  — . 
fx  2 

The  radial  force  may  now  be  written 

P  =  -  R  {K  +  LA  cos  (ms  +  a)  -  MB  cos  (ms  +  /3)  \ 

M 

The  tangential  force  may  be  calculated  in  the  same  way,  it  is 


T  =  -  R  \MA  sin  {ms  +  a)  +  NB  sin  (ins  +  (i)  | 


•      •     (9). 


(10). 


The  normal  force  is 


Z  =  -  -  RJC  cos  (ms  +  y)        (II). 


6.      We   have  found   the   expressions  for   the   forces  which   act   upon   each   member  of  a 

system  of  equal  satellites  which  originally  formed   a   uniform   ring,  but  are  now  affected  with 

displacements   depending   on    circular  functions.      If  these  displacements  can   be  propagated 

m 
round  the  ring  in  the  form  of  waves  with  the  velocity  — ,  the  quantities  a,  /3  and  y  will  depend 

on  t,  and  the  complete  expressions  will  be 

p  =  A  cos  (ms  +  nt  +  o)  | 

a  =  B  sin  (ms  +  nt  +  /3)  > (12). 

£  =  Ccos  (ms  +  n't  +  y)  I 


.1/ 


•  The  following  values  of  several  quantities  which  enter  into 
these  investigations  are  calculated  for  a  ring  of  36  satellites. 
A' =  245 
sin2m6cos20         cos2m6 
sin3  &        '  "    sin  t) 
m=    0  0  43  -43        0  0 

m=    1  32  32  -16       16        37 

m=   2  107  28  26       25       115 

m=   3  212  25  81       28      221 

m=  4  401  24  177      32      411 

m=  9  975  20  468      30      986 

m  =  18  1569  18  767        0    1582 


When  ji  is  very  great, 

-      Z  =  .5259  when   hi  =  - , 

M  I  2 

=  .4342  m  =  T,, 
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Let  us  find  in  what  cases  expressions  such  as  these  will  be  true,  and  what  will  be  the 
result  when  they  are  not  true. 

Let  the  position  of  a  satellite  at  any  time  be  determined  by  the  values  of  r,  0  and  T, 
where  r  is  the  radius  vector  reduced  to  the  plane  of  reference,  0  the  angle  of  position  measured 
on  that  plane,  and  £  the  distance  from  it.     The  equations  of  motion  will  be 


m- 


d*r      „  l 
dt  r 


dr  dd>         d'd> 

— - — (-  r =  T 

dt  dt  dt" 


dX  X. 

dt2  r* 


(13). 


If  we  substitute  the  value  of  ^  in  the  third  equation   and  remember   that   r  is  nearly  =  1, 
we  find 


n2  =  S  +  -RJ 


(14). 


As  this  expression  is  necessarily  positive,  the  value  of  ri  is  always  real,  and  the  disturb- 
ances normal  to  the  plane  of  the  ring  can  always  be  propagated  as  waves,  and  therefore  can 
never  be  the  cause  of  instability.  We  therefore  confine  our  attention  to  the  motion  in  the 
plane  of  the  ring  as  deduced  from  the  two  former  equations. 

Putting  r  =  1  +  p  and  0  =  cot  +  s  +  a,  and  omitting  powers  and  products  of  p,  a  and  their 
differential  coefficients, 


da-  d2p 

r  dt  dt  r 

dp  d2<j      m 

dt  dt- 

Substituting  the  values  of  p  and  cr  as  given  above,  these  equations  become 

1                                    1      ^ 
to'  -  S RK  +  («r  +  25" RL  +  nq)  A  cos  (ms  +  nt  +  a) 

M  fx 

+  (2wn  +  -  RM)  B  cos  (ins  +  nt  +  fi)  =  0     .     . 

M 

(2ton  +  -  RM)  A  sin  (ms  +nt  +  a)  +  (n2  +  -  RN)  B  sin  (ms  +  nt  +  (Z)  =  0 


(15). 


06), 

07). 


Putting  for  (ms  +  nt)  any  two  different  values,  we  find  from  the  second  equation 

«  =  /? (IS), 

and  (own  +  I  RM)  A  +  (ri*  +  -RN)  B  =  0    . (19); 

and  from  the  first 


(to'  +  2S  -  -  RL  +  n)  A  +  (2w?i  +  -  RM)  B  =  0 
M  /x 


(20), 
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and  w*-  S  --RK  =  0 (21). 

M 

Eliminating  A  and  B  from  these  equations,  we  get 

n*  -  {»»s  -  2.9  +  -  R  (L  -  N)\  n'  -  \w  -RMn 

fX  fX 

+  («'  +  2S  -  -RL)-RN  -  —RM*  =  0  .     .     .     .     (22). 

MM  M' 

a  biquadratic  equation  to  determine  «. 

For  every  reaZ  value  of  n  there  are  terms  in  the  expressions  for  p  and  a  of  the  form 

A  cos  (ms  +  n£  +  a). 
For  every  pure  impossible  root  of  the  form  ±  v  —  In'  there  are  terms  of  the  forms 

^6*" 'cos  (wis  +  a). 

Although  the  negative  exponential  coefficient  indicates  a  continually  diminishing  displacement 
which  is  consistent  with  stability,  the  positive  value  which  necessarily  accompanies  it  indicates 
a  continually  increasing  disturbance,  which  would  completely  derange  the  system  in  course 
of  time. 

For  every  mixed  root  of  the  form  ±v/-l«'  +  »i|  there  are  terms  of  the  form 
Ae~nt  cos  {ms  +  nt  +  a). 

If  we  take  the  positive  exponential,  we  have  a  series  of  m  waves  travelling  with  velocity  — 

and  increasing  in  amplitude  with  the  coefficient  e+n(.  The  negative  exponential  gives  us  a 
series  of  m  waves  gradually  dying  away,  but  the  negative  exponential  cannot  exist  without  the 
possibility  of  the  positive  one  having  a  finite  coefficient,  so  that  it  is  necessary  for  the  stability 
of  the  motion  that  the  four  values  of  n  be  all  real,  and  none  of  them  either  impossible  quan- 
tities or  the  sums  of  possible  and  impossible  quantities. 

We  have  therefore  to  determine   the  relations  among  the  quantities  KLMNRS,  that  the 
equation 

w4  -  { S  +  - R  (3K  +  L  -  N)}  n2  -  4w-RMn 


+  {3S  +  -  R  {K  -  L)\-  RN  -  —R'M1  =  ZZ=  0 

IX  IX  [X~ 


may  have  four  real  roots. 


7.     In  the  first  place,  U  is  positive,  when  n  is  a  large  enough  quantity,  whether  positive 
or  negative. 

It  is  also  positive  when  n  =  0,  provided  S  be  large,   as  it   must   be,  compared  with  -  RL, 

IX 

-RM  and  -RN. 

M  M 
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If  we   can    now  find  a  positive  and   a  negative  value  of  n  for  which  U  is  negative,  there 
must  be  four  real  values  of  n  for  which  U  =  0,  and  the  four  roots  will  be  real. 


Now  if  we  put  n  =  ±  \/  -  v  S, 


U  =  -  -  S*  +  -  ~  R  (7N  ±  4  +/2  M  -  L  -  3K)  S  +  -R'  (KN  -  LN  -  AT), 

4  2  M  P- 

which  is  negative  if  S  be  large  compared  to  R. 

So  that  a  ring  of  satellites  can  always  be  rendered  stable  by  increasing  the  mass  of  the 
central  body  and  the  angular  velocity  of  the  ring. 

The  values  of  L,  M and  N  depend  on  m,  the  number  of  undulations  in   the  ring.     When 

m  =  — .  the  values  of  L  and  Ar  will  be  at  their   maximum   and  M  =  0.      If  we  determine  the 

2' 

relation  between  S  and  R  in  this  case   so  that   the  system   may  be  stable,  the  stability  of  the 
system  for  every  other  displacement  will  be  secured. 

8.  To  find  the  mass  which  must  be  given  to  the  central  body  in  order  that  a  ring  of 
satellites  may  permanently  revolve  round  it. 

We  have  seen  that  when  the  attraction  of  the  central  body  is  sufficiently  great  compared 
with  the  forces  arising  from  the  mutual  action  of  the  satellites,  a  permanent  ring  is  possible. 
Now  the  forces  between  the  satellites  depend  on  the  manner  in  which  the  displacement  of  each 
satellite  takes  place.  The  conception  of  a  perfectly  arbitrary  displacement  of  all  the  satellites 
may  be  rendered  manageable  by  separating  it  into  a  number  of  partial  displacements  depending 
on  periodic  functions.  The  motions  arising  from  these  small  displacements  will  all  take  place 
independently,  so  that  we  have  to  consider  only  one  at  a  time. 

Of  all  these  displacements,  that  which  produces   the  greatest  disturbing  forces  is   that  in 

which  consecutive  satellites  are  oppositely  displaced,  that  is,  when  m  =  — ,  for  then  the  nearest 

satellites  are  displaced  so  as  to  increase  as  much  as  possible  the  effects  of  the  displacement  of 
the  satellite  between  them.      If  we  make  fj.  a  large  quantity,  we  shall  have 

„  sin2»»6>  cos2  6       u3  .  ,  ,       „     .       (U1  .    n„   . 

2  -    -r-^ =  ~  (1  +  3"3  +  5"3  +  &C.)  =  -3  (V0518), 

Sin*  V  7T!  7T 

L  =  —„  -5259,  M  =  0,  N  =  2L,  K  very  small. 

w 

Let  -RL  =  x,  then  the  equation  of  motion  will  be 

W4  _  (#  -  x)  n2  +  2w  (35  -  x)  =  U  =  0 (23). 

The  conditions  of  this  equation  having  real  roots  are 

S>x (24), 

(5-  -  x)z  >  8x  (3S-x) (2.5). 
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The  last  condition  gives  the  equation 

S3  -  26Sx  +  9a>*  >  0, 

whence  S  >  25-649#,     or  S  <  0-351.1? (26). 

The  last  solution  is  inadmissible  because  iS  must  be  greater  than  sn,  so  that  the  true 
condition  is 

£  >  25-649.1?, 

>  25-649-  .R— -'5259, 

S  >  -4352,liIR (27). 

So  that  if  there  were  100  satellites  in  the  ring,  then 

8  >  4352  R 
is  the  condition  which  must  be  fulfilled  in  order  that  the  motion  arising  from  every  conceivable 
displacement  may  be  periodic. 

If  this  condition  be  not  fulfilled,  and  if  S  be  not  sufficient  to  render  the  motion  perfectly 
stable,  then  although  the  motion  depending  upon  long  undulations  may  remain  stable,  the 
short  undulations  will  increase  in  amplitude  till  some  of  the  neighbouring  satellites  are  brought 
into  collision. 

9.  To  determine  the  nature  of  the  motion  when  the  system  of  satellites  is  of  small  mass 
compared  with  the  central  body. 

The  equation  for  the  determination  of  n  is 

U  =  »4  -  {a,2  +  -  R(SK  +  L  -  N)\  n-  -  4-w-RMn 

+  {3w*--R(2K  +  L)\-RN-  -,RSJ\P  =  0     .      .     .      (28). 
M  M  /x 

When   R  is  very  small  we  may  approximate   to  the  values  of  n  by  assuming  that  two  of 
them  are  nearly  ±  w,  and  that  the  other  two  are  small. 
If  we  put  n  =  ±  to, 

U  =  -  -R  (2A'  +  L  ±  4,M  -  4JV)  w2  +  &c, 

M 

dU 

—  =  ±  2ft>3  +  &c. 
an 

Therefore  the  corrected  values  of  n  are 

<  1  1 

n  =  ±  \w  +  —  R  (2A'+  L  -  4JV")}  +  —  RM   .     .  .     .     (29). 

fXlti  '  fXW  v       ' 


The  small  values  of  n  are  nearly  i  V  3  -  RN,  correcting  them  in  the  same  way,  we  find 
the  approximate  values 

n  =  ±  sj 3-RN-2—RM (30). 

fl  fXW 

E 


26  ON   THE   STABILITY   OP  THE   MOTION 

The  four  values  of  n  are  therefore 


Wl  =  -  w R  (2JT  +L-*M  -  4JV) 

fxw 


-  V3  -  RN  -  —  .Kilt/ 


(31), 


rc3  =  +  Vs-  RN  -  —  RM 

nA  =  +  w  +  —  R  (2K  +  L  +  4M  -  iN) 

/Alt) 

and  the  complete  expression  for  p,  so  far  as  it  depends  on  terms  containing  ins,  is  therefore 
p  =  Ai  cos  (ms  +  nxt  +  at)  +  A2  cos  (ms  +  n2t  +  a2) 

+  A3  cos  (ms  +  n3t  +  a3)  +  Ai  cos  (ms  +  ntt  +  a4) (32), 

and  there  will  be  other  systems,  of  four  terms  each,  for  every  value  of  m  in  the  expansion  of 
the  original  disturbance. 

We    are  now   able   to   determine   the   value  of  a  from   equations   (18),    (20),  by  putting 
/3  =  a,  and 


B=  - 


2con  +  -  RM 
M_ 

ns  +  -  RN 


(33). 


So  that  for  every  term  of  p  of  the  form 

p  =  A  cos  (ms  +  nt  +  a) (34), 

there  is  a  corresponding  term  in  <r, 

2u>n  +  -  RM 

a  =  -  — ^^^^—  A  sin  (ms  +  nt  +  a) (35). 

n*  +-RN 

M 

10.  Let  us  now  fix  our  attention  on  the  motion  of  a  single  satellite,  and  determine  its 
motion  by  tracing  the  changes  of  p  and  <x  while  t  varies  and  s  is  constant,  and  equal  to  the 
value  of  s  corresponding  to  the  satellite  in  question. 

We  must  recollect  that  p  and  <t  are  measured  outwards  and  forwards  from  an  imaginary 
point  revolving  at  distance  1  and  velocity  w,  so  that  the  motions  we  consider  are  not  the 
absolute  motions  of  the  satellite,  but  its  motions  relative  to  a  point  fixed  in  a  i-evolving  plane. 
This  being  understood,  we  may  describe  the  motion  as  elliptic,  the  major  axis  being  in  the 

tangential  direction,  and  the  ratio  of  the  axes  being  nearly  2  - ,  which  is  nearly  2  for  nx  and  nt 

and  is  very  large  for  n-,  and  raa. 
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2tj- 
The  time  of  revolution  is    —  ,  or  if  we  take  a  revolution  of  the  ring  as  the  unit  of  time, 
n 

the  time  of  a  revolution  of  the  satellite  about  its  mean  position  is  — . 

n 

The  direction  of  revolution  of  the  satellite  about  its  mean  position  is  in  every  case  opposite 
to  that  of  the  motion  of  the  ring. 

11.  The  absolute  motion  of  a  satellite  may  be  found  from  its  motion  relative  to  the  ring 
by  writing 

r=l+p  =  l+^cos  (ms  +  nt  +  a), 

9  =  tot  +  s+a  =  wt+s  —  2  —  ^4  sin  (ms  +  nt  +  a), 
n 

When  n  is  nearly  equal  to  ±  w,  the  motion  of  each  satellite  in  space  is  nearly  elliptic. 
The  excentricity  is  A,  the  longitude  at  epoch  s,  and  the  longitude  when  at  the  greatest  dis- 
tance from  Saturn  is  for  the  negative  value  nx 

-  R  (2A'  +  L  -  4,M  -  4JV)  t  +  ms  +  a, 

JULO) 

and  for  the  positive  value  nt 

R  (fiK  +  L  +  4M  -  4>N)  t  -ms-  a. 

fXW 

We  must  recollect  that  in  all  cases  the  quantity  within  brackets  is  negative,  so  that  the 
major  axis  of  the  ellipse  travels  forwards  in  both  cases.  The  chief  difference  between  the  two 
cases  lies  in  the  arrangement  of  the  major  axes  of  the  ellipses  of  different  satellites.  In  the 
first  case  as  we  pass  from  one  satellite  to  the  next  in  front  the  axes  of  the  two  ellipses  lie  in 
the  same  order.  In  the  second  case  the  particle  in  front  has  its  major  axis  behind  that  of  the 
other.  In  the  cases  in  which  n  is  small  the  radius  vector  of  each  satellite  increases  and 
diminishes  during  a  periodic  time  of  several  revolutions.  This  gives  rise  to  an  inequality,  in 
which  the  tangential  displacement  far  exceeds  the  radial,  as  in  the  case  of  the  annual 
equation  of  the  Moon. 

12.  Let  us  next  examine  the  condition  of  the  ring  of  satellites  at  a  given  instant.  We 
must  therefore  fix  on  a  particular  value  of  t  and  trace  the  changes  of  p  and  cr  for  different 
values  of  s. 

From  the  expression  for  p  we  learn  that  the  satellites  form  a  wavy  line,  which  is  furthest 
from  the  centre  when  (ms  +  nt  +  a)  is  a  multiple  of  (2ir),  and  nearest  to  the  centre  for  inter- 
mediate values. 

From  the  expression  for  <j  we  learn  that  the  satellites  are  sometimes  in  advance  and  some- 
times in  the  rear  of  their  mean  position,  so  that  there  are  places  where  the  satellites  are  crowded 
together,  and  others  where  they  are  drawn  asunder.  When  n  is  positive,  B  is  of  the  opposite 
sign   to   J,  and  the  crowding  of  the  satellites   takes  place  when  they   are  furthest   from   the 

e  2 
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centre.  When  n  is  negative,  the  satellites  are  separated  most  when  furthest  from  the  centre, 
and  crowded  together  when  they  approach  it. 

The  form  of  the  ring  at  any  instant  is  therefore  that  of  a  string  of  beads  forming  a 
re-entering  curve,  nearly  circular,  but  with  a  small  variation  of  distance  from  the  centre 
recurring  to  times,  and  forming  m  regular  waves  of  transverse  displacement  at  equal  intervals 
round  the  circle.  Besides  these,  there  are  waves  of  condensation  and  rarefaction,  the  effect  of 
longitudinal  displacement.  When  n  is  positive  the  points  of  greatest  distance  from  the  centre 
are  points  of  greatest  condensation,  and  when  n  is  negative  they  are  points  of  greatest  rare- 
faction. 

13.  We  have  next  to  determine  the  velocity  with  which  these  waves  of  disturbance  are 
propagated  round  the  ring.  We  fixed  our  attention  on  a  particular  satellite  by  making  s 
constant,  and  on  a  particular  instant  by  making  t  constant,  and  thus  we  determined  the 
motion  of  a  satellite  and  the  form  of  the  ring.  We  must  now  fix  our  attention  on  a  phase  of 
the  motion,  and  this  we  do  by  making  p  or  a  constant.      This  implies 

ms  +  nt  +  a  =  constant, 

ds  n 

dt  m 

So  that  the  particular  phase  of  the  disturbance  travels  round  the  ring  with  an  angular  velocity 

= relative  to  the  ring  itself.      Now  the  ring  is  revolving  in  space  with  the  velocity  w,  so 

m 

that  the  angular  velocity  of  the  wave  in  space  is 

■&  =  <»-- (36). 

TO 

Thus  each  satellite  moves  in  an  ellipse,  while  the  general  aspect  of  the  ring  is  that  of  a 
curve  of  to  waves  revolving  with  velocity  •&.  This,  however,  is  only  the  part  of  the  whole 
motion,  which  depends  on  a  single  term  of  the  solution.  In  order  to  understand  the  general 
solution  we  must  shew  how  to  determine  the  whole  motion  from  the  state  of  the  ring  at  a 
given  instant. 

14.  Given  the  position  and  motion  of  every  satellite  at  any  one  time,  to  calculate 
the  position  and  motion  of  every  satellite  at  any  other  time,  provided  that  the  condition  of 
stability  is  fulfilled. 

The  position  of  any  satellite  may  be  denoted  by  the  values  of  p  and  a  for  that  satellite, 

and  its  velocity  and  direction  of  motion  are  then  indicated  by  the  values  of  -/-  and  — -  at  the 

J  J  dt  dt 

given  instant. 

These  four  quantities  may  have  for  each  satellite  any  four  arbitrary  values,  as  the  position 
and  motion  of  each  satellite  are  independent  of  the  rest,  at  the  beginning  of  the  motion. 

Each  of  these  quantities  is  therefore  a  perfectly  arbitrary  function  of  s,  the  mean  angular 
position  of  the  satellite  in  the  ring. 
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But  any  function  of  s  from  s  =  0  to  s  =  2tt,  however  arbitrary  or  discontinuous,  can  be 
expanded  in  a  series  of  terms  of  the  form  A  cos  (s  +  a)  +  A  cos  (2s  +  a')  +  &c.      See  §  3. 

dp         da 
Let  each  of  the  four  quantities  p,  -—  ,  a,  —r-  be  expressed  in  terms  of  such  a  series,  and  let 

at  at 

the  terms  in  each  involving  ms  be 

p  =  E  cos  (ms  +  e) (37), 

dp 

j-  =  F  cos  (ms  +/) (38), 

<r  =  G  cos  (ms  +  g) (39), 

da 

—  =  H  cos  (ms  +  h) (40). 

These  are  the  parts  of  the  values  of  each  of  the  four  quantities  which  are  capable  of  being 
expressed  in  the  form  of  periodic  functions  of  ms.  It  is  evident  that  the  eight  quantities 
E,  F,  G,  H,  e,  f,  g,  h,  are  all   independent  and  arbitrary. 

The  next  operation  is  to  find  the  values  of  L,  M,  N,  belonging  to  disturbances  in  the  ring 
whose  index  is  m,  [see  equation  (8)],  to  introduce  these  values  into  equation  (28),  and  to 
determine  the  four  values  of  n,  (w,  n2  n3  n4). 

This  being  done,  the  expression  for  p  is  that  given  in  equation  (32),  which  contains  eight 
arbitrary  quantities  (At  A2  A.t  At  a,  a2  a3  a4). 

Giving  t  its  original  value  in  this  expression,  and  equating  it  to  E  cos  (ms  +  e),  we  get  an 
equation  which  is  equivalent  to  two.      For,  putting  ms  =  0,  we  have 

At  coscti  +  -^2  cos  a2  +  A3cosas  +  ^4cosa4  =  E  cose (*!)• 

And  putting  ms  =  — ,  we  have  another  equation 

At  sin  a,  +  A*  sin  a2  +  At  sin  a3  +  -44  sin  a4  =  E  sin  e (42). 

Differentiating  (32)  with  respect  to  t,  we  get  two  other  equations 

-  Ax  nx  sin  a  —  &c.  =  F  cos/ (4S), 

Al  n^  cos  a  +  &c.  =  F  sin/ (44)- 

Bearing  in  mind  that  Bt  B2  &c.  are  connected  with  At  A2  &c.  by  equation  (33),  and  that 
B  is  therefore  proportional  to  A,  we  may  write  B  =  Afi,  where 

2wTl  H RM 

0= £ — . 

n"  +  -NR 
/3  being  thus  a  function  of  n  and  a  known  quantity. 
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The  value  of  a  then  becomes  at  the  epoch 

a  =  Ax  /3j  sin  (ms  +  at)  +  &c.  =  G  cos  {ms  +  g), 
from  which  we  obtain  the  two  equations 

Al(i1smal  +  &c.  =  G  cosg (45), 

A1  /3j  cos  a,  +  &c.  =  -  G  sin  g (46). 

Differentiating  with  respect  to  t,  we  get  the  remaining  equations 

A,  fit  nx  cosai  +  &c.  =  H  cos/t (47)> 

Al  filn1sin  at  +  &c.  =  i7sin  h (48). 

We  have  thus  found  eight  equations  to  determine  the  eight  quantities  Al  &c.  and  c^  &c: 
To  solve  them,  we  may  take  the  four  in  which  Ax  cos  e^  &c.  occur,  and  treat  them  as  simple 
equations,  so  as  to  find  Ax  cos  ax  &c.  Then  taking  those  in  which  Ax  sin  ax  &c.  occur,  and 
determining  the  values  of  those  quantities,  we  can  easily  deduce  the  values  of  Al  and  of  a, 
&c,  from  these. 

We  now  know  the  amplitude  and  phase  of  each  of  the  four  waves  whose  index  is  m. 
All  other  systems  of  waves  belonging  to  any  other  index  must  be  treated  in  the  same  way,  and 
since  the  original  disturbance,  however  arbitrary,  can  be  broken  up  into  periodic  functions  of 
the  form  of  equations  (37 — 40),  our  solution  is  perfectly  general,  and  applicable  to  every 
possible  disturbance  of  a  ring  fulfilling  the  condition  of  stability  (27). 

15.  We  come  next  to  consider  the  effect  of  an  external  disturbing  force,  due  either  to 
the  irregularities  of  the  planet,  the  attraction  of  satellites,  or  the  motion  of  waves  in  other  rings. 

All  disturbing  forces  of  this  kind  may  be  expressed  in  series  of  which  the  general  term  is 

A  cos  (vt  +  ms+a), 
where  v  is  an  angular  velocity  and  m  a  whole  number. 

Let  Pcos  (ms  +  vt+p)  be  the  central  part  of  the  force,  acting  inwards,  and  Q  sin  (ms  +  vt+q) 
the  tangential  part,  acting  forwards.  Let  p  =  A  cos  (ms  +  vt  +  a)  and  <r  =  B  sin  (ms  +  vt  +  /3), 
be  the  terms  of  p  and  cr  which  depend  on  the  external  disturbing  force.  These  will  simply 
be  added  to  the  terms  depending  on  the  original  disturbance  which  we  have  already  investigated, 
so  that  the  complete  expressions  for  p  and  a  will  be  as  general  as  before.  In  consequence  of 
the  additional  forces  and  displacements,  we  must  add  to  equations  (16)  and  (17),  respectively, 
the  following  terms  : 

\3w* R(ZK+L)  +  v*\  A  cos  (ms  +  vt+a) 

+  ('2(ou  +  -RM)  B  cos  (ms  +  vt  +  /3)-P  cos  (ms  +  vt  +p)  =  0   .      .     .     (49). 

M 

(Qoov  +  -  RM)  A  sin  (ms  +  vt  +  a) 
M 

+  (v2  +  -  RN)  B  sin  (ms  +  vt  +  (3)  +  Q  sin  (ms  +  vt  +  q)  =  0  .      .      .      (50). 
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Making  (ms  +  vt)  =  0  in  the  first  equation  and  -  in  the  second, 

\3m2-  -B{2K+L)  +  v2}  Acosa  +  {2u>v  +  - BM)  B  cos  /3  -  P  cos  p  =  0  .     .      .      (51). 
M  M 

(2wv  +  -  BM)  A  cos  a  +  (v*  +  -  BN)  B  cos  /3  +  Q  cos  q  =  0  .      .      .      (52). 
M  M 

Then  if  we  put 

U'  =  v'  -  W  +  -  B  (2A'  +  L  -  N)  }v*  -  4  -  i?J/v 

+  53o)2--JB(2/r+JL)}-i?^-  -^  B°M2.     .     .     .      (53), 
|K  MM 

we  shall  find  the  value  of  ^  cos  a  and  5  cos  /3  ; 

v-  +  -  BN                    2u>v  +  -BM 
A  cos  a  = jt, P  cos  p  +  p Q  cos  9         .      .      .      (54). 

Zwv  +  -  BM                   vs  +  3a)2  -  -B  (K  +  L) 
£008/3= =a Pcosp -^ Qcos9     .      .      (55). 

Substituting  sines  for  cosines  in  equations  (5l)  (52),  we  may  find  the  values  of  A  sin  a 
and  B  sin  /3. 

Now  U'  is  precisely  the  same  function  of  v  that  U  is  of  n,  so  that  if  v  coincides  with  one 
of  the  four  values  of  n,  U'  will  vanish,  the  coefficients  A  and  B  will  become  infinite,  and  the 
ring  will  be  destroyed.       The  disturbing  force  is  supposed  to  arise  from  a  revolving  body,  or 

an   undulation  of  any   kind  which  has  an  angular  velocity relatively  to   the  ring,   and 

'  m 

v 

therefore  an  absolute  angular  velocity  =  u> . 

m 

If  then  the  absolute  angular  velocity  of  the  disturbing  body  is  exactly  or  nearly  equal  to 
the  absolute  angular  velocity  of  any  of  the  free  waves  of  the  ring,  that  wave  will  increase  till 
the  ring  be  destroyed. 

The  velocities  of  the  free  waves  are  nearly 

«(l  +-),  w+  -  V  3  -BN,  co--  \/3-BN,  and  w{\  - -)       .     .      .      (56). 
\       m)  m  n  m  n  \       ml 

When  the  angular  velocity  of  the  disturbing  body  is  greater  than  that  of  the  first  wave, 
between  those  of  the  second  and  third,  or  less  than  that  of  the  fourth,  U'  is  positive.  When 
it  is  between  the  first  and  second,  or  between  the  third  and  fourth,  U'  is  negative. 
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Let  us  now  simplify  our  conception  of  the  disturbance  by  attending  to  the  central  force 
only,  and  let  us  put  p  =  0,  so  that  P  is  a  maximum  when  ms  +  vt  is  a  multiple  of  Zir.  We 
find  in  this  case  a  =  0,  and  (i  =  0.      Also 

v2  +  -  RN 
A= %r—  P (57), 

2uv  +  -  RM 
B= £ P (58). 

When  U  is  positive,  A  will  be  of  the  same  sign  as  P,  that  is,  the  parts  of  the  ring  will 
be  furthest  from  the  centre  where  the  disturbing  force  towards  the  centre  is  greatest.  When 
U  is  negative,  the  contrary  will  be  the  case. 

When  v  is  positive,  B  will  be  of  the  opposite  sign  to  A,  and  the  parts  of  the  ring  furthest 
from  the  centre  will  be  most  crowded.      When  v  is  negative  the  contrary  will  be  the  case. 

Let  us  now  attend  only  to  the  tangential  force,  and  let  us  put  q  =  0.  We  find  in  this  case 
also  a  =  0,  (Z  =  0, 

2o>«  +  -  RM 
*=  JJT Q (59), 

v°-  +  3u,2--R(K+L) 
B= ^ (60). 

The  tangential  displacement  is  here  in  the  same  or  in  the  opposite  direction  to  the 
tangential  force,  according  as  U'  is  negative  or  positive'.  The  crowding  of  satellites  is  at  the 
points  farthest  from  or  nearest  to  Saturn  according  as  v  is  positive  or  negative. 

16.  The  effect  of  any  disturbing  force  is  to  be  determined  in  the  following  manner. 
The  disturbing  force,  whether  radial  or  tangential,  acting  on  the  ring  may  be  conceived  to 
vary  from  one  satellite  to  another,  and  to  be  different  at  different  times.  It  is  therefore  a 
perfectly  arbitrary  function  of  s  and  t. 

Let  Fourier's  method  be  applied  to  the  general  disturbing  force  so  as  to  divide  it  up  into 
terms  depending  on  periodic  functions  of  s,  so  that  each  term  is  of  the  form  F  (t)  cos  (ms  +  a), 
where  the  function  of  t  is  still  perfectly  arbitrary. 

But  it  appears  from  the  general  theory  of  the  permanent  motions  of  the  heavenly  bodies 
that  they  may  all  be  expressed  by  periodic  functions  of  t  arranged  in  series.  Let  vt  be  the 
argument  of  one  of  these  terms,  then  the  corresponding  term  of  the  disturbance  will  be  of  the 
form  P  cos  (ms  +  vt  +  a). 

This  term  of  the  disturbing  force  indicates  an  alternately  positive  and  negative  action, 
disposed  in  m  waves  round  the  ring,   completing  its  period  relatively  to  each  particle  in  the 

27T 

time  —  ,  and  travelling  as  a   wave  among  the  particles  with  an  angular  velocity  -  — ,    the 
v  *>, 
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V 

angular  velocity  relative  to  fixed  space  being  of  course  w .      The  whole  disturbing  force 

may  be  split  up  into  terms  of  this  kind. 

17.  Each  of  these  elementary  disturbances  will  produce  its  own  wave  in  the  ring, 
independent  of  those  which  belong  to  the  ring  itself.  This  new  wave,  due  to  external 
disturbance,  and  following  different  laws  from  the  natural  waves  of  the  ring,  is  called  the 
forced  wave.  The  angular  velocity  of  the  forced  wave  is  the  same  as  that  of  the  disturbing 
force,  and  its  maxima  and  minima  coincide  with  those  of  the  force,  but  the  extent  of  the 
disturbance  and  its  direction  depend  on  the  comparative  velocities  of  the  forced  wave  and  the 
four  natural  waves. 

When  the  velocity  of  the  forced  wave  lies  between  the  velocities  of  the  two  middle  free 
waves,  or  is  greater  than  that  of  the  swiftest,  or  less  than  that  of  the  slowest,  then  the  radial 
displacement  due  to  a  radial  disturbing  force  is  in  the  same  direction  as  the  force,  but  the 
tangential  displacement  due  to  a  tangential  disturbing  force  is  in  the  opposite  direction  to  the 
force. 

The  radial  force  therefore  in  this  case  produces  a  positive  forced  wave,  and  the  tangential 
force  a  negative  forced  wave. 

When  the  velocity  of  the  forced  wave  is  either  between  the  velocities  of  the  first  and  second 
free  waves,  or  between  those  of  the  third  and  fourth,  then  the  radial  disturbance  produces 
a  forced  wave  in  the  contrary  direction  to  that  in  which  it  acts,  or  a  negative  wave,  and  the 
tangential  force  produces  a  positive  wave. 

The  coefficient  of  the  forced  wave  changes  sign  whenever  its  velocity  passes  through  the 
value  of  any  of  the  velocities  of  the  free  waves,  but  it  does  so  by  becoming  infinite,  and  not  by 
vanishing,  so  that  when  the  angular  velocity  very  nearly  coincides  with  that  of  a  free  wave, 
the  forced  wave  becomes  very  great,  and  if  the  velocity  of  the  disturbing  force  were  made 
exactly  equal  to  that  of  a  free  wave,  the  coefficient  of  the  forced  wave  would  become  infinite. 
In  such  a  case  we  should  have  to  readjust  our  approximations,  and  to  find  whether  such  a 
coincidence  might  involve  a  physical  impossibility. 

The  forced  wave  which  we  have  just  investigated  is  that  which  would  maintain  itself  in 
the  ring,  supposing  that  it  had  been  set  agoing  at  the  commencement  of  the  motion.  It  is  in 
fact  the  form  of  dynamical  equilibrium  of  the  ring  under  the  influence  of  the  given  forces.  In 
order  to  find  the  actual  motion  of  the  ring  we  must  combine  this  forced  wave  with  all  the  free 
waves,  which  go  on  independently  of  it,  and  in  this  way  the  solution  of  the  problem  becomes 
perfectly  complete,  and  we  can  determine  the  whole  motion  under  any  given  initial  circum- 
stances, as  we  did  in  the  case  where  no  disturbing  force  acted. 

For  instance,  if  the  ring  were  perfectly  uniform  and  circular  at  the  instant  when  the 
disturbing  force  began  to  act,  we  should  have  to  combine  with  the  constant  forced  wave  a 
system  of  four  free  waves  so  disposed,  that  at  the  given  epoch,  the  displacements  due  to  them 
should  exactly  neutralize  those  due  to  the  forced  wave.  By  the  combined  effect  of  these  four 
free  waves  and  the  forced  one  the  whole  motion  of  the  ring  would  be  accounted  for,  beginning 
from  its  undisturbed  state. 

F 
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The  disturbances  which  are  of  most  importance  in  the  theory  of  Saturn's  rings  are  those 
which  are  produced  in  one  ring  by  the  action  of  attractive  forces  arising  from  waves  belonging 
to  another  ring. 

The  effect  of  this  kind  of  action  is  to  produce  in  each  ring,  besides  its  own  four  free  waves, 
four  forced  waves  corresponding  to  the  free  waves  of  the  other  ring.  There  will  thus  be 
eight  waves  in  each  ring,  and  the  corresponding  waves  in  the  two  rings  will  act  and  react  on 
each  other,  so  that,  strictly  speaking,  every  one  of  the  waves  will  be  in  some  measure  a 
forced  wave,  although  the  system  of  eight  waves  will  be  the  free  motion  of  the  two  rings  taken 
together.  The  theory  of  the  mutual  disturbance  and  combined  motion  of  two  concentric  rings 
of  satellites  requires  special  consideration. 

18.  On  the  motion  of  a  ring  of  satellites  when  the  conditions  of  stability  are  not 
fulfilled. 

We  have  hitherto  been  occupied  with  the  case  of  a  ring  of  satellites,  the  stability  of  which 
was  ensured  by  the  smallness  of  mass  of  the  satellites  compared  with  that  of  the  central  body. 
We  have  seen  that  the  statically  unstable  condition  of  each  satellite  between  its  two  immediate 
neighbours  may  be  compensated  by  the  dynamical  effect  of  its  revolution  round  the  planet,  and 
a  planet  of  sufficient  mass  can  not  only  direct  the  motion  of  such  satellites  round  its  own  body, 
but  can  likewise  exercise  an  influence  over  their  relations  to  each  other,  so  as  to  overrule  their 
natural  tendency  to  crowd  together,  and  distribute  and  preserve  them  in  the  form  of  a  ring. 

We  have  traced  the  motion  of  each  satellite,  the  general  shape  of  the  disturbed  ring,  and 
the  motion  of  the  various  waves  of  disturbance  round  the  ring,  and  determined  the  laws  both 
of  the  natural  or  free  waves  of  the  ring,  and  of  the  forced  waves,  due  to  extraneous  disturbing 
forces. 

We  have  now  to  consider  the  cases  in  which  such  a  permanent  motion  of  the  ring  is 
impossible,  and  to  determine  the  mode  in  which  a  ring,  originally  regular,  will  break  up,  in 
the  different  cases  of  instability. 

The  equation  from  which  we  deduce  the  conditions  of  stability  is — 

U=  w4  -  loo"  +  - R  (2K  +  L  -  N)\  n*  -  4w-RMn  + 

[3o/--R(2K+L)-Rn]  ~^RsM*~o. 

[ft.  M  J  M 

The  quantity,  which,  in  the  critical  cases,  determines  the  nature  of  the  roots  of  this 
equation,  is  N.  The  quantity  M  in  the  third  term  is  always  small  compared  with  L  and  N 
when  m  is  large,  that  is,  in  the  case  of  the  dangerous  short  waves.  We  may  therefore  begin 
our  study  of  the  critical  cases  by  leaving  out  the  third  term.  The  equation  then  becomes  a 
quadratic  in  ns,  and  in  order  that  all  the  values  of  n  may  be  real,  both  values  of  n2  must  be 
real  and  positive. 

The  condition  of  the  values  of  n2  being  real  is 

wi  +  w*-R(*K+2L-uN)  +  -Ri('2K +  L  +  Nf>0     ....     (6l), 
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which  shows  that  a>2  must  either  be  about  14  times  at  least  smaller,  or  about  14  times  at  least 
greater,  than  quantities  like  —  RN. 

That  both  values  of  w2  may  be  positive,  we  must  have, 

w"+-R  (2A'  +  L  -  N)  >  0 

«*)• 


{3a,2  -   -  R  (°K  +  L)\  -RN>0 


We  must  therefore  take  the  larger  values  of  o>2,  and  also  add  the  condition  that  N  be 
positive. 

RN 

We  may  therefore  state  roughly,  that,  to  ensure  stability, ,  the  coefficient  of  tangential 

M 
attraction,   must  lie  between   zero  and   -jL  a>2.      If  the  quantity   be  negative,  the  two  small 
values  of  n  will  become  pure  impossible  quantities.     If  it  exceed  ^  a>2,  all  the  values  of  n 
will  take  the  form  of  mixed  impossible  quantities. 

If  we  write  x  for  —  RN,  and  omit  the  other  disturbing  forces,  the  equation  becomes 

U  =  n*  -{<1)~-x)n2  +  3w2a!  =  0    ..." (63), 


whence  n2  =  1  (to2  -  x)  ±  1  \Zw*  -14a)2*  +  x2 (64). 

If  x  be  small,  two  of  the  values  of  n  are  nearly  ±  a,,  and  the  others  are  small  quantities, 
real  when  .r  is  positive  and  impossible  when  x  is  negative. 

t 

If  x  be  greater  than  (7  —  v  48)  to2,  or  —  nearly,  the  term  under  the  radical  becomes  negative, 

14 

and  the  value  of  n  becomes 

n  =  ±  1  s/y/lZoSx  +  w"  -  X  ±  1  y/~l  >/\/l2w*X  -a,2  +  X      .       .       .       (65), 
where  one  of  the  terms  is  a  real  quantity,  and   the  other  impossible.      Every  solution  may  be 
put  under  the  form 

n—p^^Z-lq (^)> 

where  q  =  0  for  the  case  of  stability,  p  =  0  for   the  pure  impossible  roots,  and  p  and  q  finite 
for  the  mixed  roots. 

Let  us  now  adopt  this  general  solution  of  the  equation  for  n,  and  determine  its  mechanical 
significance  by  substituting  for  the  impossible  circular  functions  their  equivalent  real  ex- 
ponential functions. 

Substituting  the  general  value  of  n  in  equations  (34),  (35), 

p  =  A  [cos  {ms  +  (p  +  */ -\q)  t  +  a]  +  cos  {ms  +  (p  -  y/-lq)t  +  a}]      .      (67), 

**>(p  +  y/~  1  q)    .  — 

; / A2 sin  \ms  +  (p  +  V  -  1  q)  t  +  a\ 

(p  +  V-lq)2+x         '  ' 

2ai  (p  -  \/-  1  q)  

-^-V^7<7)2+*sin^  +  (p-^-l9)'  +  a| 


=  -A 


(68). 
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Introducing  the  exponential  notation,  these  values  become 

p  =  A  (eqi  +  £"«')  cos  (ms  +  pt  +  a) {69), 

2wA  \     p(p~  +  q"  +  x)  (eql  4-  e~")  sin  (ms  +  pt  +  a)  \ 


f     p(p2  +  q  +  x)  (eql  +  e~qi)  sin  (ms  +  pt  +  a)  1 
(p'~  +  q*y+  2(p2  -  q9)x  +  x"  \  +  q(p2  +  q2  -  w)  (eqt  -  e~*1)  cos  (ms  +  pt  +  a)  )' 

We  have  now  obtained  a   solution   free  from   impossible  quantities,  and   applicable  to  every 
case. 

When  q  =  0,  the  case  becomes  that  of  real  roots,  which  we  have  already  discussed.  When 
p  =  0,  we  have  the  case  of  pure  impossible  roots  arising  from  the  negative  values  of  n2.  The 
solutions  corresponding  to  these  roots  are 

p  =  A  («**  +  e"")  cos  (ms  +  a) (71). 

ff=  _?^   (««'  -  e-««)  cos  (m«  +  a) (72). 

q*-x 

The  part  of  the  coefficient  depending  on  e"'1  diminishes  indefinitely  as  the  time  increases, 
and  produces  no  marked  effect.  The  other  part,  depending  on  eqt,  increases  in  a  geometrical 
proportion  as  the  time  increases  arithmetically,  and  so  breaks  up  the  ring.  In  the  case  of  .r 
being  a  small  negative  quantity,  q~  is  nearly  Sx,  so  that  the  coefficient  of  a  becomes 

-  3- A. 
9 

It  appears  therefore  that   the  motion  of    each   particle  is  either  outwards  and   backwards 

or  inwards  and   forwards,   but   that  the  tangential  part  of    the    motion    greatly  exceeds    the 

normal  part. 

It  may  seem  paradoxical  that  a  tangential  force,  acting  towards  a  position  of  equilibrium, 
should  produce  instability,  while  a  small  tangential  force  from  that  position  ensures  stability, 
but  it  is  easy  to  trace  the  destructive  tendency  of  this  apparently  conservative  force. 

Suppose  a  particle  slightly  in  front  of  a  crowded  part  of  the  ring,  then  if  w  is  negative 
there  will  be  a  tangential  force  pushing  it  forwards,  and  this  force  will  cause  its  distance  from 
the  planet  to  increase,  its  angular  velocity  to  diminish,  and  the  particle  itself  to  fall  back  on 
the  crowded  part,  thereby  increasing  the  irregularity  of  the  ring,  till  the  whole  ring  is  broken 
up.  In  the  same  way  it  may  be  shown  that  a  particle  behind  a  crowded  part  will  be  pushed 
into  it.  The  only  force  which  could  preserve  the  ring  from  the  effect  of  this  action,  is  one 
which  would  prevent  the  particle  from  receding  from  the  planet  under  the  influence  of  the  tan- 
gential force,  or  at  least  prevent  the  diminution  of  angular  velocity.  The  transversal  force  of 
attraction  of  the  ring  is  of  this  kind,  and  acts  in  the  right  direction,  but  it  can  never  be  of  suffi- 
cient magnitude  to  have  the  required  effect.  In  fact  the  thing  to  be  done  is  to  render  the  least 
term  of  the  equation  in  n?  positive  when  iVis  negative,  which  requires 

-  RCiK  +  L)  >  3m2, 
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and  this  condition  is  quite  inconsistent  with  any  constitution  of  the  ring  which  fulfils  the  other 
condition  of  stability  which  we  shall  arrive  at  presently. 

We  may  observe  that  the  waves  belonging  to  the  two  real  values  of  n,  ±  w,  must  be 
conceived  to  be  travelling  round  the  ring  during  the  whole  time  of  its  breaking  up,  and 
conducting  themselves  like  ordinary  waves,  till  the  excessive  irregularities  of  the  ring  become 
inconsistent  with  their  uniform  propagation. 

The  irregularities  which  depend  on  the  exponential  solutions  do  not  travel  round  the  rino- 
by  propagation  among  the  satellites,  but  remain  among  the  same  satellites  which  first  began  to 
move  irregularly. 

We  have  seen  the  fate  of  the  ring  when  x  is  negative.  When  w  is  small  we  have  two 
small  and  two  large  values  of  n,  which  indicate  regular  waves,  as  we  have  already  shown.  As 
x  increases,  the  small  values  of  n  increase,  and  the  large  values  diminish,  till  they  meet  and 
form  a  pair  of  positive  and  a  pair  of  negative  equal  roots,  having  values  nearly  ±  •6Sai.  When 
x  becomes  greater  than  about  y±u?,  then  all  the  values  of  n  become  impossible,  of  the  form 
p  +  V    -  1  q,  q  being  small  when  as  first  begins  to  exceed  its  limits,  and  p  being  nearly  ±  ^Sco. 

The  values  of  p  and  a  indicate  periodic  inequalities  having  the  period  — ,  but  increasing 

P 
in  amplitude  at  a  rate  depending  on  the  exponential  eqi.  At  the  beginning  of  the  motion  the 
oscillations  of  the  particles  are  in  ellipses  as  in  the  case  of  stability,  having  the  ratio  of  the 
axes  about  1  in  the  normal  direction  to  3  in  the  tangential  direction.  As  the  motion  continues, 
these  ellipses  increase  in  magnitude,  and  another  motion  depending  on  the  second  term  of  a  is 
combined  with  the  former,  so  as  to  increase  the  ellipticity  of  the  oscillations  and  to  turn  the 
major  axis  into  an  inclined  position,  so  that  its  fore  end  points  a  little  inwards,  and  its  hinder 
end  a  little  outwards.  The  oscillations  of  each  particle  round  its  mean  position  are  therefore 
in  ellipses,  of  which  both  axes  increase  continually  while  the  eccentricity  increases,  and  the  major 
axis  becomes  slightly  inclined  to  the  tangent,  and  this  goes  on  till  the  ring  is  destroyed.  In 
the  mean  time  the  irregularities  of  the  ring  do  not  remain  among  the  same  set  of  particles  as  in 

...  p 

the  former  case,  but   travel  round   the  ring   with  a  relative  angular   velocity .      Of  these 

m 

waves  there  are  four,  two  travelling  forwards  among  the  satellites,  and  two  travelling  backwards. 

One  of  each  of  these  pairs  depends  on   a  negative  value  of  q,  and   consists  of  a   wave  whose 

amplitude  continually   decreases.       The  other  depends  on  a  positive  value  of  q,  and   is   the 

destructive  wave  whose  character  we  have  just  described. 

19.  We  have  taken  the  case  of  a  ring  composed  of  equal  satellites,  as  that  with  which  we 
may  compare  other  cases  in  which  the  ring  is  constructed  of  loose  materials  differently  arranged. 

In  the  first  place  let  us  consider  what  will  be  the  conditions  of  a  ring  composed  of  satellites 
of  unequal  mass.  We  shall  find  that  the  motion  is  of  the  same  kind  as  when  the  satellites  are 
equal. 
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For  by  arranging  the  satellites  so  that  the  smaller  satellites  are  closer  together  than  the 
larger  ones,  we  may  form  a  ring  which  will  revolve  uniformly  about  Saturn,  the  resultant 
force  on  each  satellite  being  just  sufficient  to  keep  it  in  its  orbit. 

To  determine  the  stability  of  this  kind  of  motion,  we  must  calculate  the  disturbing  forces 
due  to  any  given  displacement  of  the  ring.  This  calculation  will  be  more  complicated  than  in 
the  former  case,  but  will  lead  to  results  of  the  same  general  character.  Placing  these  forces  in 
the  equations  of  motion,  we  shall  find  a  solution  of  the  same  general  character  as  in  the  former 
case,  only  instead  of  regular  waves  of  displacement  travelling  round  the  ring,  each  wave  will 
be  split  and  reflected  when  it  comes  to  irregularities  in  the  chain  of  satellites.  But  if  the 
condition  of  stability  for  every  kind  of  wave  be  fulfilled,  the  motion  of  each  satellite  will  consist 
of  small  oscillations  about  its  position  of  dynamical  equilibrium,  and  thus,  on  the  whole,  the 
rino-  will  of  itself  assume  the  arrangement  necessary  for  the  continuance  of  its  motion,  if  it  be 
originally  in  a  state  not  very  different  from  that  of  equilibrium. 

20.  We  now  pass  to  the  case  of  a  ring  of  an  entirely  different  construction.  It  is  possible 
to  conceive  of  a  quantity  of  matter,  either  solid  or  liquid,  not  collected  into  a  continuous  mass, 
but  scattered  thinly  over  a  great  extent  of  space,  and  having  its  motion  regulated  by  the 
gravitation  of  its  parts  to  each  other,  or  towards  some  dominant  body.  A  shower  of  rain,  hail, 
or  cinders  is  a  familiar  illustration  of  a  number  of  unconnected  particles  in  motion  ;  the  visible 
stars,  the  milky  way,  and  the  resolved  nebulas,  give  us  instances  of  a  similar  scattering  of  bodies 
on  a  larger  scale.  In  the  terrestrial  instances  we  see  the  motion  plainly,  but  it  is  governed  by 
the  attraction  of  the  earth,  and  retarded  by  the  resistance  of  the  air,  so  that  the  mutual  attrac- 
tion of  the  parts  is  completely  masked.  In  the  celestial  cases  the  distances  are  so  enormous, 
and  the  time  during  which  they  have  been  observed  so  short,  that  we  can  perceive  no  motion 
at  all.  Still  we  are  perfectly  able  to  conceive  of  a  collection  of  particles  of  small  size  compared 
with  the  distances  between  them,  acting  upon  one  another  only  by  the  attraction  of  gravitation, 
and  revolving  round  a  central  body.  The  average  density  of  such  a  system  may  be  smaller 
than  that  of  the  rarest  gas,  while  the  particles  themselves  may  be  of  great  density  ;  and  the 
appearance  from  a  distance  will  be  that  of  a  cloud  of  vapour,  with  this  difference,  that  as  the 
space  between  the  particles  is  empty,  the  rays  of  light  will  pass  through  the  system  without 
being  refracted,  as  they  would  have  been  if  the  system  had  been  gaseous. 

Such  a  system  will  have  an  average  density  which  may  be  greater  in  some  places  than 
others.  The  resultant  attraction  will  be  towards  places  of  greater  average  density,  and  thus 
the  density  of  those  places  will  be  increased  so  as  to  increase  the  irregularities  of  density. 
The  system  will  therefore  be  statically  unstable,  and  nothing  but  motion  of  some  kind  can 
prevent  the  particles  from  forming  agglomerations,  and  these  uniting,  till  all  are  reduced  to  one 
solid  mass. 

We  have  already  seen  how  dynamical  stability  can  exist  where  there  is  statical  instability 
in  the  case  of  a  row  of  particles  revolving  round  a  central  body.  Let  us  now  conceive  a  cloud 
of  particles  forming  a  ring  of  nearly  uniform  density  revolving  about  a  central  body.  There 
will  be  a  primary   effect  of    inequalities  in    density  tending    to    draw   particles    towards    the 
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denser  parts  of  the  ring,  and  this  will  elicit  a  secondary  effect,  due  to  the  motion  of  revolution, 
tending  in  the  contrary  direction,  so  as  to  restore  the  rings  to  uniformity.  The  relative  magni- 
tude of  these  two  opposing  forces  determines  the  destruction  or  preservation  of  the  ring. 

To  calculate  these  effects  we  must  begin  with  the  statical  problem  : — To  determine  the 
forces  arising  from  the  given  displacements  of  the  ring. 

The  longitudinal  force  arising  from  longitudinal  displacements  is  that  which  has  most  effect 
in  determining  the  stability  of  the  ring.  In  order  to  estimate  its  limiting  value  we  shall  solve 
a  problem  of  a  simpler  form. 

21.  An  infinite  mass,  originally  of  uniform  density  k,  has  its  particles  displaced  by  a 
quantity  £  parallel  to  the  axis  of  x,  so  that  £  =  A  cos  mx,  to  determine  the  attraction  on  each 
particle  due  to  this  displacement. 

The  density  at  any  point  will  differ  from  the  original  density  by  a  quantity  k',  so  that 

(k  +  k')  (dx  +  d%)  =  kdx (73). 

df 

k'  =  —  k  — —  =  Akm  sin  mx (74). 

ax 

The  potential  at  any  point  will  be  V  +  Y ,  where  V  is  the  original  potential,  and  V  depends  on 
the  displacement  only,  so  that 

d?r      d*V      d*V  ,,  .    „ 

—FT  +   j  ,   +  -T^r  +  4?rA;  =  0 (75). 

ax'         dy-        d%~ 

Now  V'  is  a  function  of  x  only,  and  therefore, 

V  =  4-irAk  —  sin  mx (76), 

m 

and  the  longitudinal  force  is  found  by  differentiating  V'  with  respect  to  x. 

dV 
X  =  ——  =  4-nkA  cos  mx  =  4-n-k^ (77). 

Now  let  us  suppose  this  mass  not  of  infinite  extent,  but  of  finite  section  parallel  to  the  plane 
of  yz.  This  change  amounts  to  cutting  off  all  portions  of  the  mass  beyond  a  certain  boundary. 
Now  the  effect  of  the  portion  so  cut  off  upon  the  longitudinal  force  depends  on  the  value  of  m. 
When  m  is  large,  so  that  the  wave-length  is  small,  the  effect  of  the  external  portion  is  insensible, 
so  that  the  longitudinal  force  due  to  short  waves  is  not  diminished  by  cutting  off  a  great  portion 
of  the  mass. 

22.  Applying  this  result  to  the  case  of  a  ring,  and  putting  s  for  x,  and  a  for  £,  we  have 

a  =  A  cos  ms,  and  T  =  4>irkA  cos  ms, 

so  that  -  RN  =  4nrk, 

when  m  is  very  large,  and  this  is  the  greatest  value  of  N. 
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The  value  of  L  has  little  effect  on  the  condition  of  stability.       If  L  and   M  are   both 

neglected,  that  condition  is  to2  >  27'856  (2nk) (78), 

and  if  L  be  as  much  as  AA7,  then  o>8  >  25649  (2irk) (79), 

so  that  it  is  not  important  whether  we  calculate  the  value  of  L  or  not. 

The  condition  of  stability  is,  that  the  average  density  must  not  exceed  a  certain  value. 
Let  us  ascertain  the  relation  between  the  maximum  density  of  the  ring  and  that  of  the  planet. 

Let  b  be  the  radius  of  the  planet,  that  of  the  ring  being  unity,  then  the  mass  of  Saturn  is 
4r7r63&'  =  to2  if  k'  be  the  density  of  the  planet.  If  we  assume  that  the  radius  of  (he  ring  is 
twice  that  of  the  planet,  as  Laplace  has  done,  then  b  =  -i  and 

k' 

-  =  334-2  to  3077 (80), 

so  that  the  density  of  the  ring  cannot  exceed  g^  of  that  of  the  planet.  Now  Laplace  has 
shown  that  if  the  outer  and  inner  parts  of  the  ring  have  the  same  angular  velocity,  the  ring  will 
not  hold  together  if  the  ratio  of  the  density  of  the  planet  to  that  of  the  ring  exceeds  1*3,  so  that 
in  the  first  place,  our  ring  cannot  have  uniform  angular  velocity,  and  in  the  second  place, 
Laplace's  ring  cannot  preserve  its  form,  if  it  is  composed  of  loose  materials  acting  on  each  other 
only  by  the  attraction  of  gravitation,  and  moving  with  the  same  angular  velocity  throughout. 

23.  On  the  forces  arising  from  inequalities  of  thickness  in  a  thin  stratum  of  fluid  of  in- 
definite extent. 

The  forces  which  act  on  any  portion  of  a  continuous  fluid  are  of  two  kinds,  the  pressures 
of  contiguous  portions  of  fluid,  and  the  attractions  of  all  portions  of  the  fluid  whether  near  or 
distant.  In  the  case  of  a  thin  stratum  of  fluid,  not  acted  on  by  any  external  forces,  the 
pressures  are  due  mainly  to  the  component  of  the  attraction  which  is  perpendicular  to  the 
plane  of  the  stratum.  It  is  easy  to  show  that  a  fluid  acted  on  by  such  a  force  will  tend  to 
assume  a  position  of  equilibrium,  in  which  its  free  surface  is  plane;  and  that  any  irregularities 
will  tend  to  equalise  themselves,  so  that  the  plane  surface  will  be  one  of  stable  equilibrium. 

It  is  also  evident,  that  if  we  consider  only  that  part  of  the  attraction  which  is  parallel  to 
the  plane  of  the  stratum,  we  shall  find  it  always  directed  towards  the  thicker  parts,  so  that  the 
effect  of  this  force  is  to  draw  the  fluid  from  thinner  to  thicker  parts,  and  so  to  increase  irre- 
gularities and  destroy  equilibrium. 

The  normal  attraction  therefore  tends  to  preserve  the  stability  of  equilibrium,  while  the 
tangential  attraction  tends  to  render  equilibrium  unstable. 

According  to  the  nature  of  the  irregularities  one  or  other  of  these  forces  will  prevail,  so 
that  if  the  extent  of  the  irregularities  is  small,  the  normal  forces  will  ensure  stability,  while, 
if  the  inequalities  cover  much  space,  the  tangential  forces  will  render  equilibrium  unstable,  and 
break  up  the  stratum  into  beads. 
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To  fix  our  ideas,  let  us  conceive  the  irregularities  of  the  stratum  split  up  into  the  form 
of  a  number  of  systems  of  waves  superposed  on  one  another,  then,  by  what  we  have  just 
said,  it  appears,  that  very  short  waves  will  disappear  of  themselves,  and  be  consistent  with 
stability,  while  very  long  waves  will  tend  to  increase  in  height,  and  will  destroy  the  form  of 
the  stratum. 

In  order  to  determine  the  law  according  to  which  these  opposite  effects  take  place,  we 
must  subject  the  case  to  mathematical  investigation. 

Let  us  suppose  the  fluid  incompressible,  and  of  the  density  k,  and  let  it  be  originally  con- 
tained between  two  parallel  planes,  at  distances  +  c  and  -  c  from  that  of  (ay),  and  extending 
to  infinity.  Let  us  next  conceive  a  series  of  imaginary  planes,  parallel  to  the  plane  of  (?/«),  to 
be  plunged  into  the  fluid  stratum  at  infinitesimal  distances  from  one  another,  so  as  to  divide 
the  fluid  into  imaginary  slices  perpendicular  to  the  plane  of  the  stratum. 

Next  let  these  planes  be  displaced  parallel  to  the  axis  of  eo  according  to  this  law — that  if 
w  be  the  original  distance  of  the  plane  from  the  origin,  and  £  its  displacement  in  the  direction 
of  x, 

£  =  A  cos  m.v (81). 

According  to  this  law  of  displacement,  certain  alterations  will  take  place  in  the  distances 
between  consecutive  planes;  but  since  the  fluid  is  incompressible,  and  of  indefinite  extent  in 
the  direction  of  y,  the  change  of  dimension  must  occur  in  the  direction  of  z.  The  original 
thickness  of  the  stratum  was  2c.  Let  its  thickness  at  any  point  after  displacement  be  2c  +  2£, 
then  we  must  have 

(2c  +  2p  (l  +  gj   =  2c (82), 

df 

or  t  =  —  c  —  =  cmA  sin  mm (83). 

5  dx 

Let  us  assume  that  the  increase  of  thickness  2^  is  due  to  an  increase  of  £  at  each  surface; 
this  is  necessary  for  the  equilibrium  of  the  fluid  between  the  imaginary  planes. 

We  have  now  produced  artificially,  by  means  of  these  planes,  a  system  of  waves  of  longitu- 

.    2tt  . 

dinal  displacement  whose  length  is  —  and  amplitude  A;  and  we  have  found  that  this  has  pro- 
duced a  system  of  waves  of  normal  displacement  on  each  surface,  having  the  same  length,  with 
a  height  =  cmA. 

In  order  to  determine  the  forces  arising  from  these  displacements,  we  must,  in  the  first 
place,  determine  the  potential  function  at  any  point  of  space,  and  this  depends  partly  on  the 
state  of  the  fluid  before  displacement,  and  partly  on  the  displacement  itself.  We  have,  in  all 
cases — 

d*V      d?V      d2V  ,    „ 

M  +  df  +  d*?  =  -4irp (84)- 

G 


42  ON   THE   STABILITY   OF   THE   MOTION 

Within  the  fluid,  p  =  k;  beyond  it,  p  =  0. 
Before  displacement,  the  equation  is  reduced  to 

drV 

37-~** (85)- 

Instead  of  assuming  V  =  0  at  infinity,  we  shall  assume  V  =  0    at   the  origin,  and  since  in 
this  case  all  is  symmetrical,  we  have 

dV  , 

within    the   fluid  V,  =  -  27r&#2,  —  =  -  4>irkz 

dz 


at    the   bounding    planes  V  =  -  2tt&c2    — —  =  =p  47T&C    \-         ....  (86) ; 

beyond   them  V2  =  27T&C  (=p  2»  -  c)  —  =  =p  47T&C 

the  upper  sign  being  understood  to  refer  to  the  boundary  at  distance  +  c,  and  the  lower  to  the 
boundary  at  distance  -  c  from  the  origin. 

Havino-  ascertained  the  potential  of  the  undisturbed  stratum,  we  find  that  of  the  disturbance 
by  calculating  the  effect  of  a  stratum  of  density  k  and  thickness  £,  spread  over  each  surface 
according  to  the  law  of  thickness  already  found.  By  supposing  the  coefficient  A  small  enough, 
(as  we  may  do  in  calculating  the  displacements  on  which  stability  depends,)  we  may  diminish 
the  absolute  thickness  indefinitely,  and  reduce  the  case  to  that  of  a  mere  "superficial  density," 
such  as  is  treated  of  in  the  theory  of  electricity.  We  have  here,  too,  to  regard  some  parts  as 
of  negative  density ;  but  we  must  recollect  that  we  are  dealing  with  the  difference  between  a 
disturbed  and  an  undisturbed  system,  which  may  be  positive  or  negative,  though  no  real  mass 
can  be  negative. 

Let  us  for  an  instant  conceive  only  one  of  these  surfaces  to  exist,  and  let  us  transfer  the 
origin  to  it.      Then  the  law  of  thickness  is 

£=  mcA  sin  moo (83), 

and  we  know  that  the  normal  component  of  attraction  at  the  surface  is  the  same  as  if  the  thick- 
ness had  been  uniform  throughout,  so  that 

dV  ,Y 

oar  s 

on  the  positive  side  of  the  surface. 
Also,  the  solution  of  the  equation 

<PF     d'V 

—  +  —  =o, 

dtv2       dz2 
consists  of  a  series  of  terms  of  the  form  Ce"  sin  ix 

Of  these  the  only  one  with  which  we  have  to  do  is  that  in  which  i  =  -  m.  Applying  the 
condition  as  to  the  normal  force  at  the  surface,  we  get 
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V  =  2irkce~mzA&mmx (87) 

for  the  potential  on  the  positive  side  of  the  surface,  and 

V  =  2irkcem'A  sin  mx (88). 

on  the  negative  side. 

Calculating  the  potentials  of  a  pair  of  such  surfaces  at  distances  +  c  and  -  c  from  the 
plane  of  xy,  and  calling  V  the  sum  of  their  potentials,  we  have  for  the  space  between  these 
planes 

Vi  =  2trkcA  sin  mxe~mc  (emz  +  e-"'2) 


(89); 
beyond    them  V2  =  2irkcA  sin  nixe*"1'  (emc  +  e 

the  upper  or  lower  sign  of  the  index  being  taken  according  as  as  is  positive  or  negative. 

These  potentials  must  be  added  to  those  formerly  obtained,  to  get  the  potential  at  any 
point  after  displacement. 

We  have  next  to  calculate  the  pressure  of  the  fluid  at  any  point,  on  the  supposition  that 
the  imaginary  planes  protect  each  slice  of  the  fluid  from  the  pressure  of  the  adjacent  slices,  so 
that  it  is  in  equilibrium  under  the  action  of  the  forces  of  attraction,  and  the  pressure  of  these 
planes  on  each  side.  Now  in  a  fluid  of  density  k,  in  equilibrium  under  forces  whose  potential 
is  V,  we  have  always — 

dp  -k- 

so  that  if  we  know  that  the  value  of  p  is  p0  where  that  of  V  is  V0,  then  at  any  other  point 

p=Po  +  k(V-  Va). 
Now,  at  the  free  surface  of  the  fluid,  p  =  0,  and  the  distance  from  the  free  surface  of  the 
disturbed  fluid  to  the  plane  of  the  original  surface  is  £,  a  small  quantity.  The  attraction 
which  acts  on  this  stratum  of  fluid  is,  in  the  first  place,  that  of  the  undisturbed  stratum,  and 
this  is  equal  to  4nrkc,  towards  that  stratum.  The  pressure  due  to  this  cause  at  the  level  of 
the  original  surface  will  be  4ttA;26'£',  and  the  pressure  arising  from  the  attractive  forces  due  to 
the  displacements  upon  this  thin  layer  of  fluid,  will  be  small  quantities  of  the  second  order, 
which  we  neglect.      We  thus  find  the  pressure  when  z  =  c  to  be, 

p0  =  4-Trk'2c*mJ  sin  mx. 
The  potential  of  the  undisturbed  mass  when  z  =  c  is 

V0  =  -  2wkc"; 
and  the  potential  of  the  disturbance  itself  for  the  same  value  of  z,  is 

V0'  =  ZnkcA  sin  mx  (l  +  e-2""'). 
So  that  we  find  the  general  value  ofp  at  any  other  point  to  be 

p  =  2ttZ;-  (c2  -  %2)  +  27rk2cA  sin  mx  {2cm  -  1  -  e-8mf  +  e  "">  (e'"z  +  e''"')  \       (90). 
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This  expression  gives  the  pressure  of  the  fluid  at  any  point,  as  depending  on  the  state  of 
constraint  produced  by  the  displacement  of  the  imaginary  planes.  The  accelerating  effect  of 
these  pressures  on  any  particle,  if  it  were  allowed  to  move  parallel  to  x,  instead  of  being  confined 
by  the  planes,  would  be, 

1  dp 

k  dx  ' 

The  accelerating  effect  of  the  attractions  in  the  same  direction  is 

dV 
dx 

so  that  the  whole  acceleration  parallel  to  x  is 

X  =  -  SirkmcA  cosw.r  (2mc  -  e~*mc  -  1) (91). 

It  is  to  be  observed,  that  this  quantity  is  independent  of  z,  so  that  every  particle  in  the 
slice,  by  the  combined  effect  of  pressure  and  attraction,  is  urged  with  the  same  force,  and,  if 
the  imaginary  planes  were  removed,  each  slice  would  move  parallel  to  itself  without  distortion, 
as  long  as  the  absolute  displacements  remained  small.  We  have  now  to  consider  the  direction 
of  the  resultant  force  X,  and  its  changes  of  magnitude. 

We  must  remember  that  the  original  displacement  is  A  cos  mx,  if  therefore  (2me  -  e~imc-  l) 
be  positive,  X  will  be  opposed  to  the  displacement,  and  the  equilibrium  will  be  stable,  whereas 
if  that  quantity  be  negative,  X  will  act  along  with  the  displacement  and  increase  it,  and  so 
constitute  an  unstable  condition. 

It  may  be  seen  that  large  values  of  mc  give  positive  results  and  small  ones  negative.  The 
sign  changes  when 

2mc  =  1.147 (92), 

which  corresponds  to  a  wave-length 

2w 

\  =  2c- -  =  2c  (5.471) (93). 

1.147 

The  length  of  the  complete  wave  in  the  critical  case  is  5471  times  the  thickness  of  the 
stratum.      Waves  shorter  than  this  are  stable,  longer  waves  are  unstable. 

The  quantity  2mc  (2mc  -  e"2mc  -  l), 

has  a  minimum  when  2mc  =  '607 (94), 

and  the  wave  length  is  10-353  times  the  thickness  of  the  stratum. 

In   this   case  2roc  (2mc  -  e  -•"■  -  1)  =  -  '5°9 (95), 

and  X  =  -509TrkA  cosmx (96). 
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24.  Let  us  now  conceive  that  the  stratum  of  fluid,  instead  of  being  infinite  in  extent,  is 
limited  in  breadth  to  about  100  times  its  thickness.  The  pressures  and  attractions  will  not  be 
much  altered  by  this  removal  of  a  distant  part  of  the  stratum.  Let  us  also  suppose  that  this 
thin  but  broad  strip  is  bent  round  in  its  own  plane  into  a  circular  ring  whose  radius  is  more 
than  ten  times  the  breadth  of  the  strip,  and  that  the  waves,  instead  of  being  exactly  parallel  to 
each  other,  have  their  ridges  in  the  direction  of  radii  of  the  ring.  We  shall  then  have  trans- 
formed our  stratum  into  one  of  Saturn's  Rings,  if  we  suppose  those  rings  to  be  liquid,  and 
that  a  considerable  breadth  of  the  ring  has  the  same  angular  velocity. 

Let  us  now  investigate  the  conditions  of  stability  by  putting 

x  =  -  2irkmc  (2mc  -  e~2mc  -  1) 

2 

into  the  equation  for  n.      We  know   that  at  must  lie  between  0  and    —  to  ensure  stability. 

13-9  } 

Now  the  greatest  value  cf  x  in  the  fluid  stratum  is  'oOQirk.      Taking  Laplace's  ratio  of  the 

diameter  of  the  ring  to  that  of  the  planet,  this  gives  42-5  as  the  minimum  value  of  the  density 

of  the  planet  divided  by  that  of  the  fluid  of  the  ring. 

Now  Laplace  has  shown  that  any  value  of  this  ratio  greater  than  1"3  is  inconsistent  with 
the  rotation  of  any  considerable  breadth  of  the  fluid  at  the  same  angular  velocity,  so  that 
our  hypothesis  of  a  broad  ring  with  uniform  velocity  is  untenable. 

But  the  stability  of  such  a  ring  is  impossible  for  another  reason,  namely,  that  for  waves  in 
which  2mc>l  .147,  oc  is  negative,  and  the  ring  will  be  destroyed  by  these  short  waves  in  the 
manner  described  at  page  (36). 

When  the  fluid  ring  is  treated,  not  as  a  broad  strip,  but  as  a  filament  of  circular  or  elliptic- 
section,  the  mathematical  difficulties  are  very  much  increased,  but  it  may  be  shown  that  in  this 
case  also  there  will  be  a  maximum  value  of  x,  which  will  require  the  density  of  the  planet  to  be 
several  times  that  of  the  ring,  and  that  in  all  cases  short  waves  will  give  rise  to  negative  values 
of  x,  inconsistent  with  the  stability  of  the  ring. 

It  appears,  therefore,  that  a  ring  composed  of  a  continuous  liquid  mass  cannot  revolve  about 
a  central  body  without  being  broken  up,  but  that  the  parts  of  such  a  broken  ring  may,  under 
certain  conditions,  form  a  permanent  ring  of  satellites. 

On  the  Mutual  Perturbations  of  Two  Rings. 

25.  We  shall  assume  that  the  difference  of  the  mean  radii  of  the  rings  is  small  compared 
with  the  radii  themselves,  but  large  compared  with  the  distance  of  consecutive  satellites  of  the 
same  ring.     We  shall  also  assume  that  each  ring  separately  satisfies  the  conditions  of  stability. 

We  have  seen  that  the  effect  of  a  disturbing  force  on  a  ring  is  to  produce  a  series  of  waves 
whose  number  and  period  correspond  with  those  of  the  disturbing  force  which  produces 
them,  so  that  we  have  only  to  calculate  the  coefficient  belonging  to  the  wave  from  that  of  the 
disturbing  force. 
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Hence  in  investigating  the  simultaneous  motions  of  two  rings,  we  may  assume  that  the 
mutually  disturbing  waves  travel  with  the  same  absolute  angular  velocity,  and  that  a  maximum 
in  one  corresponds  either  to  a  maximum  or  a  minimum  of  the  other,  according  as  the  coefficients 
have  the  same  or  opposite  signs. 

Since  the  motions  of  the  particles  of  each  ring  are  affected  by  the  disturbance  of  the  other 
ring,  as  well  as  of  that  to  which  they  belong,  the  equations  of  motion  of  the  two  rings  will  be 
involved  in  each  other,  and  the  final  equation  for  determining  the  wave-velocity  will  have  eight 
roots  instead  of  four.  But  as  each  of  the  rings  has  four  free  waves,  we  may  suppose  these  to 
originate  forced  waves  in  the  other  ring,  so  that  we  may  consider  the  eight  waves  of  each  ring 
as  consisting  of  four  free  waves  and  four  forced  ones. 

In  strictness,  however,  the  wave-velocity  of  the  "free"  waves  will  be  affected  by  the  existence 
of  the  forced  waves  which  they  produce  in  the  other  ring,  so  that  none  of  the  waves  are  really 
"free'''  in  either  ring  independently,  though  the  whole  motion  of  the  system  of  two  rings  as  a 
whole  is  free. 

We  shall  find,  however,  that  it  is  best  to  consider  the  waves  first  as  free,  and  then  to 
determine  the  reaction  of  the  other  ring  upon  them,  which  is  such  as  to  alter  the  wave-velocity 
of  both,  as  we  shall  see. 

The  forces  due  to  the  second  ring  may  be  separated  into  three  parts. 
1st.      The  constant  attraction  when  both  rings  are  at  rest. 

2nd.      The  variation  of  the  attraction  on  the  first  ring,  due  to  its  own  disturbances. 
3rd.      The  variation  of  the  attraction  due  to  the  disturbances  of  the  second  ring. 

The  first  of  these  affects  only  the  angular  velocity.  The  second  affects  the  waves  of  each 
ring  independently,  and  the  mutual  action  of  the  waves  depends  entirely  on  the  third  class  of 
forces. 

26.      To  determine  the  attractions  between  two  rings. 

Let  R  and  a  be  the  mass  and  radius  of  the  exterior  ring,  R'  and  a  those  of  the  interior,  and 
let  all  quantities  belonging  to  the  interior  ring  be  marked  with  accented  letters.    (Fig  5.) 

1st.      Attraction  between  the  Rings  when  at  rest. 

Since  the  rings  are  at  a  distance  small  compared  with  their  radii,  we  may  calculate  the 
attraction  on  a  particle  of  the  first  ring  as  if  the  second  were  an  infinite  straight  line  at  distance 
a'  —  a  from  the  first. 

The  mass  of  unit  of  length  of  the  second  ring  is  -, ,  and  the  accelerating  effect  of  theattrac- 

2tto 

tion  of  such  a  filament  on  an  element  of  the  first  ring  is 

7?' 

— — ^  inwards (97). 

ira  (a  —  a) 
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The  attraction  of  the  first  ring  on  the  second  may  be  found  by  transposing  accented  and 
unaccented  letters. 

In  consequence  of  these  forces,  the  outer  ring  will  revolve  faster,  and  the  inner  ring  slower 
than  would  otherwise  be  the  case.  These  forces  enter  into  the  constant  terms  of  the  equations 
of  motion,  and  may  be  included  in  the  value  of  K. 

2nd.      Variation  due  to  disturbance  of  first  ring. 

If  we  put  a(l  +  p)  for  a  in  the  last  expression,  we  get  the  attraction  when  the  first  ring  is 
displaced.      The  part  depending  on  p  is 

R'a 


ira'(a  -  a'Y ' 


inwards (98)- 


This  is  the  only  variation  of  force  arising  from  the  displacement  of  the  first  ring.  It  affects 
the  value  of  L  in  the  equations  of  motion. 

3rd.       Variation  due  to  waves  in  the  second  ring. 

On  account  of  the  waves,  the  second  ring  varies  in  distance  from  the  first,  and  also  in  mass  of 
unit  of  length,  and  each  of  these  alterations  produces  variations  both  in  the  radial  and  tangential 
force,  so  that  there  are  four  things  to  be  calculated  : 

1st.  Radial  force  due  to  radial  displacement. 

2nd.  Radial  force  due  to  tangential  displacement. 

3rd.  Tangential  force  due  to  radial  displacement. 

4th.  Tangential  force  due  to  tangential  displacement. 

1st.      Put  a'(l  +  p')  for  a,  and  we  get  the  term  in  p 

R'      (2a  -  a)     .  . 

— >      ~r~, T*  P  mwards  =\p,  say (99)- 

■wa       («'  -  a)2  r  r        j 


2nd.     By  the  tangential  displacement  of  the  second   ring   the  section  is    reduced    in   the 
proportion  of  1  to  1  —  — ■,  and  therefore  there  is  an  alteration  of  the  radial  force  equal  to 

R'  da   .  .  i  da'  . 

—-,  inwards  =    -/a    — ^say, (100). 


ird{a  —  a)      ds'  ds', 

3rd.  By  the  radial  displacement  of  the  second  ring  the  direction  of  the  filament  near  the 
part  in  question  is  altered,  so  that  the  attraction  is  no  longer  radial  but  forwards,  and  the 
tangential  part  of  the  force  is 

R  dp  ,  dp 

—  =   +  fx  ->—>  forwards 001)- 


ira\a  —  a')       ds'  ds' 

4th.      By  the  tangential  displacement  of  the  second  ring  a  tangential  force  arises,  depending 
on  the  relation  between  the  length  of  the  waves  and  the  distance  between  the  rings. 
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a  -  a  .         />+  °°  «■»  sin  px  d.v 

If  we  make  m- — ; —  =  p,  and  m ^rj-  =  Hi 

a  J-n    (1  +  or)* 

the  tangential  force  is 

■  „      ;.anff=i>v (102). 

7r«  (a  -  a  )■* 

We  may  now  write  down  the  values  of  X,  n  and  v  by  transposing  accented  and  unaccented 
letters. 

R       (2a -a)  i?  R  n  fin*\ 

X  = —     v ix  = — "  = — ; r;D'     •     •      (103). 

Tra       (a  -  a')2  7ra(a  -  a)  ■7ra(a  -  a) 

Comparing  these  values  with  those  of  \'  fx'  and  v,  it  will  be  seen  that  the  following  relations 
are  approximately  true  when  a  is  nearly  equal  to  a  : 

X'  /u       v       R  a 

\  fj.       v      Ra 


(104). 


27.     To  form  the  equation  of  motion. 
The  original  equations  were 

w*  +  o»V  +  zw d4-  -  5r  =  p  -  s  +  K  -  &s  -  L>  AP  -  MBP  +  x'p  -*'%> 

r  dt       at  as 

2w  —  +  — -„  =  Q  =  MA<s  +  NBa  +  ^  ~rr  +  v'a. 
dt      dt-  ds 

Putting  p  =  A  cos  (ms  +  nt),    a  =  B  sin  (ms  +  nt), 

p'  =  A'  cos  (ms  +  nt),  a  =  B'  sin  (ms  +  nt), 
then  o>"  =  S  +  K 

(a.2  +  2S  +  n"  -  Z.)  A  +  (2wn  +  M)  B  -  X'A'  +  fx'mB'  =  0  1 

(2wn  +  M)  A  +  (n2  +  N)  B  -  /x'mA'  +  v'B'  =  0  j  (      '' 

The  corresponding  equations  for  the  second  ring  may  be  found  by  transposing  accented 
and  unaccented  letters.  We  should  then  have  four  equations  to  determine  the  ratios  of  JB, 
A'B',  and  a  resultant  equation  of  the  eighth  degree  to  determine  n.  But  we  may  make  use 
of  a  more  convenient  method,  since  X',  //,  and  v  are  small.      Eliminating  B  we  find 

An*  -  A(w"  +  2K  +  L-  N)  n2  -  iAwMn  +  AN  (3a,2)  1  _ 

(-  X'A'  +  ix'mB')  n2  +  (ix'mA'  -  v'B')2wn  j  '' 

Putting  B  =  (3A,  A'  =  w  A,  B"  =  (5' A'  =  (Z'xA, 

n*  -  \w2  (+  2K)  +  L  -  iV\  n2  -  \wMn  +  3w2N  )  _ 

+  (-  X'  +  fx'mfi')  ri*w  +  (fx'm  -  i//3')  2wnx  J  '' 

dU 

— -  =  4w3  -  2wsn  +  &c (108), 

dn 

dU 

—  =  -  X  n2  +  m  mfi'n*  +  2(x  main  -  2v  /3'o)« (109), 


we  have 
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,  dn      X'n  -  ii'mfi'n  -  Qn'mu  +  2i//3'w  .      . 

whence  7-= -„ (110)- 

dx  4?r  —  2a> 

28.  If  we  were  to  solve  the  equation  for  n,  leaving  out  the  terms  involving  x,,  we  should 
find  the  wave-velocities  of  the  four  free  waves  of  the  first  ring,  supposing  the  second  ring  to 
be  prevented  from  being  disturbed.  But  in  reality  the  waves  in  the  first  ring  produce  a 
disturbance  in  the  second,  and  these  in  turn  react  upon  the  first  ring,  so  that  the  wave-velocity 
is  somewhat  different  from  that  which  it  would  be  in  the  supposed  case.  Now  if  x  be  the 
ratio  of  the  radial  amplitude  of  displacement  in  the  second  ring  to  that  in  the  first,  and  if  n 
be  a  value  of  n  supposing  x  =  0,  then  by  Maclaurin's  theorem, 

—     dn 
n  =  n  +-—  x  .      ...........     (111). 

ax 

The   wave-velocity  relative  to  the   ring  is ,  and  the  absolute  angular  velocity  of  the 

m 

wave  in  space  is 

n  n        1  dn 

•ar  =  w =  u) —x (112), 

m  m      mdx 


=  p-qx (113), 


,  n  l  dn 

where  p  =  w ,   and  q  = . 

m  m  dx 


Similarly  in  the  second  ring  we  should  have 

t*'  =  P  -q  - (114); 

x 

and  since  the  corresponding  waves  in  the  two  rings  must  have  the  same  absolute  angular 
velocity, 

'  -  1 

•sr=sr,  or  p-qx=p  -a- (115). 

x 

This  is  a  quadratic  equation  in  x,  the  roots  of  which  are  real  when 

(p  -  p'y  +  499' 

is  positive.  When  this  condition  is  not  fulfilled,  the  roots  are  impossible,  and  the  general 
solution  of  the  equations  of  motion  will  contain  exponential  factors,  indicating  destructive 
oscillations  in  the  rings. 

Since  q  and  q  are  small  quantities,  the  solution  is  always  real  whenever  p  and  p  are 
considerably  different.  The  absolute  angular  velocities  of  the  two  pairs  of  reacting  waves,  are 
then  nearly 

w'         1   ,       qq 

p  +  ,  ,   and  p   - 


p-p  p-p 

instead  of  p  and  p,  as  they  would  have  been  if  there  had  been  no  reaction  of  the  forced  wave 
upon  the  free  wave  which  produces  it. 

H 
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When  p  and  p  are  equal  or  nearly  equal,  the  character  of  the  solution  will  depend  on  the 
sign  of  qq.      We  must  therefore  determine  the  signs  of  q  and  q  in  such  cases. 

Putting  /3  =      ,  - ,  we  may  write  the  values  of  q  and  q 


. ,  ,      /to       o>\        .  ,  to  ft)     -\ 

X  +  2/x m  {—;  —  -     -  4i>  — — 


4*1"  -  2  ft)2 


9  = 


/ft)         ft)  \  ,  0)    w 

,    X  +  2/xm ,      -  4k  —  - 

^  \n      n  I  n  n 


(116). 


4ra'a  -  2a)/z  ^ 

Referring  to  the  values  of  the  disturbing  forces,  we  find  that 

X'  f/       v       R'a 

X  i>.       v      Ra 

„  an  W  -  2tt>'2    iZ'a 

Hence  — ,  =  — . 017} 

q       ri  4rr  -  2o)2     Ra  v       ' 

Since  qq'  is  of  the  same  sign  as  ■—,  we  have  only  to  determine  whether  2n ,  and  2n 

q  n  n 

are  of  the  same  or  of  different  signs.      If  these  quantities  are  of  the  same  sign,  qq'  is  positive, 

if  of  different  signs,  qq'  is  negative. 

2 

Now  there  are  four  values  of  n,  which  give  four  corresponding  values  of  2w : 

n 

nt  =  —  o)  +  &c,  2w, is  negative, 

2 

w2  =  —  a  small  quantity,  2na is  positive, 

n, 

ii  •  w~  • 

n3  =  +  a  small  quantity,  2% is  negative, 

n3 

o?  . 

n4  =  ft)  -  &c,  2n4 is  positive. 

w4 

The  quantity  with  which  we  have  to  do  is  therefore  positive  for  the  even  orders  of  waves 
and  negative  for  the  odd  ones,  and  the  corresponding  quantity  in  the  other  ring  obeys  the  same 
law.  Hence  when  the  waves  which  act  upon  each  other  are  either  both  of  even  or  both  of  odd 
names,  qq  will  be  positive,  but  when  one  belongs  to  an  even  series,  and  the  other  to  an  odd 
series,  qq  is  negative. 

29.     The  values  of  p  and  p  are,  roughly, 

ft)  <a 

p1=w  +  -  -  &c,    p2  =  w  +  &c,  p3  =  w  -  &c.,      pi  =  w  -  —  +  &c, 

(118). 

,,tt>  ,,  ,  ,  ,  ,         G>  „ 

pt  =  (0   + Sec,  p2  =  <*>  +  &c-j  P3  =  w  —  &c->  p4  =<"-—+  &c. 
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w  is  Greater  than  a;,  so  that  jo,'  is  the  greatest,  and  pt  the  least  of  these  values,  and  of  those  of 
the  same  order,  the  accented  is  greater  than  the  unaccented.  The  following  cases  of  equality 
are  therefore  possible  under  suitable  circumstances : 

Pi  =  Ps,  Pi  =  Ps\ 

Pi  =  Pi,  Pi  =  Pi, 

Pt  =  Pi  (when  m  =  l),     p%  =  p3', 

Pz  =  Pi- 

In  the  cases  in  the  first  column  qq  will  be  positive,  in  those  in  the  second  column  qq  will 
be  negative. 

30.  Now  each  of  the  four  values  of  p  is  a  function  of  m,  the  number  of  undulations  in 
the  ring,  and  of  a  the  radius  of  the  ring,  varying  nearly  as  a"§.  Hence  m  being  given,  we 
may  alter  the  radius  of  the  ring  till  any  one  of  the  four  values  of  p  becomes  equal  to  a  given 
quantity,  say  a  given  value  of  p',  so  that  if  an  indefinite  number  of  rings  coexisted,  so  as  to 
form  a  sheet  of  rings,  it  would  be  always  possible  to  discover  instances  of  the  equality  of  p  and 
p'  among  them.  If  such  a  case  of  equality  belongs  to  the  first  column  given  above,  two 
constant  waves  will  arise  in  both  rings,  one  travelling  a  little  faster,  and  the  other  a  little  slower 
than  the  free  waves.  If  the  case  belongs  to  the  second  column,  two  waves  will  also  arise  in 
each  ring,  but  the  one  pair  will  gradually  die  away,  and  the  other  pair  will  increase  in 
amplitude  indefinitely,  the  one  wave  strengthening  the  other  till  at  last  both  rings  are  thrown 
into  confusion. 

The  only  way  in  which  such  an  occurrence  can  be  avoided  is  by  placing  the  rings  at  such 
a  distance  that  no  value  of  m  shall  give  coincident  values  of  p  and  p.  For  instance,  if 
m  >  2a),  but  u>'  <  So),  no  such  coincidence  is  possible.  For  p,  is  always  less  than  p2',  it  is 
greater  than  p±  when  m  =  1  or  2,  and  less  than  pt  when  m  is  3  or  a  greater  number.  There 
are  of  course  an  infinite  number  of  ways  in  which  this  noncoincidence  might  be  secured,  but  it 
is  plain  that  if  a  number  of  concentric  rings  were  placed  at  small  intervals  from  each  other, 
such  coincidences  must  occur  accurately  or  approximately  between  some  pairs  of  rings,  and  if 
the  value  of  (p  -  j))2  is  brought  lower  than  -  4-qq,  there  will  be  destructive  interference. 

This  investigation  is  applicable  to  any  number  of  concentric  rings,  for,  by  the  principle  of 
superposition  of  small  displacements,  the  reciprocal  actions  of  any  pair  of  rings  are  independent 
of  all  the  rest. 

31.  On  the  effect  of  long -continued  disturbances  on  a  system  of  rings. 

The  result  of  our  previous  investigations  has  been  to  point  out  several  ways  in  which 
disturbances  may  accumulate  till  collisions  of  the  different  particles  of  the  rings  take  place. 
After  such  a  collision  the  particles  will  still  continue  to  revolve  about  the  planet,  but  there  will 
be  a  loss  of  energy  in  the  system  during  the  collision  which  can  never  be  restored.  Such 
collisions  however  will  not  affect  what  is  called  the  Angular  Momentum  of  the  system  about 
the  planet,  which  will  therefore  remain  constant. 

h2 
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Let  M  be  the  mass  of  the  system  of  rings,  and  §m  that  of  one  ring  whose  radius  is  r,  and 
angular  velocity  w  =  S?  r"$.      The  angular  momentum  of  the  ring  is 

wr^lm  =  S^r^m, 
half  its  vis  viva  is  ^toVlm  =  I5r-15m. 

The  potential  energy  due  to  Saturn's  attraction  on  the  ring  is 

— -  dm  —  —  Sr~l$m. 
r 

The  angular  momentum  of  the  whole  system  is  invariable,  and  is 

Sf2(ri$m)  =  A (119). 

The  whole  energy  of  the  system  is  the  sum  of  half  the  vis  viva  and  the  potential  energy, 
and  is 

- %S2  (r-^ro)  -  E (120). 

A  is  invariable,  while  E  necessarily  diminishes.  We  shall  find  that  as  E  diminishes,  the 
distribution  of  the  rings  must  be  altered,  some  of  the  outer  rings  moving  outwards,  while 
the  inner  rings  move  inwards,  so  as  either  to  spread  out  the  whole  system  more,  both  on  the 
outer  and  on  the  inner  edge  of  the  system,  or,  without  affecting  the  extreme  rings,  to  diminish 
the  density  or  number  of  the  rings  at  the  mean  distance,  and  increase  it  at  or  near  the  inner 
and  outer  edges. 

Let  us  put  x  =  ri,  then  A  =  <$2  (xdm)  is  constant. 

2  (xdm) 
Now  let  x}  =  , 

2  (dm) 

and  x  =  «!  +  x, 

then  we  may  write 

2E 
-  -=  =  2  (r-^m)  =  2  (x~2  dm), 
S 

=  2dm  (a?,"*  -  2  —  +  3  —  -  &c), 

a;,3         .V 

=  —  2  {dm)  -  —. -  2  (xdm)  +  —d  2  (.r'^m)  -  &c.     .     ( 121). 
#,  at  1  x* 

Now  2  (  dm)  =  M  a  constant,  2  {xdm)  =  0,  and  2  (a/2<V)  is  a  quantity  which  increases 
when  the  rings  are  spread  out  from  the  mean  distance  either  way,  x'  being  subject  only  to  the 
restriction  2  {xdm)  =  0.  But  2  (x'2dm)  may  increase  without  the  extreme  values  of  x  being 
increased,  provided  some  other  values  be  increased. 

32.  In  fact,  if  we  consider  the  very  innermost  particle  as  moving  in  an  ellipse,  and  at  the 
further  apse  of  its  orbit  encountering  another  particle  belonging  to   a  larger  orbit,  we  know 
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that  the  second  particle,  when  at  the  same  distance  from  the  planet,  moves  the  faster.  The 
result  is,  that  the  interior  satellite  will  receive  a  forward  impulse  at  its  further  apse,  and  will 
move  in  a  larger  and  less  eccentric  orbit  than  before.  In  the  same  way  one  of  the  outermost 
particles  may  receive  a  backward  impulse  at  its  nearer  apse,  and  so  be  made  to  move  in  a 
smaller  and  less  eccentric  orbit  than  before.  When  we  come  to  deal  with  collisions  among 
bodies  of  unknown  number,  size,  and  shape,  we  can  no  longer  trace  the  mathematical  laws  of 
their  motion  with  any  distinctness.  All  we  can  now  do  is  to  collect  the  results  of  our  investi- 
gations and  to  make  the  best  use  we  can  of  them  in  forming  an  opinion  as  to  the  constitution 
of  the  actual  rings  of  Saturn  which  are  still  in  existence  and  apparently  in  stead)'  motion, 
whatever  catastrophes  may  be  indicated  by  the  various  theories  we  have  attempted. 

33.  To  find  the  Loss  of  Energy  due  to  internal  friction  in  a  broad  Fluid  Ring,  the 
parts  of  which  revolve  about  the  Planet,  each  with  the  velocity  of  a  satellite  at  the  same 
distance. 

Conceive  a  fluid,  the  particles  of  which  move  parallel  to  the  axis  of  x  with  a  velocity  u, 
m  being  a  function  of  sr,  then  there  will  be  a  tangential  pressure  on  a  plane  parallel  to  xy 

du 

=  fi  — —  on  unit  of  area 
dz 

due  to  the  relative  sliding  of  the  parts  of  the  fluid  over  each  other. 

In  the  case  of  the  ring  we  have 

w  =  S'-r-i. 

The  absolute  velocity  of  any  particle  is  wr.      That  of  a  particle  at  distance  (r  +  Sr)  is 

wr  +  —  (wr)  Sr. 
dr 

If  the  angular  velocity  had  been  uniform,  there  would  have  been  no  sliding,  and  the  velocity 
would  have  been 

wr  +  wSr. 

The  sliding  is  therefore 

dw  . 
dr 

,  .  ,  ,,     ,  .  dw 

and  the  friction  on  unit  of  area  perpendicular  tor  is  y}1  —  . 
r  dr 

The  loss  of  Energy,  per  unit  of  area,  is  the  product  of  the  sliding  by  the  friction, 

dw}2 
or  M,.s  —     Sr  in  unit  of  time. 

dr 
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The  loss  of  Energy  in  a  part  of  the  Ring  whose  radius  is  r,  breadth  3r,  and  thickness  c, 

2 
3r. 


,     dw 

Zirrcij.  — - 

dr 


g 
In  the  case  before  us  it  is  -  TraScr'^r. 

2 

If  the  thickness  of  the  ring  is  uniform  between  r  =  a  and  r  =  b,  the  whole  loss  of  Energy 

is  -irfxSci-  -  -  )  , 

2  \6      a) 

in  unit  of  time. 

Now  half  the  vis  viva  of  an  elementary  ring  is 

irpcrSr  r2w~  =  irpcSlr, 

and  this  between  the  limits  r  =  a  and  r  =  b  gives 

7rpcS(a  —  b). 

The  potential  due  to   the  attraction  of  .Sis  twice  this  quantity  with    the   sign  changed,  so 
that 

E  =  -  wpcS(a  -b). 

dE    9     _yi     i\ 

and  — -  =  -  tt/uS , 

dt       2  \b       a) 

l    dE  _       9  m    l 
£    d£  2  p  ab' 

A 
Now    Professor  Stokes  finds  \/  —  =  0'0564  for  water, 

P 

and  =  0-116    for  air, 

taking  the  unit  of  space  one  English  inch,  and  the  unit  of  time  one  second.  We  may  take 
a  =  88,209  miles,  and  b  =  77,636  for  the  ring  A;  and  a  =  75,845  and  b  =  58,660  for  the  ring  b. 
We  may  also  take  one  year  as  the  unit  of  time.  The  quantity  representing  the  ratio  of  the  loss 
of  energy  in  a  year  to  the  whole  energy  is 

1    dE  1 


E    dt        60,880,000,000,000 


for  the  ring  A, 


and  — ■ for  the  ring  B, 

39,540,000,000,000 


showing  that  the  effect  of  internal  friction  in  a  ring  of  water  moving  with  steady  motion  is 
inappreciably  small.  It  cannot  be  from  this  cause  therefore  that  any  decay  can  take  place  in 
the  motion  of  the  ring,  provided  that  no  waves  arise  to  disturb  the  motion. 
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Recapitulation  of  the  Theory  of  the  Motion  of  a  Rigid  Ring. 

The  position  of  the  ring  relative  to  Saturn  at  any  given  instant  is  defined  by  three  variable 
quantities. 

1st.  The  distance  between  the  centre  of  gravity  of  Saturn  and  the  centre  of  gravity  of 
the  ring.     This  distance  we  denote  by  r. 

2nd.  The  angle  which  the  line  r  makes  with  a  fixed  line  in  the  plane  of  the  motion  of 
the  ring.     This  angle  is  called  6. 

3rd.  The  angle  between  the  line  r  and  a  line  fixed  with  respect  to  the  ring  so  that  it 
coincides  with  r  when  the  ring  is  in  its  mean  position.      This  is  the  angle  (p. 

The  values  of  these  three  quantities  determine  the  position  of  the  ring  so  far  as  its  motion 
in  its  own  plane  is  concerned.  They  may  be  referred  to  as  the  radius  vector,  longitude,  and 
angle  of  libration  of  the  ring. 

The  forces  which  act  between  the  ring  and  the  planet  depend  entirely  upon  their 
relative  position.  The  method  adopted  above  consists  in  determining  the  potential  (V)  of 
the  ring  at  the  centre  of  the  planet  in  terms  of  r  and  (p.  Then  the  work  done  by  any 
displacement  of  the  system  is  measured  by  the  change  of  VS  during  that  displacement.     The 

dV 
attraction  between  the  centre  of  gravity  of  the  Ring  and  that  of  the  planet  is  -  S  — — ,  and 

dV 
the  moment  of   the  couple  tending  to  turn  the  ring  about  its  centre  of  gravity  is  S—r-  . 

d(p 

It  is  proved  in  Problem  V.  that  if  a  be  the  radius  of  a  circular  ring,  r0  =  af  the  distance 
of  its  centre  of  gravity  from  the  centre  of  the  circle,  and  R  the  mass  of  the  ring,  then,  at   the 

„  ,       .        dV  R  „  dV 

centre  of  the  ring,  —  =-—-/,  — —  =  o. 
8    dr  a*J    dcp 

t      ,  ,_       d°  V      ,  R ,  V       ,  •  ,    • 

It  also  appears  that  — —  =  £  —  (1  +  g),  which  is  positive  when  g  >  —  1, 
dr~  a 

and  that  — -  =  A  —  f2  (3  —  g),  which  is  positive  when  g<3. 
dtp  a 

d?V 
If  — —   is  positive,  then    the    attraction    between    the   centres  decreases    as    the   distance 
dr2         * 

increases,  so  that,  if  the  two  centres  were  kept  at  rest  at  a  given  distance  by  a  constant  force, 

d2V 
the  equilibrium  would  be  unstable.      If  — :   is  positive,  then  the  forces   tend   to  increase  the 

angle  of  libration,  in  whichever  direction  the  libration  takes  place,  so  that  if  the  ring  were 
fixed  by  an  axis  through  its  centre  of  gravity,  its  equilibrium  round  that  axis  would  be 
unstable. 

In  the  case  of  the  uniform  ring  with  a  heavy  particle  on  its  circumference  whose  weight 
=  "82  of  the  whole,  the  direction  of  the  whole  attractive  force  of  the  ring  near  the  centre  will 
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pass  through  a  point  lying  in  the  same  radius  as  the  centre  of  gravity,  but  at  a  distance   from 
the  centre  =  f  «•      (Fig.  6.) 

If  we  call  this  point  0,  the  line  SO  will  indicate  the  direction  and  position  of  the  force 
acting  on  the  ring,  which  we  may  call  F. 

It  is  evident  that  the  force  F,  acting  on  the  ring  in  the  line  OS,  will  tend  to  turn  it 
round  its  centre  of  gravity  R  and  to  increase  the  angle  of  libration  KRO.  The  direct  action 
of  this  force  can  never  reduce  the  angle  of  libration  to  zero  again.  To  understand  the  indirect 
action  of  the  force,  we  must  recollect  that  the  centre  of  gravity  (R)  of  the  ring  is  revolving 
about  Saturn  in  the  direction  of  the  arrows,  and  that  the  ring  is  revolving  about  its  centre  of 
gravity  with  nearly  the  same  velocity.  If  the  angular  velocity  of  the  centre  of  gravity  about 
Saturn  were  always  equal  to  the  rotatory  velocity  of  the  ring,  there  would  be  no  libration. 

Now  suppose  that  the  angle  of  rotation  of  the  ring  is  in  advance  of  the  longitude  of  its 
centre  of  gravity,  so  that  the  line  RO  has  got  in  advance  of  SRK  by  the  angle  of  libration 
KRO.  The  attraction  between  the  planet  and  the  ring  is  a  force  F  acting  in  SO.  We 
resolve  this  force  into  a  couple,  whose  moment  is  FRN,  and  a  force  F  acting  through  R  the 
centre  of  gravity  of  the  ring. 

The  couple  affects  the  rotation  of  the  ring,  but  not  the  position  of  its  centre  of  gravity, 
and  the  force  RF  acts  on  the  centre  of  gravity  without  affecting  the  rotation. 

Now  the  couple,  in  the  case  represented  in  the  figure,  acts  in  the  positive  direction,  so  as 
to  increase  the  angular  velocity  of  the  ring,  which  was  already  greater  than  the  velocity  of 
revolution  of  R  about  S,  so  that  the  angle  of  libration  would  increase,  and  never  be  reduced 
to  zero. 

The  force  RF  does  not  act  in  the  direction  of  S,  but  behind  it,  so  that  it  becomes  a 
retarding  force  acting  upon  the  centre  of  gravity  of  the  ring.  Now  the  effect  of  a  retarding 
force  is  to  cause  the  distance  of  the  revolving  body  to  decrease  and  the  angular  velocity 
to  increase,  so  that  a  retarding  force  increases  the  angular  velocity  of  R  about  S. 

The  effect  of  the  attraction  along  SO  in  the  case  of  the  figure  is,  first,  to  increase  the  rate 
of  rotation  of  the  ring  round  R,  and  secondly,  to  increase  the  angular  velocity  of  R  about  S. 
If  the  second  effect  is  greater  than  the  first,  then,  although  the  line  RO  increases  its  angular 
velocity,  SR  will  increase  its  angular  velocity  more,  and  will  overtake  RO,  and  restore  the  ring 
to  its  original  position,  so  that  SRO  will  be  made  a  straight  line  as  at  first.  If  this  accelerat- 
ing effect  is  not  greater  than  the  acceleration  of  rotation  about  R  due  to  the  couple,  then  no 
compensation  will  take  place,  and  the  motion  will  be  essentially  unstable. 

If  in  the  figure  we  had  drawn  (p  negative  instead  of  positive,  then  the  couple  would  have 
been  negative,  the  tangential  force  on  R  accelerative,  r  would  have  increased,  and  in  the  cases 
of  stability  the  retardation  of  9  would  be  greater  than  that  of  (9  +  (f>),  and  the  normal  position 
would  be  restored,  as  before. 

The  object  of  the  investigation  is  to  find  the  conditions  under  which  this  compensation 
is  possible. 
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It  is  evident  that  when  SRO  becomes  straight,  there  is  still  a  difference  of  angular 
velocities  between  the  rotation  of  the  ring  and  the  revolution  of  the  centre  of  gravity,  so  that 
there  will  be  an  oscillation  on  the  other  side,  and  the  motion  will  proceed  by  alternate 
oscillations  without  limit. 

If  we  begin  with  r  at  its  mean  value,  and  cf>  negative,  then  the  rotation  of  the  ring  will 
be  retarded,  r  will  be  increased,  the  revolution  of  r  will  be  more  retarded,  and  thus  rt>  will  be 
reduced  to  zero.  The  next  part  of  the  motion  will  reduce  r  to  its  mean  value,  and  bring  cp 
to  its  greatest  positive  value.  Then  r  will  diminish  to  its  least  value,  and  <p  will  vanish. 
Lastly  r  will  return  to  the  mean  value,  and  <p  to  the  greatest  negative  value. 

It  appears  from  the  calculations,  that  there  are,  in  general,  two  different  ways  in  which 
this  kind  of  motion  may  take  place,  and  that  these  may  have  different  periods,  phases,  and 
amplitudes.  The  mental  exertion  required  in  following  out  the  results  of  a  combined  motion 
of  this  kind,  with  all  the  variations  of  force  and  velocity  during  a  complete  cycle,  would  be 
very  great  in  proportion  to  the  additional  knowledge  we  should  derive  from  the  exercise. 

The  result  of  this  theory  of  a  rigid  ring  shows  not  only  that  a  perfectly  uniform  ring 
cannot  revolve  permanently  about  the  planet,  but  that  the  irregularity  of  a  permanently 
revolving  ring  must  be  a  very  observable  quantity,  the  distance  between  the  centre  of  the  ring 
and  the  centre  of  gravity  being  between  -8158  and  '8279  of  the  radius.  As  there  is  no 
appearance  about  the  rings  justifying  a  belief  in  so  great  an  irregularity,  the  theory  of  the 
solidity  of  the  rings  becomes  very  improbable. 

When  we  come  to  consider  the  additional  difficulty  of  the  tendency  of  the  fluid  or  loose 
parts  of  the  ring  to  accumulate  at  the  thicker  parts,  and  thus  to  destroy  that  nice  adjustment 
of  the  load  on  which  stability  depends,  we  have  another  powerful  argument  against  solidity. 

And  when  we  consider  the  immense  size  of  the  rings,  and  their  comparative  thinness,  the 
absurdity  of  treating  them  as  rigid  bodies  becomes  self-evident.  An  iron  ring  of  such  a  size 
would  be  not  only  plastic  but  semifluid  under  the  forces  which  it  would  experience,  and  we 
have  no  reason  to  believe  these  rings  to  be  artificially  strengthened  with  any  material  unknown 
on  this  earth. 

Recapitulation  of  the    Theory  of  a  Ring  of  equal  Satellites. 

In  attempting  to  conceive  of  the  disturbed  motion  of  a  ring  of  unconnected  satellites,  we 
have,  in  the  first  place,  to  devise  a  method  of  identifying  each  satellite  at  any  given  time,  and 
in  the  second  place,  to  express  the  motion  of  every  satellite  under  the  same  general  formula, 
in  order  that  the  mathematical  methods  may  embrace  the  whole  system  of  bodies  at  once. 

By  conceiving  the  ring  of  satellites  arranged  regularly  in  a  circle,  we  may  easily  identify 
any  satellite,  by  stating  the  angular  distance  between  it  and  a  known  satellite  when  so 
arranged.  If  the  motion  of  the  ring  were  undisturbed,  this  angle  would  remain  unchanged  during 
the  motion,  but  in  reality,  the  satellite  has  its  position  altered  in  three  ways  :  1st,  it  may  be 
further  from  or  nearer  to  Saturn;  2ndly,  it  may  be  in  advance  or  in  the  rear  of  the  position  it 
would  have  had  if  undisturbed;  3rdly,  it  may  be  on  one  side  or  other  of  the  mean  plane  of  the 

I 


58  ON   THE   STABILITY    OF   THE   MOTION 

ring.  Each  of  these  displacements  may  vary  in  any  way  whatever  as  we  pass  from  one  satellite 
to  another,  so  that  it  is  impossible  to  assign  beforehand  the  place  of  any  satellite  by  knowing 
the  places  of  the  rest.      ^  2. 

The  formula,  therefore,  by  which  we  are  enabled  to  predict  the  place  of  every  satellite  at 
any  given  time,  must  be  such  as  to  allow  the  initial  position  of  every  satellite  to  be  independent 
of  the  rest,  and  must  express  all  future  positions  of  that  satellite  by  inserting  the  corresponding 
value  of  the  quantity  denoting  time,  and  those  of  every  other  satellite  by  inserting  the  value 
of  the  angular  distance  of  the  given  satellite  from  the  point  of  reference.  The  three 
displacements  of  the  satellite  will  therefore  be  functions  of  two  variables — the  angular  position 
of  the  satellite,  and  the  time.  When  the  time  alone  is  made  to  vary,  we  trace  the  complete 
motion  of  a  single  satellite;  and  when  the  time  is  made  constant,  and  the  angle  is  made  to 
vary,  we  trace  the  form  of  the  ring  at  a  given  time. 

It  is  evident  that  the  form  of  this  function,  in  so  far  as  it  indicates  the  state  of  the  whole 
ring  at  a  given  instant,  must  be  wholly  arbitrary,  for  the  form  of  the  ring  and  its  motion  at 
starting  are  limited  only  by  the  condition  that  the  irregularities  must  be  small.  We  have, 
however,  the  means  of  breaking  up  any  function,  however  complicated,  into  a  series  of  simple 
functions,  so  that  the  value  of  the  function  between  certain  limits  may  be  accurately  expressed 
as  the  sum  of  a  series  of  sines  and  cosines  of  multiples  of  the  variable.  This  method,  due 
to  Fourier,  is  peculiarly  applicable  to  the  case  of  a  ring  returning  into  itself,  for  the  value  of 
Fourier's  series  is  necessarily  periodic.  We  now  regard  the  form  of  the  disturbed  ring  at  any 
instant  as  the  result  of  the  superposition  of  a  number  of  separate  disturbances,  each  of  which 
is  of  the  nature  of  a  series  of  equal  waves  regularly  arranged  round  the  ring.  Each  of  these 
elementary  disturbances  is  characterised  by  the  number  of  undulations  in  it,  by  their  amplitude, 
and  by  the  position  of  the  first  maximum  in  the  ring.      Q  3. 

When  we  know  the  form  of  each  elementary  disturbance,  we  may  calculate  the  attraction  of 
the  disturbed  ring  on  any  given  particle  in  terms  of  the  constants  belonging  to  that  disturbance,  so 
that  as  the  actual  displacement  is  the  resultant  of  the  elementary  displacements,  the  actual  attrac- 
tion will  be  the  resultant  of  the  corresponding  elementary  attractions,  and  therefore  the  actual 
motion  will  be  the  resultant  of  all  the  motions  arising  from  the  elementary  disturbances.  We  have 
therefore  only  to  investigate  the  elementary  disturbances  one  by  one,  and  having  established  the 
theory  of  these,  we  calculate  the  actual  motion  by  combining  the  series  of  motions  so  obtained. 

Assuming  the  motion  of  the  satellites  in  one  of  the  elementary  disturbances  to  be  that  of 
oscillation  about  a  mean  position,  and  the  whole  motion  to  be  that  of  a  uniformly  revolving 
series  of  undulations,  we  find  our  supposition  to  be  correct,  provided  a  certain  biquadratic 
equation  is  satisfied  by  the  quantity  denoting  the  rate  of  oscillation.      ^  6. 

When  the  four  roots  of  this  equation  are  all  real,  the  motion  of  each  satellite  is  compounded 
of  four  different  oscillations  of  different  amplitudes  and  periods,  and  the  motion  of  the  whole  ring 
consists  of  four  series  of  undulations,  travelling  round  the  ring  with  different  velocities.  When 
any  of  these  roots  are  impossible,  the  motion  is  no  longer  oscillatory,  but  tends  to  the  rapid 
destruction  of  the  ring. 
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To  determine  whether  the  motion  of  the  ring  is  permanent,  we  must  assure  ourselves  that 
the  four  roots  of  this  equation  are  real,  whatever  be  the  number  of  undulations  in  the  ring; 
for  if  any  one  of  the  possible  elementary  disturbances  should  lead  to  destructive  oscillations,  that 
disturbance  might  sooner  or  later  commence,  and  the  ring  would  be  destroyed. 

Now  the  number  of  undulations  in  the  ring  may  be  any  whole  number  from  one  up  to  half 
the  number  of  satellites.  The  forces  from  which  danger  is  to  be  apprehended  are  greatest  when 
the  number  of  undulations  is  greatest,  and  by  taking  that  number  equal  to  half  the  number  of 
satellites,  we  find  the  condition  of  stability  to  be 

S  >  .  4352  (SR, 
where  S  is  the  mass  of  the  central  body,  R  that  of  the  ring,  and  n  the  number  of  satellites  of 
which  it  is  composed.  $  8.  If  the  number  of  satellites  be  too  great,  destructive  oscillations  will 
commence,  and  finally  some  of  the  satellites  will  come  into  collision  with  each  other  and  unite, 
so  that  the  number  of  independent  satellites  will  be  reduced  to  that  which  the  central  body 
can  retain  and  keep  in  discipline.  When  this  has  taken  place,  the  satellites  will  not  only  be  kept 
at  the  proper  distance  from  the  primary,  but  will  be  prevented  by  its  preponderating  mass  from 
interfering  with  each  other. 

We  next  considered  more  carefully  the  case  in  which  the  mass  of  the  ring  is  very  small,  so 
that  the  forces  arising  from  the  attraction  of  the  ring  are  small  compared  with  that  due  to  the 
central  body.  In  this  case  the  values  of  the  roots  of  the  biquadratic  are  all  real,  and  easilv 
estimated.      fi  9- 

If  we  consider  the  motion  of  any  satellite  about  its  mean  position,  as  referred  to  axes 
fixed  in  the  plane  of  the  ring,  we  shall  find  that  it  describes  an  ellipse  in  the  direction  opposite 
to  that  of  the  revolution  of  the  ring,  the  periodic  time  being  to  that  of  the  ring  as  to  to  n,  and 
the  tangential  amplitude  of  oscillation  being  to  the  radial  as  2w  to  n.     §  10. 

The  absolute  motion  of  each  satellite  in  space  is  nearly  elliptic  for  the  large  values  of  n,  the 
axis  of  the  ellipse  always  advancing  slowly  in  the  direction  of  rotation.  The  path  of  a  satellite 
corresponding  to  one  of  the  small  values  of  n  is  nearly  circular,  but  the  radius  slowly  increases 
and  diminishes  during  a  period  of  many  revolutions,      fi  11. 

The  form  of  the  ring  at  any  instant  is  that  of  a  re-entering  curve,  having  m  alternations  of 
distance  from  the  centre,  symmetrically  arranged,  and  m  points  of  condensation,  or  crowding  of 
the  satellites,  which  coincide  with  the  points  of  greatest  distance  when  n  is  positive,  and  with  the 
points  nearest  the  centre  when  n  is  negative.      5  12. 

This  system   of  undulations  travels  with  an  angular  velocity relative  to  the  ring,  and 

m 
n   . 
<v  -  —  in  space,  so  that  during  each  oscillation  of  a  satellite  a  complete  wave  passes  over  it.  fi  14. 

To  exhibit  the  movements  of  the  satellites,  I  have  made  an  arrangement  by  which  36  little 
ivory  balls  are  made  to  go  through  the  motions  belonging  to  the  first  or  fourth  series  of  waves. 
(Figs.  7,  8.) 

The  instrument  stands  on  a  pillar  A,  in  the  upper  part  of  which  turnsthe  cranked  axle  CC. 

i  2 
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On  the  parallel  parts  of  this  axle  are  placed  two  wheels,  RR  and  TT,  each  of  which  has  36  holes 
at  equal  distances  in  a  circle  near  its  circumference.  The  two  circles  are  connected  by  36  small 
cranks  of  the  form  KK,  the  extremities  of  which  turn  in  the  corresponding  holes  of  the  two 
wheels.  That  axle  of  the  crank  K  which  passes  through  the  hole  in  the  wheel  S  is  bored,  so  as 
to  hold  the  end  of  the  bent  wire  which  carries  the  satellite  S.  This  wire  may  be  turned  in  the  hole 
so  as  to  place  the  bent  part  carrying  the  satellite  at  any  angle  with  the  crank.  A  pin  P,  which 
passes  through  the  top  of  the  pillar,  serves  to  prevent  the  cranked  axle  from  turning;  and  a  pin 
Q,  passing  through  the  pillar  horizontally,  may  be  made  to  fix  the  wheel  R,  by  inserting  it  in  a 
hole  in  one  of  the  spokes  of  that  wheel.  There  is  also  a  handle  H,  which  is  in  one  piece  with 
the  wheel  T,  and  serves  to  turn  the  axle. 

Now  suppose  the  pin  P  taken  out,  so  as  to  allow  the  cranked  axle  to  turn,  and  the  pin  Q 
inserted  in  its  hole,  so  as  to  prevent  the  wheel  R  from  revolving  ;  then  if  the  crank  C  be  turned 
by  means  of  the  handle  H,  the  wheel  T  will  have  its  centre  carried  round  in  a  vertical  circle,  but 
will  remain  parallel  to  itself  during  the  whole  motion,  so  that  every  point  in  its  plane  will  describe 
an  equal  circle,  and  all  the  cranks  K  will  be  made  to  revolve  exactly  as  the  large  crank  C  does. 
Each  satellite  will  therefore  revolve  in  a  small  circular  orbit,  in  the  same  time  with  the  handle  H, 
but  the  position  of  each  satellite  in  that  orbit  may  be  arranged  as  we  please,  according  as  we  turn 
the  wire  which  supports  it  in  the  end  of  the  crank. 

In  fig  8,  which  gives  a  front  view  of  the  instrument,  the  satellites  are  so  placed  that  each  is 
turned  60°  further  round  in  its  socket  than  the  one  behind  it.  As  there  are  $6  satellites,  this 
process  will  bring  us  back  to  our  starting-point  after  6  revolutions  of  the  direction  of  the  arm  of 
the  satellite  ;  and  therefore  as  we  have  gone  round  the  ring  once  in  the  same  direction,  the  arm 
of  the  satellite  will  have  overtaken  the  radius  of  the  ring  five  times. 

Hence  there  will  be  five  places  where  the  satellites  are  beyond  their  mean  distance  from  the 
centre  of  the  ring,  and  five  where  they  are  within  it,  so  that  we  have  here  a  series  of  five 
undulations  round  the  circumference  of  the  ring.  In  this  case  the  satellites  are  crowded  together 
when  nearest  to  the  centre,  so  that  the  case  is  that  of  the  first  series  of  waves,  when  m  =  5. 

Now  suppose  the  cranked  axle  C  to  be  turned,  and  all  the  small  cranks  K  to  turn  with  it, 
as  before  explained,  every  satellite  will  then  be  carried  round  on  its  own  arm  in  the  same 
direction  ;  but,  since  the  direction  of  the  arms  of  different  satellites  is  different,  their  phases  of 
revolution  will  preserve  the  same  difference,  and  the  system  of  satellites  will  still  be  arranged  in 
five  undulations,  only  the  undulations  will  be  propagated  round  the  ring  in  the  direction 
opposite  to  that  of  the  revolution  of  the  satellites. 

To  understand  the  motion  better,  let  us  conceive  the  centres  of  the  orbits  of  the  satellites  to 
be  arranged  in  a  straight  line  instead  of  a  circle,  as  in  fig  10.  Each  satellite  is  here  represented 
in  a  different  phase  of  its  orbit,  so  that  as  we  pass  from  one  to  anoth  er  from  left  to  right,  we  find 
the  position  of  the  satellite  in  its  orbit  altering  in  the  direction  opposite  to  that  of  the  hands  of 
a  watch.  The  satellites  all  lie  in  a  trochoidal  curve,  indicated  by  the  line  through  them  in  the 
figure.  Now  conceive  every  satellite  to  move  in  its  orbit  through  a  certain  angle  in  the  direction 
of  the  arrows.      The-  satellites  will  then  lie  in  the  dotted  line,  the  form  of  which  is  the  same  as 
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that  of  the  former  curve,  only  shifted  in  the  direction  of  the  large  arrow.  It  appears,  therefore, 
that  as  the  satellites  revolve,  the  undulation  travels,  so  that  any  part  of  it  reaches  successively 
each  satellite  as  it  comes  into  the  same  phase  of  rotation.  It  therefore  travels  from  those  satellites 
which  are  most  advanced  in  phase  to  those  which  are  less  so,  and  passes  over  a  complete  wave- 
length in  the  time  of  one  revolution  of  a  satellite. 

Now  if  the  satellites  be  arranged  as  in  fig.  8,  where  each  is  more  advanced  in  phase  as  we 
go  round  the  ring  in  the  direction  of  rotation,  the  wave  will  travel  in  the  direction  opposite  to 
that  of  rotation,  but  if  they  are  arranged  as  in  fig.  12,  where  each  satellite  is  less  advanced  in 
phase  as  we  go  round  the  ring,  the  wave  will  travel  in  the  direction  of  rotation.  Fig.  8 
represents  the  first  series  of  waves  where  m  =  5,  and  fig.  12  represents  the  fourth  series  where 
m  =  7.  By  arranging  the  satellites  in  their  sockets  before  starting,  we  might  make  m  equal 
to  any  whole  number,  from  1  to  18.  If  we  chose  any  number  above  18  the  result  would  be 
the  same  as  if  we  had  taken  a  number  as  much  below  1 8  and  changed  the  arrangement  from 
the  first  wave  to  the  fourth. 

In  this  wav  we  can  exhibit  the  motions  of  the  satellites  in  the  first  and  fourth  waves.  In 
reality  they  ought  to  move  in  ellipses,  the  major  axes  being  twice  the  minor,  whereas  in  the 
machine  they  move  in  circles :  but  the  character  of  the  motion  is  the  same,  though  the  form  of 
the  orbit  is  different. 

We  may  now  show  these  motions  of  the  satellites  among  each  other,  combined  with  the 
motion  of  rotation  of  the  whole  ring.  For  this  purpose  we  put  in  the  pin  P,  so  as  to  prevent 
the  crank  axle  from  turning,  and  take  out  the  pin  Q  so  as  to  allow  the  wheel  R  to  turn.  If 
we  then  turn  the  wheel  T,  all  the  small  cranks  will  remain  parallel  to  the  fixed  crank,  and  the 
wheel  R  will  revolve  at  the  same  rate  as  T.  The  arm  of  each  satellite  will  continue  parallel 
to  itself  during  the  motion,  so  that  the  satellite  will  describe  a  circle  whose  centre  is  at  a 
distance  from  the  centre  of  R,  equal  to  the  arm  of  the  satellite,  and  measured  in  the  same 
direction.  In  our  theory  of  real  satellites,  each  moves  in  an  ellipse,  having  the  central  body 
in  its  focus,  but  this  motion  in  an  eccentric  circle  is  sufficiently  near  for  illustration.  The 
motion  of  the  waves  relative  to  the  ring  is  the  same  as  before.  The  waves  of  the  first  kind 
travel  faster  than  the  ring  itself,  and  overtake  the  satellites,  those  of  the  fourth  kind  travel 
slower,  and  are  overtaken  by  them. 

In  fig.  11  we  have  an  exaggerated  representation  of  a  ring  of  twelve  satellites  affected  by  a 
wave  of  the  fourth  kind  where  m  =  2.  The  satellites  here  lie  in  an  ellipse  at  any  given  instant, 
and  as  each  moves  round  in  its  circle  about  its  mean  position,  the  ellipse  also  moves  round  in 
the  same  direction  with  half  their  angular  velocity.  In  the  figure  the  dotted  line  represents 
the  position  of  the  ellipse  when  each  satellite  has  moved  forward  into  the  position  represented 
by  a  dot. 

Fig.  13  represents  a  wave  of  the  first  kind  where  m  =  2.  The  satellites  at  any  instant  lie 
in  an  epitrochoid,  which,  as  the  satellites  revolve  about  their  mean  positions,  revolves  in  the 
opposite  direction  with  half  their  angular  velocity,  so  that  when  the  satellites  come  into  the 
positions  represented  by  the  dots,  the  curve  in  which  they  lie  turns  round  in  the  opposite 
direction  and  forms  the  dotted  curve. 
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In  fig.  o  we  have  the  same  case  as  in  fig.  13,  only  that  the  absolute  orbits  of  the  satellites 
in  space  are  given,  instead  of  their  orbits  about  their  mean  positions  in  the  ring.  Here  each 
moves  about  the  central  body  in  an  eccentric  circle,  which  in  strictness  ought  to  be  an  ellipse 
not  differing  much  from  the  circle. 

As  the  satellites  move  in  their  orbits  in  the  direction  of  the  arrows,  the  curve  which  they 
form  revolves  in  the  same  direction  with  a  velocity  iX  times  that  of  the  ring. 

By  considering  these  figures,  and  still  more  by  watching  the  actual  motion  of  the  ivory 
balls  in  the  model,  we  may  form  a  distinct  notion  of  the  motions  of  the  particles  of  a  discon- 
tinuous ring,  although  the  motions  of  the  model  are  circular  and  not  elliptic.  The  model, 
represented  on  a  scale  of  one-third  in  figs.  7  and  8,  was  made  in  brass  by  Messrs.  Smith  and 
Ramage  of  Aberdeen. 

We  are  now  able  to  understand  the  mechanical  principle,  on  account  of  which  a  massive 
central  body  is  enabled  to  govern  a  numerous  assemblage  of  satellites,  and  to  space  them  out 
into  a  regular  ring;  while  a  smaller  central  body  would  allow  disturbances  to  arise  among  the 
individual  satellites,  and  collisions  to  take  place. 

When  we  calculated  the  attractions  among  the  satellites  composing  the  ring,  we  found  that 
if  any  satellite  be  displaced  tangentially,  the  resultant  attraction  will  draw  it  away  from  its 
mean  position,  for  the  attraction  of  the  satellites  it  approaches  will  increase,  while  that  of 
those  it  recedes  from  will  diminish,  so  that  its  equilibrium  when  in  the  mean  position  is 
unstable  with  respect  to  tangential  displacements ;  and  therefore,  since  every  satellite  of  the  ring 
is  statically  unstable  between  its  neighbours,  the  slightest  disturbance  would  tend  to  produce 
collisions  among  the  satellites,  and  to  break  up  the  ring  into  groups  of  conglomerated  satellites. 

But  if  we  consider  the  dynamics  of  the  problem,  we  shall  find  that  this  effect  need  not 
necessarily  take  place,  and  that  this  very  force  which  tends  towards  destruction  may  become 
the  condition  of  the  preservation  of  the  ring.  Suppose  the  whole  ring  to  be  revolving  round 
a  central  body,  and  that  one  satellite  gets  in  advance  of  its  mean  position.  It  will  then  be 
attracted  forwards,  its  path  will  become  less  concave  towards  the  attracting  body,  so  that  its 
distance  from  that  body  will  increase.  At  this  increased  distance  its  angular  velocity  will  be 
less,  so  that  instead  of  overtaking  those  in  front,  it  may  by  this  means  be  made  to  fall  back- 
to  its  original  position.  Whether  it  does  so  or  not  must  depend  on  the  actual  values  of  the 
attractive  forces  and  on  the  angular  velocity  of  the  ring.  When  the  angular  velocity  is  great 
and  the  attractive  forces  small,  the  compensating  process  will  go  on  vigorously,  and  the  ring 
will  be  preserved.  When  the  angular  velocity  is  small  and  the  attractive  forces  of  the  ring 
great,  the  dynamical  effect  will  not  compensate  for  the  disturbing  action  of  the  forces,  and  the 
ring  will  be  destroyed. 

If  the  satellite,  instead  of  being  displaced  forwards,  had  been  originally  behind  its  mean 
position  in  the  ring,  the  forces  would  have  pulled  it  backwards,  its  path  would  have  become 
more  concave  towards  the  centre,  its  distance  from  the  centre  would  diminish,  its  angular 
velocity  would  increase,  and  it  would  gain  upon  the  rest  of  the  ring  till  it  got  in  front  of  its 
mean  position.  This  effect  is  of  course  dependent  on  the  very  same  conditions  as  in  the 
former  case,  and  the  actual  effect  on  a  disturbed   satellite    would    be   to   make    it  describe  an 
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orbit  about  its  mean  position  in  the  ring,  so  that  if  in  advance  of  its  mean  position,  it  first 
recedes  from  the  centre,  then  falls  behind  its  mean  position  in  the  ring,  then  approaches 
the  centre  within  the  mean  distance,  then  advances  beyond  its  mean  position,  and,  lastly, 
recedes  from  the  centre  till  it  reaches  its  starting-point,  after  which  the  process  is  repeated 
indefinitely,  the  orbit  being  always  described  in  the  direction  opposite  to  that  of  the  revolution 
of  the  ring. 

We  now  understand  what  would  happen  to  a  disturbed  satellite,  if  all  the  others  were 
preserved  from  disturbance.  But,  since  all  the  satellites  are  equally  free,  the  motion  of  one 
will  produce  changes  in  the  forces  acting  on  the  rest,  and  this  will  set  them  in  motion,  and 
this  motion  will  be  propagated  from  one  satellite  to  another  round  the  ring.  Now  propagated 
disturbances  constitute  waves,  and  all  waves,  however  complicated,  may  be  reduced  to  combi- 
nations of  simple  and  regular  waves;  and  therefore  all  the  disturbances  of  the  ring  may  be  con- 
sidered as  the  resultant  of  many  series  of  waves,  of  different  lengths,  and  travelling  with  different 
velocities.  The  investigation  of  the  relation  between  the  length  and  velocity  ,of  these  waves  forms 
the  essential  part  of  the  problem,  after  which  we  have  only  to  split  up  the  original  disturbance 
into  its  simple  elements,  to  calculate  the  effect  of  each  of  these  separately,  and  then  to  combine 
the  results.  The  solution  thus  obtained  will  be  perfectly  general,  and  quite  independent  of 
the  particular  form  of  the  ring,  whether  regular  or  irregular  at  starting.      §  14. 

We  next  investigated  the  effect  upon  the  ring  of  an  external  disturbing  force.  Having 
split  up  the  disturbing  force  into  components  of  the  same  type  with  the  waves  of  the  ring 
(an  operation  which  is  always  possible),  we  found  that  each  term  of  the  disturbing  force 
generates  a  "forced  wave"  travelling  with  its  own  angular  velocity.  The  magnitude  of  the 
forced  wave  depends  not  only  on  that  of  the  disturbing  force,  but  on  the  angular  velocity 
with  which  the  disturbance  travels  round  the  ring,  being  greater  in  proportion  as  this  velocity 
more  nearly  coincides  with  that  of  one  of  the  "  free  waves  "  of  the  ring.  We  also  found  that 
the  displacement  of  the  satellites  was  sometimes  in  the  direction  of  the  disturbing  force,  and 
sometimes  in  the  opposite  direction,  according  to  the  relative  position  of  the  forced  wave  among 
the  four  natural  ones,  producing  in  the  one  case  positive,  and  in  the  other  negative  forced  waves. 
In  treating  the  problem  generally,  we  must  determine  the  forced  waves  belonging  to  every 
term  of  the  disturbing  force,  and  combine  these  with  such  a  system  of  free  waves  as  shall 
reproduce  the  initial  state  of  the  ring.  The  subsequent  motion  of  the  ring  is  that  which 
would  result  from  the  free  waves  and  forced  waves  together.  The  most  important  class  of 
forced  waves  are  those  which  are  produced  by  waves  in  neighbouring  rings.     §  15. 

We  concluded  the  theory  of  a  ring  of  satellites  by  tracing  the  process  by  which  the  ring 
would  be  destroyed  if  the  conditions  of  stability  were  not  fulfilled.  We  found  two  cases  of 
instability,  depending  on  the  nature  of  the  tangential  force  due  to  tangential  displacement.  If 
this  force  be  in  the  direction  opposite  to  the  displacement,  that  is,  if  the  parts  of  the  ring  are 
statically  stable,  the  ring  will  be  destroyed,  the  irregularities  becoming  larger  and  larger 
without  being  propagated  round  the  ring.  When  the  tangential  force  is  in  the  direction  of  the 
tangential  displacement,  if  it  is  below  a  certain  value,  the  disturbances  will  be  propagated 
round  the  ring  without   becoming  larger,  and  we  have  the  case  of  stability  treated  of  at  large. 
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If  the   force   exceed   this    value,  the  disturbances   will  still    travel   round   the   ring,   but  they 
will  increase  in  amplitude  continually  till  the  ring  falls  into  confusion.      $  18. 

We  then  proceeded  to  extend  our  method  to  the  case  of  rings  of  different  constitutions. 
The  first  case  was  that  of  a  ring  of  satellites  of  unequal  size.  If  the  central  body  be  of 
sufficient  mass,  such  a  ring  will  be  spaced  out,  so  that  the  larger  satellites  will  be  at  wider  intervals 
than  the  smaller  ones,  and  the  waves  of  disturbance  will  be  propagated  as  before,  except  that 
there  may  be  reflected  waves  when  a  wave  reaches  a  part  of  the  ring  where  there  is  a  change 
in  the  average  size  of  the  satellites.      §  19. 

The  next  case  was  that  of  an  annular  cloud  of  meteoric  stones,  revolving  uniformly  about 
the  planet.  The  average  density  of  the  space  through  which  these  small  bodies  are  scattered 
will  vary  with  every  irregularity  of  the  motion,  and  this  variation  of  density  will  produce 
variations  in  the  forces  acting  upon  the  other  parts  of  the  cloud,  and  so  disturbances  will  be 
propagated  in  this  ring,  as  in  a  ring  of  a  finite  number  of  satellites.  The  condition  that  such  a 
rino-  should  be  free  from  destructive  oscillations  is,  that  the  density  of  the  planet  should  be 
more  than  three  hundred  times  that  of  the  ring.  This  would  make  the  ring  much  rarer  than 
common  air,  as  regards  its  average  density,  though  the  density  of  the  particles  of  which  it  is 
composed  may  be  great.  Comparing  this  result  with  Laplace's  minimum  density  of  a  ring 
revolving  as  a  whole,  we  find  that  such  a  ring  cannot  revolve  as  a  whole,  but  that  the  inner 
parts  must  have  a  greater  angular  velocity  than  the  outer  parts.      §  20. 

We  next  took  up  the  case  of  a  flattened  ring,  composed  of  incompressible  fluid,  and  moving 
with  uniform  angular  velocity.  The  internal  forces  here  arise  partly  from  attraction  and 
partly  from  fluid  pressure.  We  began  by  taking  the  case  of  an  infinite  stratum  of  fluid 
affected  by  regular  waves,  and  found  the  accurate  values  of  the  forces  in  this  case.  For  long 
waves  the  resultant  force  is  in  the  same  direction  as  the  displacement,  reaching  a  maximum 
for  waves  whose  length  is  about  ten  times  the  thickness  of  the  str.atum.  For  waves  about  five 
times  as  long  as  the  stratum  is  thick  there  is  no  resultant  force,  and  for  shorter  waves  the 
force  is  in  the  opposite  direction  to  the  displacement.      §  23. 

Applying  these  results  to  the  case  of  the  ring,  we  find  that  it  will  be  destroyed  by  the 
long  waves  unless  the  fluid  is  less  than  J?  of  the  density  of  the  planet,  and  that  in  all  cases 
the  short  waves  will  break  up  the  ring  into  small  satellites. 

Passing  to  the  case  of  narrow  rings,  we  should  find  a  somewhat  larger  maximum  density, 
but  we  should  still  find  that  very  short  waves  produce  forces  in  the  direction  opposite  to  the 
displacement,  and  that  therefore,  as  already  explained  (page  36),  these  short  undulations  would 
increase  in  magnitude  without  being  propagated  along  the  ring,  till  they  had  broken  up  the 
fluid  filament  into  drops.  These  drops  may  or  may  not  fulfil  the  condition  formerly  given 
for  the  stability  of  a  ring  of  equal  satellites.  If  they  fulfil  it,  they  will  move  as  a  permanent 
ring.  If  they  do  not,  short  waves  will  arise  and  be  propagated  among  the  satellites,  with  ever 
increasing  magnitude,  till  a  sufficient  number  of  drops  have  been  brought  into  collision,  so  as 
to  unite  and  form  a  smaller  number  of  larger  drops,  which  may  be  capable  of  revolving  as  a 
permanent  ring. 
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We  have  already  investigated  the  disturbances  produced  by  an  external  force  independent 
of  the  ring;  but  the  special  case  of  the  mutual  perturbations  of  two  concentric  rings  is 
considerably  more  complex,  because  the  existence  of  a  double  system  of  waves  changes  the 
character  of  both,  and  the  waves  produced  react  on  those  that  produced  them. 

We  determined  the  attraction  of  a  ring  upon  a  particle  of  a  concentric  ring,  first,  when 
both  rings  are  in  their  undisturbed  state ;  secondly,  when  the  particle  is  disturbed  ;  and, 
thirdly,  when  the  attracting  ring  is  disturbed  by  a  series  of  waves.     $26. 

We  then  formed  the  equations  of  motion  of  one  of  the  rings,  taking  in  the  disturbing 
forces  arising  from  the  existence  of  a  wave  in  the  other  ring,  and  found  the  small  variation  of 
the  velocity  of  a  wave  in  the  first  ring  as  dependent  on  the  magnitude  of  the  wave  in  the  second 
ring,  which  travels  with  it.      §  27- 

The  forced  wave  in  the  second  ring  must  have  the  same  absolute  angular  velocity  as  the 
free  wave  of  the  first  which  produces  it,  but  this  velocity  of  the  free  wave  is  slightly  altered  by 
the  reaction  of  the  forced  wave  upon  it.  We  find  that  if  a  free  wave  of  the  first  ring  has  an 
absolute  angular  velocity  not  very  different  from  that  of  a  free  wave  of  the  second  ring,  then 
if  both  free  waves  be  of  even  orders  (that  is,  of  the  second  or  fourth  varieties  of  waves),  or 
both  of  odd  orders  (that  is,  of  the  first  or  third),  then  the  swifter  of  the  two  free  waves  has  its 
velocity  increased  by  the  forced  wave  which  it  produces,  and  the  slower  free  wave  is  rendered 
still  slower  by  its  forced  wave ;  and  even  when  the  two  free  waves  have  the  same  angular 
velocity,  their  mutual  action  will  make  them  both  split  into  two,  one  wave  in  each  ring 
travelling  faster,  and  the  other  wave  in  each  ring  travelling  slower,  than  the  rate  with  which 
they  would  move  if  they  had  not  acted  on  each  other. 

But  if  one  of  the  free  waves  be  of  an  even  order  and  the  other  of  an  odd  order,  the 
swifter  free  wave  will  travel  slower,  and  the  slower  free  wave  will  travel  swifter,  on  account 
of  the  reaction  of  their  respective  forced  waves.  If  the  two  free  waves  have  naturally  a 
certain  small  difference  of  velocities,  they  will  be  made  to  travel  together,  but  if  the 
difference  is  less  than  this,  they  will  again  split  into  two  pairs  of  waves,  one  pair  continually 
increasing  in  magnitude  without  limit,  and  the  other  continually  diminishing,  so  that  one 
of  the  waves  in  each  ring  will  increase  in  violence  till  it  has  thrown  the  ring  into  a  state  of 
coufusion. 

There  are  four  cases  in  which  this  may  happen.  The  first  wave  of  the  outer  ring  may 
conspire  with  the  second  or  the  fourth  of  the  inner  ring,  the  second  of  the  outer  with  the 
third  of  the  inner,  or  the  third  of  the  outer  with  the  fourth  of  the  inner.  That  two 
rings  may  revolve  permanently,  their  distances  must  be  arranged  so  that  none  of  these 
conspiracies  may  arise  between  odd  and  even  waves,  whatever  be  the  value  of  m.  The  number 
of  conditions  to  be  fulfilled  is  therefore  very  great,  especially  when  the  rings  are  near  together 
and  have  nearly  the  same  angular  velocity,  because  then  there  are  a  greater  number  of 
dangerous  values  of  m  to  be  provided  for. 

In  the  case  of  a  large  number  of  concentric  rings,  the  stability  of  each  pair  must  be 
investigated  separately,  and  if  in  the  case  of  any  two,  whether  consecutive  rings   or   not,  there 
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are  a  pair  of  conspiring  waves,  those  two  rings  will  be  agitated  more  and  more,  till  waves 
of  that  kind  are  rendered  impossible  by  the  breaking  up  of  those  rings  into  some  different 
arrangement.  The  presence  of  the  other  rings  cannot  prevent  the  mutual  destruction  of  any 
pair  which  bear  such  relations  to  each  other. 

It  appears,  therefore,  that  in  a  system  of  many  concentric  rings  there  will  be  continually 
new  cases  of  mutual  interference  between  different  pairs  of  rings.  The  forces  which  excite 
these  disturbances  being  very  small,  they  will  be  slow  of  growth,  and  it  is  possible  that  by 
the  irregularities  of  each  of  the  rings  the  waves  may  be  so  broken  and  confused  (see  $  10),  as 
to  be  incapable  of  mounting  up  to  the  height  at  which  they  would  begin  to  destroy  the 
arrangement  of  the  ring.  In  this  way  it  may  be  conceived  to  be  possible  that  the  gradual 
disarrangement  of  the  system  may  be  retarded  or  indefinitely  postponed. 

But  supposing  that  these  waves  mount  up  so  as  to  produce  collisions  among  the  particles, 
then  we  may  deduce  the  result  upon  the  system  from  general  dynamical  principles.  There 
will  be  a  tendency  among  the  exterior  rings  to  remove  further  from  the  planet,  and  among 
the  interior  rings  to  approach  the  planet,  and  this  either  by  the  extreme  interior  and  exterior 
rings  diverging  from  each  other,  or  by  intermediate  parts  of  the  system  moving  away  from 
the  mean  ring.  If  the  interior  rings  are  observed  to  approach  the  planet,  while  it  is  known 
that  none  of  the  other  rings  have  expanded,  then  the  cause  of  the  change  cannot  be  the  mutual 
action  of  the  parts  of  the  system,  but  the  resistance  of  some  medium  in  which  the  rings 
revolve.       $31. 

There  is  another  cause  which  would  gradually  act  upon  a  broad  fluid  ring  of  which  the 
parts  revolve  each  with  the  angular  velocity  due  to  its  distance  from  the  planet,  namely,  the 
internal  friction  produced  by  the  slipping  of  the  concentric  rings  with  different  angular 
velocities.  It  appears,  however  (fi  33),  that  the  effect  of  fluid  friction  would  be  insensible  if 
the  motion  were  regular. 

Let  us  now  gather  together  the  conclusions  we  have  been  able  to  draw  from  the 
mathematical  theory  of  various  kinds  of  conceivable  rings. 

We  found  that  the  stability  of  the  motion  of  a  solid  ring  depended  on  so  delicate  an 
adjustment,  and  at  the  same  time  so  unsymmetrical  a  distribution  of  mass,  that  even  if  the 
exact  condition  were  fulfilled,  it  could  scarcely  last  long,  and  if  it  did,  the  immense  pre- 
ponderance of  one  side  of  the  ring  would  be  easily  observed,  contrary  to  experience.  These 
considerations,  with  others  derived  from  the  mechanical  structure  of  so  vast  a  body,  compel  us 
to  abandon  any  theory  of  solid  rings. 

We  next  examined  the  motion  of  a  ring  of  equal  satellites,  and  found  that  if  the  mass  of 
the  planet  is  sufficient,  any  disturbances  produced  in  the  arrangement  of  the  ring  will  be 
propagated  round  it  in  the  form  of  waves,  and  will  not  introduce  dangerous  confusion.  If  the 
satellites  are  unequal,  the  propagation  of  the  waves  will  no  longer  be  regular,  but  disturbances 
of  the  ring  will  in  this,  as  in  the  former  case,  produce  only  waves,  and  not  growing  confusion. 
Supposing  the  ring  to  consist,  not  of  a  single  row  of  large  satellites,  but  of  a  cloud  of  evenly 
distributed  unconnected  particles,  we  found  that  such   a  cloud  must  have  a  very  small  density 
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in  order  to  be  permanent,  and  that  this  is  inconsistent  with  its  outer  and  inner  parts  moving 
with  the  same  angular  velocity.  Supposing  the  ring  to  be  fluid  and  continuous,  we  found 
that  it  will  be  necessarily  broken  up  into  small  portions. 

We  conclude,  therefore,  that  the  rings  must  consist  of  disconnected  particles ;  these  may 
be  either  solid  or  liquid,  but  they  must  be  independent.  The  entire  system  of  rings  must 
therefore  consist  either  of  a  series  of  many  concentric  rings,  each  moving  with  its  own  velocity, 
and  having  its  own  systems  of  waves,  or  else  of  a  confused  multitude  of  revolving  particles,  not 
arranged  in  rings,  and  continually  coming  into  collision  with  each  other. 

Taking  the  first  case,  we  found  that  in  an  indefinite  number  of  possible  cases  the  mutual 
perturbations  of  two  rings,  stable  in  themselves,  might  mount  up  in  time  to  a  destructive 
magnitude,  and  that  such  cases  must  continually  occur  in  an  extensive  system  like  that  of 
Saturn,  the  only  retarding  cause  being  the  possible  irregularity  of  the  rings. 

The  result  of  long-continued  disturbance  was  found  to  be  the  spreading  out  of  the  rings 
in  breadth,  the  outer  rings  pressing  outwards,  while  the  inner  rings  press  inwards. 

The  final  result,  therefore,  of  the  mechanical  theory  is,  that  the  only  system  of  rings  which 
can  exist  is  one  composed  of  an  indefinite  number  of  unconnected  particles,  revolving  round 
the  planet  with  different  velocities  according  to  their  respective  distances.  These  particles  may 
be  arranged  in  series  of  narrow  rings,  or  they  may  move  through  each  other  irregularly.  In 
the  first  case  the  destruction  of  the  system  will  be  very  slow,  in  the  second  case  it  will  be  more 
rapid,  but  there  may  be  a  tendency  towards  an  arrangement  in  narrow  rings,  which  may  retard 
the  process. 

We  are  not  able  to  ascertain  by  observation  the  constitution  of  the  two  outer  divisions  of 
the  system  of  rings,  but  the  inner  ring  is  certainly  transparent,  for  the  limb  of  Saturn  has 
been  observed  through  it.  It  is  also  certain,  that  though  the  space  occupied  by  the  ring  is 
transparent,  it  is  not  through  the  material  parts  of  it  that  Saturn  was  seen,  for  his  limb  was 
observed  without  distortion  ;  which  shows  that  there  was  no  refraction,  and  therefore  that  the 
rays  did  not  pass  through  a  medium  at  all,  but  between  the  solid  or  liquid  particles  of  which 
the  ring  is  composed.  Here  then  we  have  an  optical  argument  in  favour  of  the  theory  of 
independent  particles  as  the  material  of  the  rings.  The  two  outer  rings  may  be  of  the  same 
nature,  but  not  so  exceedingly  rare  that  a  ray  of  light  can  pass  through  their  whole  thickness 
without  encountering  one  of  the  particles. 

Finally,  the  two  outer  rings  have  been  observed  for  200  years,  and  it  appears,  from  the 
careful  analysis  of  all  the  observations  by  M.  Struve,  that  the  second  ring  is  broader  than  when 
first  observed,  and  that  its  inner  edge  is  nearer  the  planet  than  formerly.  The  inner  ring  also 
is  suspected  to  be  approaching  the  planet  ever  since  its  discovery  in  1850.  These  appearances 
seem  to  indicate  the  same  slow  progress  of  the  rings  towards  separation  which  we  found  to  be 
the  result  of  theory,  and  the  remark,  that  the  inner  edge  of  the  inner  ring  is  most  distinct, 
seems  to  indicate  that  the  approach  towards  the  planet  is  less  rapid  near  the  edge,  as  we  had 
reason  to  conjecture.  As  to  the  apparent  unchangeableness  of  the  exterior  diameter  of  the 
outer  ring,  we  must  remember  that  the  outer  rings  are  certainly  far  more  dense  than  the  inner 
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one,  and  that  a  small  change  in  the  outer  rings  must  balance  a  great  change  in  the  inner  one. 
It  is  possible,  however,  that  some  of  the  observed  changes  may  be  due  to  the  existence  of  a 
resisting  medium.  If  the  changes  already  suspected  should  be  confirmed  by  repeated  observa- 
tions with  the  same  instruments,  it  will  be  worth  while  to  investigate  more  carefully  whether 
Saturn's  Rings  are  permanent  or  transitionary  elements  of  the  Solar  System,  and  whether  in 
that  part  of  the  heavens  we  see  celestial  immutability,  or  terrestrial  corruption  and  generation, 
and  the  old  order  giving  place  to  new  before  our  own  eyes. 


APPENDIX. 

On  the  Stability  of  the  Steady  Motion  of  a  Rigid  Body  about  a  Fixed  Centre 
of  Force.    By  Professor  W.  Thomson,  (communicated  in  a  letter). 

The  body  will  be  supposed  to  be  symmetrical  on  the  two  sides  of  a  certain  plane  containing 
the  centre  of  force,  and  no  motion  except  that  of  parts  of  the  body  parallel  to  the  plane  will  be 
considered.  Taking  it  as  the  plane  of  construction,  let  G  (fig.  14)  be  the  centre  of  gravity  of 
the  body,  and  O  a  point  at  which  the  resultant  attraction  of  the  body  is  in  the  line  OG  towards 
G.     Then  if  the  body  be  placed  with  O  coinciding  with  the  centre  of  force,  and  set  in  a  state 

of  rotation  about  that  point  as  an  axis,  with  an  angular  velocity  equal  to    \/  - — ,   (where  f 

denotes  the  attraction  of  the  body  on  a  unit  of  matter  at  0,  S  the  amount  of  matter  in  the 
central  body,  M  the  mass  of  the  revolving  body,  and  a  the  distance  OG),  it  will  continue, 
provided  it  be  perfectly  undisturbed,  to  revolve  uniformly  at  this  rate,  and  the  attraction  Sf 
on  the  moving  body  will  be  constantly  balanced  by  the  centrifugal  force  w2aM  of  its  motion. 
Let  us  now  suppose  the  motion  to  be  slightly  disturbed,  and  let  it  be  required  to  investigate 
the  consequences.  Let  X,  S,  Y,  be  rectangular  axes  of  reference  revolving  uniformly  with 
the  angular  velocity  a;,  round  S,  the  fixed  attracting  point.  Let  x,  y,  be  the  coordinates  of  G 
with  reference  to  these  axes,  and  let  XS,  YS  denote  the  components  of  the  whole  force  of 
attraction  of  S  on  the  rigid  body.  Then  since  this  force  is  in  the  line  through  S,  its  moment 
round  G  is 

SYx-  SXy; 
the  components  of  the  forces  on  the   moving  body  being  reckoned  as  positive  when  they  tend 
to  diminish  x  and  y  respectively.      Hence  if  k  denote  the  radius  of  gyration  of  the  body  round 
G,  and  if  <p  denote  the  angle  which  OG  makes  with  SX  (i.  e.  the  angle  GOK),  the  equations 
of  motion  are, 

Mis?  — i-  -  S  (  Yx  -  Xy )  =  0. 
dt! 

In  the  first  place  we  see  that  one  integral  of  these  equations  is 

M(xd^-yd*)+M»&  +  -f)+Mk>df  =  H. 
V    dt      J  dtl  y  dt 

This  is  the  "equation  of  angular  momentum." 
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In  considering  whether  the  motion  round  S  with  velocity  w  when  O  coincides  with  S  is 
stable  or  unstable,  we  must  find  whether  every  possible  motion  with  the  same  "angular 
momentum"  round  S  is  such  that  it  will  never  bring  O  to  more  than  an  infinitely  small  distance 
from  S:  that  is  to  say,  we  must  find  whether,  for  every  possible  solution  in  which  H  = 
M  (ffl*  +  A2)  w,  and  for  which  the  co-ordinates  of  O  are  infinitely  small  at  one  time,  these 
co-ordinates  remain  infinitely  small.  Let  these  values  at  time  t  be  denoted  thus:  SN  =  £,  and 
NO  =  rJ.  let  OG  be  at  first  infinitely  nearly  parallel  to  OX,  i.  e.  let  cp  be  infinitely  small 
(the  full  solution  will  tell  us  whether  or  not  <p  remains  infinitely  small);  then,  as  long  as  cp  is 
infinitely  small,  we  have 

x  =  a  +  £,  y  =  17  +  acp, 
and  the  equations  of  motion  have  the  forms 

M{\  ~  +  a -r?  +  2w  -5  -  w°  (ij  +  ad>)\  +  SY  =  0, 
[dt2  dt  dt  r  J 

and  we  may  write  the  equation  of  angular  momentum  instead  of  the  third  equation, 

AtUa  +  I)  (J  +  a  d-£)  -  (r,  +  «<£)  g)  +  a  (a  +  £)2  +  "  ('>  +  °<PY  +  k'  ^}  =  H 

If  now  we  suppose  £  and  ^  to  be  infinitely  small,  the  last  of  these  equations  becomes 

,  „      ,  „.  dd>  dij 

(a~  +  k)  -^  +  2wat  +  a  —  =  0 («)■ 

v  lit  s  dt 

If  p  and  q  denote  the  components  parallel  and  perpendicular  to  OG  of  the  attraction  of 
the  body  on  a  unit  of  matter  at  S,  we  have 

X  =  p  cos  cj)  -  q  sin  (p  =  p,   and  Y  =  p  sin  <p  +  q  cos  cp  =  pep  +  q, 
since  q  and  cp  are  each  infinitely  small;  and  if  we  put  V  =  potential  at  S,  and 
_  d-V         _  d'V  dxV 

°  =  dj s  '   P  ~  d^  '    7  =  rfp, ' 

then  p  =/-  a£  -  7>/,  q  =  -  fin  -  y%, 

^  =/-  «f  -  7.7,  F  =  /<£  -  /3«,  -  7?- 
If  we  make  these  substitutions  for  X  and  Y,  and  take  into  account  that 

./W«f (6), 

the  first  and  second  equations  of  motion  become 

£r*-ji-'*l-*"t-M{/**7*~* (r)' 

d*n  dE  (Feb      S    „ 

df+~wdt-^"  +  a^F-M^"  +  ^)  =  0   •    •    •        •    {d)- 
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Combining  equations  (a),  (c),  and  (d),  by  the  same  method  as  that  adopted  in  the  text, 
ire  find  that  the  differential  equation  in  £,  rj,  or  <p,  is  of  the  form 


,  d'u       _  d~u 
J  -T  +  B  — -  +  Cu  =  0, 
dtl  dtl 


where     ^  =  &% 


5  =  w2  (ftlf  -  (?)  -  -f  {k-a  +  («2  +  t*)/3|, 

C  =  wi  (F  -  3a*)  +  w»  -  «a»  +  *2)  (a  +  j3)  -  4as/3*  +  (a2  +  4s)  —  (a/3  -  7). 

In  comparing  this  result  with  that  obtained  in  the  Essay,  we  must  put 

rg  for  a, 

R  for  .af, 

R  +  S  for  & 

L   for  a, 

iVr02  for  j8, 

ilfr0  for  7. 
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